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Preface 


This  thesis  describes  the  design  and  implementation  of  a 
microprocessor  based  data  acquisition  and  processing  system  for  the  Air 
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be  performed. 
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Abstract 

In  order  to  design  a  constrained  layer  damping  treatment  for 
aircraft  components  the  Air  Force  Flight  Dynamics  Laboratory  at  Uright- 
Patterson  Air  Force  Base  Ohio  requires  a  knowledge  of  vibration  and 
tempet  iture  variations  which  a  component  encounters  during  its  life 
This  report  discusses  the  design  and  implementation  of  a  microprocessor 
based  system  to  acquire  vibration  and  temperature  information,  and 
computes  the  cumulative  damage  factor  versus  frequency  and  temperature 
in  real  tine.  The  cumulative  damage  ‘'actor  is  computed  for  one  of 
eight  user  selectable  frequency  ranges  and  for  a  temperature  range  of 
-50  to  +225  degrees  Fahrenheit  and  is  stored  in  a  non-volatile  memory- 
The  designed  system  uses  an  LSI  11  microcomputer  system  consisting 
of  an  LSI-11  processor,  an  analog-to-digital  converter,  a  real-time 
clock.  PROM,  a  non-volatile  memory,  and  a  serial  interface.  The  software 
developed  for  this  system  controls  the  data  acquisition,  the  computa¬ 
tion  of  a  fast  Fourier  transform  (FFT),  the  power  spectrum  of  vibration 
data,  and  the  subsequent  computation  of  the  cumulative  damage  ‘'actor. 

The  data  acquisition  and  FFT  modules  were  coded  in  assembly  language, 
and  the  other  modules  were  coded  in  FORTRAN-  The  output  of  the  system 
consists  of  a  table  containing  the  cumulative  damage  as  a  function  of 
frequency  and  temperature.  A  bench  nodel  system  was  constructed  an<! 


tested,  and  a  proposal  for  a  flight-worthy  system  is  made 
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DESIfiN  MS  IMPLEMENTATION  Q£  A. 
MICROPROCESSOR  FLIGHT  DA’IAOE  MONITOR 


I.  Introduction 

Aircraft  and  avionics  hardware  operate  in  an  environment  of 
varying  temperature  and  vibration.  The  resonant  fatigue  damage  caused 
by  this  environment  is  a  major  factor  in  airframe  and  avionics 
hardware  structural  failure.  In  order  to  reduce  the  resonant  fatigue 
damage,  design  engineers  must  know  what  temperatures  and  vibration 
frequencies  an  object  will  encounter  during  its  lice.  With  this 
information,  a  constrained  layer  damping  treatment  which  reduces  the 
resonant  fatigue  damage  can  be  applied  to  the  object  (Ref  10). 

The  Air  Force  Flight  Dynamics  Laboratory  (AFFDL),  Wright-Patterson 
Air  Force  Base  Ohio,  has  the  technology  to  produce  a  constrained 
layer  damping  treatment,  although  the  current  flight  data  acquisition 
and  analysis  systems  are  very  large  and  hard  to  use.  and  require  undue 
time  and  resources  (Ref  10).  In  order  to  simplify  the  data 
acquisition,  a  microprocessor  based  data  acquisition  and  processing 
system  to  acquire  the  necessary  data,  and  process  it  in  real-time,  was 
designed  (Ref  7).  This  investigation,  sponsored  by  AFFDL,  is  a 
continuation  of  that  effort  and  will  focus  on  an  expanded  design,  bench 
model  construction,  testing,  and  recommendations  for  a  flight-worthy 
system. 
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Damage  Factor  (DF) 


Background 

Research  has  determined  that  the  resonant  vibrations  which  cause 
structural  damage  occur  at  frequencies  under  6,000  Hz  and  in  a 
temperature  range  between  -50  degrees  Fahrenheit  and  +225  degrees 
Fahrenheit  (Ref  10).  The  damage  factor  (DF)  at  a  given  frequency  (f) 
and  temperature  (T)  can  be  approximated  with  the  following  equation 
(Ref  10): 

DFT(f)  -  K(A)B/(f)2B“1  ,  (1) 

where  k  ia  a  constant,  A  is  the  root  mean  square  (RMS)  value  of  the 
acceleration  at  the  specified  frequency  (f),  and  B  is  approximately 
3.3.  Figure  1  is  a  graph  of  damage  factor  vs.  frequency  computed  by 
using  Eq  1.  This  damage  factor,  when  computed  for  a  given  frequency 
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and  temperature,  and  summed  for  a  predetermined  period  of  time  (t), 
forms  the  cumulative  damage  factor  (CDF)  as  shown  in  Eq  2: 

CDF(f,T)  -  £DF  (f,T)  ,  (2) 

t 

The  cumulative  damage  factor  for  a  given  frequency  and  temperature 
forms  one  damage  coefficient  in  the  "three  dimensional"  cumulative 
damage  matrix,  as  shown  in  Figure  2  (Ref  10). 

Current  in-flight  data  acquisition  techniques  use  analog  signal 
conditioners  and  multi-channel  analog  tape  recorders  to  record  the 
frequency  and  temperature  data  necessary  to  compute  the  CDF  After  the 
data  has  been  acquired,  the  CDF  is  computed  and  plotted.  This 
technique  requires  skilled  personnel  to  operate  the  equipment,  and 
due  to  the  time  required  to  compute  the  cumulative  damage  factor,  there 
can  he  a  long  delay  between  the  ending  of  the  test  and  obtaining  the 
cumulative  damage  matrix  output  (Ref  10). 

In  order  to  simplify  this  process,  it  was  decided  to  design  a 
Microprocessor  Flight  Damage  Monitor  (MFDM)  which  would  acquire  and 
process  the  data  in  real-time.  Such  a  system  appeared  attractive,  as 
it  would  be  easy  to  install  and  remove,  would  not  require  skilled 
personnel  to  operate,  and  the  cumulative  damage  data  would  be  available 
immediately  after  the  test  (Ref  10). 

System  Concept 

The  MFDM  system  concept  consists  of  a  self-contained  configuration 
which  would  compute,  in  real  time,  the  cumulative  damage  factor  CDF 
from  analog  temperature  and  frequency  data,  and  store  the  cumulative 
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Fig  2.  Cumulative  Damage  Matrix 


damage  factor  in  non-volatile  memory .  This  concept  was  expanded  to 
include  an  interactive  capability  where  the  user  can  select  the 
frequency  range  and  resolution  for  which  the  CDF  is  to  be  computed, 
select  the  duration  of  the  test,  and  have  the  cumulative  danage  matrix 
printed  (Ref  10),  In  addition  to  the  interactive  requirements,  two 
8tart-up  modes  were  required.  The  first  mode  would  start  processing 
data  upon  system  power  initiation  and  the  second  node  would  start 
processing  data  when  given  a  start  command.  Figure  3  shows  an  activity 
diagram  (Ref  711)  which  summarises  the  system. 

Approach 

The  MFDM  (also  referred  to  as  the  "system"  )  was  developed  in  four 
phases  bassed  upon  exeepted  top  down  structured  approach.  The  four 
phases  were:  the  conceptual  phase  the  design  phase,  the  implementation 
phase,  and  the  testing  phase.  The  conceptual  phase  uses  structured 
analysis  and  design  to  define  the  detailed  requirements  placed  on  the 
system.  The  design  phase  included  the  algorithm  design  and  selection 
of  the  system  hardware.  During  the  implementation  phase,  the  system  is 
coded,  assembled,  and  given  a  functional  check  Then  the  testing  phase 
verifies  the  operation  of  the  system. 

Development  Environment 

An  LSI-11  microcomputer,  manufactured  by  Digital  Equipment 
Corporation  (DEC),  was  used  as  the  system  processor  because  of  its 
poverful  instruction  set,  processing  speed,  data  structure  and 
availability  of  support  hardware  and  software.  The  support  hardware 
used  in  the  MFDM  is  listed  in  Table  I.  In  addition  to  the  LSI-11 
microcomputer  and  support  hardware,  the  development  system  had  a 
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User  Inputs  1  J  Power  On/Off 


TABLE  I 

LSI-11  Support  Hardware 


Module 

Description 

1. 

KD11-J 

LSI-11  Processor 

2. 

MMV11-A 

4  K  by  16  bit  Core  Memory 

3. 

DLV11-J 

4  Channel  Asynchronous  Serial 
Line  Interface 

4. 

KWV11-A 

Programmable  Real-Time  Clock 

5. 

Anvil -A 

16  Channel  12  bit  Analog-to- 
Digltal  Converter 

6. 

MSV11-DB 

8  K  by  16  bit  Semiconductor 

Memory  (Used  Two  Modules) 

7. 

BCV1B 

Backplane  Expansion  Cables 

8. 

H9270 

Backplanes  (two) 

9. 

TEV11 

120  OHM  Bus  Terminator 

10. 

H780-H 

Powersupply 

Heathklt  H27  dual  floppy  disk  drive,  a  Heathkit  H14  line  printer,  and  a 
Lear  Siegler  ADM-3  terminal.  The  software  was  developed  using  DEC's  P.T- 
11  version  03-3  operating  system,  version  2.1-11  FORTRAN  IV  compiler, 
and  version  03.01  macroassembler  (Ref  8:643-647).  This  development 
system  was  located  atthe  Air  Force  Institute  of  Technology  Uright- 
Patterson  AFB,  Ohio. 

Scope  of  Thesis 

First  the  system  requirements  and  design  constraints  are  defined 
(Chap  II),  and  then  a  structured  analysis  and  design  of  the  system  is 
performed,  producing  a  structure  chart  for  the  entire  system  (Chap  III). 
Following  this,  the  hardware  selections  (Ref  7  25-36)  are  reviewed  with 
respect  to  the  new  system  design  (Chap  IV).  After  the  hardware 
considerations  the  required  software  modules  are  developed  by  modifying 
existing  nodules  (Ref  7; 75-95)  where  possible,  and  creating  new  modules 
for  functions  which  had  not  yet  been  implemented  (Chap  V).  The  hardware 
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and  the  software  were  then  Integrated  Into  a  bench  model  of  the  MFDM 
which  was  functionally  tested  (Chap  VI).  A  proposed  flight-worthy  system 
Is  then  discussed  (Chap  VII).  Recommendations  for  the  f light-worthy 
system  and  investigative  conclusions  about  this  effort  are  presented 
In  Chan  VIII. 


II.  Requirements  Definition 


( 

Introduction 

Before  any  system  can  be  designed  and  implemented,  the  system 
requirements  must  be  defined.  These  requirements  form  a  base  line  from 
which  a  system  can  be  designed,  and  must  include  the  purpose  of  the 
system,  its  input  parameters,  output  parameters,  any  time  constraints, 
and  any  hardware  constraints.  The  purpose  of  this  chapter  is  to 
accurately  define  the  system  in  terms  of  the  requirements  established 
by  the  AFFDL  (Ref  10)  and  the  design  constraints  of  the  system. 

System  Requirements 

The  purpose  of  the  MFDM  is  the  real-time  computation  of  the  CDF  as 
a  function  of  frequency  and  temperature  from  analog  frequency  and 
*  temperature  signals.  The  output  of  this  system  will  he  a  cumulative 

damage  table  which  contains  the  CDF  as  a  function  of  frequency  and 
temperature.  This  gives  a  general  definition  of  the  system,  the 
specific  requirements  of  which  are  as  follows  (Ref  7:6-22): 


-  The  input  to  the  system  will  he  an  analog  frequency 
signal  and  an  analog  temperature  signal.  Signal 
amplifiers  will  amplify  the  signals  so  as  to  make  the 
signals  compatible  with  an  analog-to-digital  converter. 

-  The  temperature  range  for  which  resonant  fatigue  damage 
usually  occurs  is  between  -50  degrees  and  +225  degrees 
Fahrenheit.  In  order  to  produce  a  damning  layer,  the 
cumulative  damage  factor  must  be  known  at  25  degree 
intervals  for  this  temperature  range. 

-  The  frequency  range  over  which  damage  occurs  depends  upon 
the  object  being  tested  as  well  as  its  environment.  For 
these  reasons,  the  frequency  range  and  interval  for  which 
the  DCF  will  be  computed  must  be  user  selectable.  Table 
II  gives  the  frequency  ranges/intervals  for  which  the  CDF 
must  be  computed  in  order  to  provide  the  user  with  the 
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capability  of  selecting  the  combination  which  beat  meets 
his  needs* 

This  system  must  be  capable  of  operating  for  the  duration 
of  a  transport  aircraft's  mission*  The  duration  of  these 
missions  can  be  as  long  as  20  hours,  so  the  MFDM  must  be 
capable  of  operating  for  20  hours,  because  the  length  of 
test  missions  vary  the  user  must  be  able  to  set  the 
desired  run  tine  before  the  mission. 

The  user  must  be  able  to  interactively,  fron  a  computer 
terminal,  configure  the  system  as  desired.  The  interac¬ 
tive  aspects  of  the  system  nust  bo  very  easy  to 
understand  since  they  may  be  used  by  people  without  a 
computer  background.  Thus,  the  user  should  be  pronpted 
for  responses  and  nust  be  able  to  respond  with  simple 
responses  such  as  one  letter  or  one  numeral.  If  a  wrong 
input  is  entered  by  the  user  an  appropriate  error  message 
will  be  written  at  the  terminal,  and  the  prompt  line 
repeated. 

The  user  must  be  able  to  enter  a  test  date  and  a  line 
of  comments  into  the  computer.  This  date  and  comnent  is 
to  be  printed  in  the  heading  of  the  CDF  output  table 
along  with  the  total  time  the  system  processed  data.  This 
will  provide  the  information  needed  for  correlating  the 
output  tables  with  the  test  mission. 

There  must  be  two  modes  of  operation:  one,  an  automatic 
node  where,  upon  application  of  system  power,  data 
processing  starts  automatically  for  the  last,  user  set, 
system  configuration;  and  a  second  manual  mode  hv  which 
the  U3er  can  start  data  processing  by  giving  a  start 
comand  from  a  terminal.  This  will  allow  the  system  to 
be  used  automatically  when  on  a  plane,  and  manually  when 
used  in  a  laboratory  - 

When  data  processing  has  been  performed  for  the  desired 
length  of  tine,  data  processing  will  automatically  stop. 
Program  control  will  return  to  the  user  if  operating  in 
the  manual  mode  and  will  halt  if  operating  in  the 
automatic  mode.  This  will  prevent  the  system  fron 
processing  data  after  a  mission  is  over,  and  will  allow 
the  user  to  control  the  length  of  time  for  which  data  is 
processed. 

When  processing  data,  the  CDF  matrix  nust  be  updated  as 
fast  as  possible.  CDF  update  rate  of  once  every  6 
seconds  is  desired  as  this  will  provide  at  least  five 
updates  during  a  plane's  3Q-3ccond  take-off  roll. 


10 


TABLE  II 


The  Required  Frequency  Ranges  and  Intervals 


Frequency  Range 
(0  -  MAXFREO) 

(K  Hz) 

Interval 

(Hz) 

0  -  6.0 

100 

0-30 

100 

0  -  3.0 

50 

0  -  3.0 

25 

0  -  1.5 

100 

0  -  1-5 

50 

0  -  1.5 

25 

-  The  dynamic  range  of  the  CDF  matrix  must  be  at  least  360 

db,  wl^ch  allows  any  two  elements  of  the  matrix  to  differ 
by  10  .  The  reason  for  the  large  dynamic  range  is  that  if 

there  is  a  large  cumulative  damage  factor  at  one 
frequency  its  large  magnitude  will  not  cause  smaller 
cumulative  danage  factors  to  be  lost  due  to  the  finite 
word  size  of  the  computer. 

-  The  system  program  and  CDF  matrix  must  be  stored  in  a  non¬ 
volatile  memory  so  program  execution  can  start 
automatically  upon  power-up,  and  system  power  can  be 
turned  off  after  the  processing  is  complete,  without  the 
loss  of  the  CDF  matrix. 

-  When  an  automatic  power-up  occurs  the  CDF  matrix  will  be 
cleared  only  if  it  is  the  first  power-up  since  a  manual 
power-up.  This  allows  data  to  be  acquired  from  multiple 
flights  if  a  manual  power-up  is  not  performed  between 
flights. 


This  completes  the  li3t  of  system  requirements  defined  for  this 
project.  Two  other  system  requirements  exist  but  are  not  requirements 
for  this  phase  of  the  project  (Ref  10).  They  are:  the  system  must  be 
battery-powered,  and  it  must  be  as  small  as  possible.  This  will  allow 
the  system  to  be  used  in  a  small  space  where  no  power  is  available 
such  as  in  the  wing  of  an  airplane.  Engineers  at  AFFDL  will  design  the 
case  for  the  system  as  well  as  the  battery  pack  required  to  power  the 


system  when  an  operational  bench  model  system  has  been  produced  and 

tested* 

Summary 

This  chapter  has  defined  the  general  system  requirements.  These 
requirements  will  now  be  used  to  design  algorithms  which  will  be  t3ed 
to  implement  the  systen.  The  design  and  implementation  of  the 
algorithm  is  the  subject  of  the  following  chapters 
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III.  Alp, or! thra  Design 


Introduction 

This  chapter  explains  the  design  of  the  algorithms  which  were  used 
to  implement  the  systen  according  to  the  requirements  defined  in 
Chapter  II.  These  algorithms  were  developed  using  a  combination  of 
structured  design  techniques  (Ref  9),  resulting  in  a 
structure  chart  which  describes  the  system.  The  structure  chart  and 
associated  module  interface  table  developed  in  this  chapter,  along  with 
the  discussion  of  the  function  of  each  module  which  nnke  up  the 
structure  chart  define  the  algorithm. 

Algorithm  Development 

The  algorithm  for  this  systen  was  developed  in  three  steps. 

First,  an  algorithm  was  developed  for  the  "data  processing"  function. 
Then  one  was  developed  for  the  interactive  user  "configuration" 
function,  and  finally  the  two  algorithms  were  integrated  into  one 
"system"  algorithm. 

Data  Processing  Algorithm 

The  basic  steps  which  must  be  performed  to  compute  the  cumulative 
damage  factor  are  as  follows  (Ref  10): 

Obtain  analog  vibration  and  temperature  signals  from 
the  item  being  tested  and  digitize  the  signals 
(subroutine  EXECUTIVE). 

Compute  the  Fourier  transform  of  the  vibrational  data 
for  the  frequency  range  and  interval  of  interest 
(subroutine  FFT). 

Compute  the  average  power  in  each  frequency  interval 
(subroutine  POUERSPECTRUM) . 
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Compute  the  damage  factor  for  each  frequency  interval 
(subroutine  DA.’iAOEFACTOR)  • 

Repeat  these  steps  as  long  as  desired,  summing  the  damage 
factor  into  a  cumulative  damage  matrix  which  contains 
the  cumulative  damage  versus  frequency  and  temperature. 

These  five  steps  are  illustrated  in  Figure  4.  The  linear 
amplifiers  shown  in  the  figure  are  used  to  amplify  the  voltage  signal 
of  the  sensors  to  a  level  required  by  the  analog-to-digital  converter. 

In  order  to  show  the  relationship  between  the  modules  discussed  in 
the  following  paragraphs,  a  dataflow  diagram  and  structure  chart  for 
the  process  just  described  is  shown  in  Figures  5  and  6  respectively 
These  two  figures  are  constructed  from  data  obtained  from  Ref  7  55-65 
which  discusses  a  basic  system  to  compute  the  cumulative  damage 
factor.  The  next  step  is  to  define  the  functions  of  each  of  the 
modules  shown  in  Figure  6.  The  Fourier  transform  nodule  (Ref  7  80- 
85),  will  be  discussed  first,  followed  by  the  power  spectrum  and  damage 
factor  modules.  After  completing  this,  the  nodule  interfaces  will  be 
discussed. 

FFT.  This  module  is  responsible  for  the  computation  of  the  fast 
Fourier  transform  (FFT)  of  the  digitized  vibration  data.  It  converts 
the  data  from  a  function  of  tine  to  a  function  of  frequency.  Since  the 
input  data  is  a  sampled  version  of  a  continuous  time  domain  signal,  a 
discrete  Fourier  transform  method  has  to  be  used.  Eq  3  gives  an 
expression  for  the  transform  of  N  sampled  points  from  the  tine  domain 
to  the  frequency  domain. 
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Frequency  Channel 


Computing  The  CDF 


Diagram  for  Data  Processing 


Fig  6.  Structure  Chart  for  Data  Processing 


Equation  3  la: 


F(n> 


k-N-1 

-  E 


k-0 


X(k)Ok 


n  -  0  1,2,  .  N-l  (3) 


where 


W  -  exp (-2  l/N)  (4) 

F(n)  Is  the  nth  point  In  tl  v  domain,  and  X(k)  Is  the  kth 

point  in  the  time  domain.  ...1:)  may  be  real  or  complex,  but  F(n)  is 

always  complex  (Ref  7:70).  Constraints  on  Eq  3  require  that  the 

frequency  function  be  band  limited  at  some  frequencv  f  .and  in  order 

c 

to  insure  that  no  aliasing  will  occur  in  the  frequency  domain,  the 
time  domain  function  oust  be  sampled  at  a  rate  at  least  twice  f^  .  The 
frequency  2f^  -1/T  is  known  as  the  Nyquist  sampling  rate,  where  T  Is 
the  sampling  interval  (Ref  3  85).  The  number  of  sampled  points  (V.) 
can  be  determined  using  Eq  5: 

NAf  -  1/T  -  2f  (5) 

c 

where  Af  is  the  desired  frequency  interval  in  the  frequency  domain  and 

l/T  Is  the  Nyquist  sampling  frequency  (Ref  3:85-96). 

Solving  Eq  3  with  a  computer  is  a  slow  and  time-consuming 

2 

process  as  it  requires  on  the  order  of  N  complex  multiplications  and 
additions  (Ref  4:316).  To  speed  up  this  calculation,  Cooley  and 
Tukey  developed  a  fast  Fourier  transform  for  N  being  an  integer  power 
of  two.  This  reduced  the  order  to  Nx  complex  multiplications  and 
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additions  whore  x  is  the  integer  power  of  two,  satisfying  N«2  (Ref 
3:174).  Since  X(k)  in  Eq  3  can  be  real  or  complex,  it  ia  possible 
to  solve  Eq  3  with  H  complex  data  points  and  have  the  H  complex 
data  points  consist  of  a  2N  point  real  data  sequence  (Ref  3:167). 

This  reduces  the  order  of  complexity  to  Nx/2.  Thus  Eq  3  can  be 
efficiently  implemented  on  a  computer  for  N-2  where  x  is  an  integer  and 
the  number  of  real  data  points  is  2M.  A  complete  explanation  of  this 
procedure  can  be  found  in  Ref  3  148-238  and  Ref  7  141-156.  Table  III 
shows  the  values  of  2M  and  l/T  which  satisfy  these  equations,  using  the 
requirements  for  the  frequency  range  and  interval  defined  in  Chapter 
II.  After  completing  the  Fourier  transform  of  the  input  data,  the 
power  within  each  frequency  interval  must  be  calculated. 

POWERS PECTRUM.  This  module  computes  the  power  in  each  of  the 
frequency  intervals  which  make  up  the  frequency  spectrum.  The 
frequency  component  F(n)  resulting  from  the  Fourier  transform  is  a 
complex  number  which  can  be  represented  as  given  in  Eq  6: 

F(n)  -  X2  +  iY2  ,n  -  0.1,2 _  N-l  (6) 

n  n 

and  the  power  within  a  frequency  interval  can  be  computed  using 
Eq  7: 

A(n)  -  (X2)  +  (Y2)  ,  n  -  0,1,2.  -  N-l  (7) 

n  n 

Because  the  input  signal  is  random,  Eq  7  represents  an 
Instantaneous  power,  and  that  which  is  needed  is  the  average 
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Table  III 

Nyquist  Sanpling  Rate  for  the  Required  Frequency  Ranges 


Frequency  Range 
(0  -  MAXFREQ  K  Hz) 

Interval 

(Hz) 

2N 

Sanpling  Rate 
(1/T  K  Hz) 

0  -  6.0 

100 

128 

12.8 

0  -  3-0 

100 

64 

6.4 

0  -  3.0 

50 

128 

6.4 

0  -  3  0 

25 

256 

6.4 

0  -  1.5 

100 

32 

3.2 

0  -  15 

50 

64 

3-2 

0  -  1.5 

25 

128 

3.2 

power  (Ref  10).  The  average  power  (A)  can  be  approximated  by  taking 
the  average  of  a  set  of  P  instantaneous  power  calculations  (Ref  7  73). 
Mathematically,  this  can  be  expressed  as: 

P-1 

A(n)  -  1/P 53  Ak(n)  ,  n  «  0,1.2  ...,N-1  (8) 

k«l 

where  the  subscript  k  represents  a  unique  set  of  N  instantaneous  power 
calculations  The  average  power  computed  using  Fq  8  can  be  interpreted 
as  the  average  (RMS)  value  of  the  acceleration  at  the  center  of  the 
frequency  interval,  which  is  the  frequency  of  F(n).  After  completing 
the  power  calculation,  the  damage  factor  can  be.  calculated. 

DAMAGRF ACTOR •  The  damage  factor  (DF)  can  be  computed  using 
Eq  1,  which  is  repeated  here: 

DF(f (n))  -  K(A(n))3'3/f (n)5'4  ,  n  -  1  2, 3, - N-l  (1) 

where  k  is  a  constant  for  all  frequency  f(n)  and  3  in  Fq  1  was 
replaced  with  its  value  of  3.3.  Since  the  damage  factor  is  a  relative 
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quantity  (Ref  10),  and  the  K  factor  Is  in  all  n  damage  factor 
computations  It  can  be  factored  out  as  In  F.q  9. 

DF(n)  oc  A3*3/f(n)5*4  n  -  1  2, 3, . . . .N-l  (9) 

The  cumulative  damage  factor  (CDF)  is  then  the  sum  of  the  damage 
factors  calculated  using  Sq  9,  for  a  period  of  time  T.  This  Is 
expressed  mathematically  as: 

CDF(f (n)  T)  -  £l>Ft(f(n))  f  n  -  1.2, 3 . N-l  (10) 

where  t  Is  the  rate  at  which  the  DF  (f  )  Is  to  be  calculated,  and  T  Is 
the  temperature  Interval  at  which  the  damage  occurred. 

This  completes  the  discussion  of  the  data  processing  nodules-  The 
remaining  module,  EXECUTIVE  is  the  module  which  controls  the  data 
processing. 

EXECUTIVE.  The  function  of  this  module  is  to  acquire  the  2N 
vibration  data  points  and  the  temperature  data,  and  coordinate  the 
computation  of  the  cumulative  damage  factor.  The  first  time  this 
nodule  is  called,  it  must  clear  the  cumulative  damage  matrix  before 
the  data  processing  can  begin.  In  addition,  it  must  monitor  the 
system  time,  and  stop  data  processing  when  the  users  set  run  time  has 
expired. 

Modules* Interfaces.  The  next  step  in  developing  the  algorithm  is 
to  establish  the  interfaces  between  the  nodules.  Table  IV  lists  the 


Table  IV 


Interfaces  for  Data  Processing  Modules 


INPUT 

OUTPUT 

ADDATA 

FFTDATA 

PWRDATA,  TEMP 

ADDATA,  TEMP 

FFTDATA 

PWRDATA 

CDFMATRIX 

input  and  output  Interfaces  corresponding  to  the  numbered  lines  in 
Figure  6.  The  executive  nodule  nust  acquire  the  2N  data  points 
(ADDATA)  and  pass  then  to  the  FFT  nodule  which  returns  the  transformed 
data  points  as  N  complex  frequency  points  (FFTDATA).  These  complex 
frequency  points  are  then  passed  to  the  power  spectrum  module  which 
computes  the  power  in  each  interval  of  the  frequency  spectrum.  The  N 
real  power  values  are  then  returned  in  the  array  HJRDATA.  The  damage 
factor  module  uses  PWRDATA  to  compute  the  damage  factor  which  is 
then  summed  into  the  cumulative  damage  factor  matrix  (CDFMATRIX) .  The 
cumulative  damage  matrix  is  the  final  product .  and  it  has  a  maximum 
size  of  12  N  elements,  where  the  12  is  the  number  of  temperature 
intervals  required  (Ref  Chap  II). 

This  completes  the  design  of  the  data  processing  algorithm 
required  to  satisfy  the  requirements  defined  in  Chapter  II.  The  next 
section  will  discuss  the  design  of  the  "configuration"  algorithm  which 
satisfies  the  interactive  requirements  defined  in  Chapter  II. 

Configuration  Algorithm 

A  "configuration"  algorithm  was  developed  which  meets  the 
interactive  requirements  defined  in  Chapter  II.  This  initializes  the 
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number  of  data  points  (N)  which  the  data  processing  algorithm 
requires. and  also  initializes  the  rate  at  which  the  analog  vibration 
signal  is  to  be  sampled*  The  requirements  defined  in  Chapter  II 
specify  that  the  user  be  able  to  interactively  do  the  following: 
configure  the  system  to  process  data  for  a  specified  frequency  rarge 
interval,  and  time  period;  initiate  data  processing;  and  have  the  CDF 
matrix  printed*  The  requirements  for  an  "interactive  system"  also 
specify  that  the  system  be  easy  to  operate  by  someone  without  any 
computer  experience,  that  the  the  user  be  prompted  for  a  response,  and 
that  the  user's  response  consist  of  a  single  letter  or  numeral. 

Because  of  these  requirements,  it  was  decided  to  design  an 
interactive  system  similar  to  that  used  in  the  UCSD  Pascal  operating 
system  (Ref  12).  This  interactive  system  was  chosen  because  the  author 
felt  it  was  the  best  systm  he  had  used  and  it  meets  the  requirements 
of  the  system.  The  user  is  prompted  with  the  allox/able  responses  in 
the  following  form:  the  first  letter  of  the  response,  a  left 
parenthesis,  and  the  remainder  of  the  word.  If  one  of  the  possible 
functions  was  to  setup  the  system  for  a  new  configuration,  one  of  the 
prompts  would  be  "S(ETUP".  If  this  was  the  desired  function,  the  user 
would  enter  an  "S"  followed  by  a  carriage  return  ("<CR>").  Thus,  by 
reading  the  prompt  lines  and  entering  the  first  letter  of  his  choice, 
the  user  can  control  which  function  is  executed.  If  an  illegal 
response  is  made,  an  error  message  is  printed  at  the  user's  terminal 
and  the  prompt  line  is  repeated,  allowing  the  user  to  enter  a  legal 
response.  A  structure  chart  of  the  algorithm  designed  to  meet  these 
requirements  is  shown  in  Figure  7,  and  a  discussion  of  the  functions  of 
each  of  the  modules  which  make  up  the  structure  chart  follows. 
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MFDM.  The  Microprocessor  Flight  Damage  Monitor  controlling 
nodule  is  the  MFDM*  This  is  the  module  which  is  executed  first  when 
the  system  is  started  in  the  "manual"  mode.  This  module  will  prompt 
the  user  to  input  the  function  which  he  desires  to  be  executed.  The 
user  will  respond  by  entering  the  fist  letter  of  the  desired  function 
-  either  S(ETUP,  E(XECUTE,  0(UTPUT,  or  Q(UIT.  Depending  on  the  input 
value,  one  of  the  subordinate  modules  SETUP,  XKCUTE,  or  CDFOUTPUT  will 
be  called,  unless  the  input  is  a  "Q"  which  will  stop  program  execution. 

If  an  illegal  entry  is  made,  the  error  module  will  be  called  to 
print  out  an  error  message  and  then  return  control  to  the  beginning  of 
the  nodule.  After  the  called  subordinate  module  completes  execution, 
program  control  will  pass  back  to  MFDM. 

SETUP.  The  functions  which  the  SETUP  module  must  perform  are 
defined  by  the  system  requirements.  The  user  must  be  able  to  select  a 
specific  frequency  range  and  frequency  interval,  enter  the  run  time  for 
which  data  is  to  be  processor,  and  enter  a  date  and  a  line  of  comments 
which  are  to  be  printed  in  the  header  of  the  CDF  table.  Tn  order  to 
make  this  algorithm  easier  to  modify  should  future  requirements 
warrant  it,  this  nodule  was  designed  to  be  an  interactive  module.  It 
first  prints  the  last  used  system  configuration  and  then  prompts  the 
user  to  enter  the  name  of  the  function  he  desires  to  change.  Then 
program  control  passes  to  that  module.  If  an  illegal  entry  is  made, 
the  error  module  will  be  called  to  print  out  an  error  message  and  then 
return  control  to  the  beginning  of  the  module.  By  using  this  method  of 
design,  future  modifications  should  only  require  the  adding,  deleting, 
or  changing  of  small  special  purpose  modules  to  make  changes  to  the 
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functions  performed  by  SETUP 


GETDATE  and  GETCOtlMENT.  These  two  nodules  read  the  data /comment 
from  the  terminal  as  the  user  inputs  it.  The  date  can  be  up  to  13 
characters  long,  and  the  conments  up  to  40  characters  long.  No  error 
checking  is  performed  since  whatever  is  input  will  be  printed  in  the 
CDF  output  table  heading  in  the  sane  format  as  it  was  input.  If  the 
input  is  longer  than  the  13  or  40  characters  allowed,  the  input  string 
will  be  truncated  to  the  first  13  or  40  characters. 

GETFREQUENCY.  If  the  frequency  range/interval  is  to  be  changed, 
this  module  is  called.  It  will  first  print  a  list  of  possible  frequency 
range/interval  combinations  as  shown  in  Table  II.  Then  it  will  prompt 
the  user  to  enter  the  combination  which  is  desired.  As  in  previous 
modules,  if  an  illegal  response  is  made  the  ERROR  module  will  be  called 
with  an  error  message.  After  control  returns  from  the  error  module, 
the  list  of  possible  combinations  will  be  printed  at  the  terminal 
followed  by  a  repeat  of  the  prompt  for  the  user's  response. 

When  a  valid  response  Is  received,  this  module  will  initialize  the 
global  variables  used  In  other  routines,  according  to  the  selected 
configuration.  The  variables  which  must  be  initialized  are:  the  number 
of  data  points  to  be  sampled  (N),  the  frequency  interval  (INTERVAL), 
the  number  of  frequency  intervals  (NFREQ),  and  the  data  sampling  rate 
(CLKRATE).  After  initializing  these  variables  control  is  returned 
to  the  SETUP  module  where  other  changes  can  be  requested  by  the  user. 

GETTIME.  This  module  procedure  will  be  used  if  the  duration  of 
the  test  is  to  be  changed.  It  will  accept  as  input  the  planned  length 
of  the  test,  in  hours,  and  check  that  the  tine  is  greater  than  zero. 

If  an  illegal  input  is  received,  procedure  ERROR  will  be  called  with  an 
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appropriate  error  message  followed  by  a  request  to  enter  a  new  time. 
After  receiving  a  valid  tine,  the  variable  "DURATION"  will  contain  the 
run  time  in  hours.  This  completes  the  discussion  of  the  modules  called 
by  SETUP  except  for  the  error  module  which  will  be  discussed  next 

ERROR.  This  module  will  be  cal’ed  by  any  module  which  receives  an 
Invalid  input  from  the  user.  The  error  message  to  be  printed  will  be 
passed,  as  an  argument,  to  this  module.  The  word  "ERROR"  will  be 
concatenated  to  the  beginning  of  the  error  message,  and  the  message 
will  be  printed  at  the  user's  terminal.  After  printing  the  message, 
program  control  returns  to  the  calling  module. 

XECUTE.  This  nodule  is  called  when  the  user  wants  to  start 
program  execution  fron  the  terminal.  This  module  will  check  the  run 
time,  and  if  it  is  positive,  a  message  indicating  the  start  of  data 
processing  will  be  printed  at  the  user's  terninal.  Program  control 
will  then  pass  to  the  EXECUTIVE  module  in  the  "process  data"  algorithm. 
If  the  run  tine  Is  less  than  or  equal  to  zero,  the  TIME  nodule  will  be 
called  so  that  the  user  can  enter  the  desired  run  time.  After 
completion  of  data  processing,  program  control  will  return  to  this 
module.  A  message  indicating  termination  of  data  processing  will  then 
be  printed  at  the  user's  terminal,  and  program  control  will  return  to 
MFDM. 

The  last  module  which  can  be  called  by  MFDM  is  the  CDFOUTPUT 
module  which  prints  out  the  CDF  matrix. 

CDFOUTPUT .  This  module  prints  the  CDF  matrix  at  a  printer 
connected  to  the  system.  The  CDF  will  be  printed  in  tabular  form  as  a 
function  of  frequency  and  temperature.  Each  page  of  the  table  will 


contain  a  heading  which  will  include  the  user-supplied  date  and 
comments  in  addition  to  the  time  in  which  data  was  actually  processed. 
After  printing  the  CDF  matrix,  program  control  will  return  to  MFDM. 

This  completes  the  development  of  the  "configuration"  algorithm, 
except  for  the  definition  of  the  interfaces  which  will  be  done  when  the 
two  algorithms  are  combined  into  one  system  algorithm. 

System  Algorithm 

The  system  algorithm  is  the  combination  of  the  two  algorithms  just 
developed.  A  structure  chart  for  the  system  algorithm  is  shown  in 
Figure  8.  The  module  interfaces  are  given  in  Table  VII,  which  is  shown 
in  Chapter  V.  The  only  requirement  that  the  two  separate  algorithms  do 
not  satisfy  is  that  of  the  "automatic  start-up"  mode.  The  PWRUPDN 
module  fulfills  this  requirement. 

PWRUPDN.  This  node  of  operation  requires  a  non-volatile  memory 
for  both  the  data  and  the  program.  When  system  power  is  applied, 
program  control  will  then  pass  to  the  beginning  of  the  EXECUTIVE  nodule 
where  the  CDF  array  is  cleared.  If  the  power  is  removed  and  then 
reapplied,  the  power-up  sequence  will  have  to  check  the  run  tine  to  see 
if  it  is  still  less  than  the  pre-set  run  time.  If  so,  program  control 
will  have  to  pass  to  the  EXECUTIVE  module  at  a  point  after  the  CDF 
array  has  been  cleared  so  data  processing  can  continue.  If  the  run  time 
has  expired,  the  program  execution  will  halt- 

Having  an  automatic  power-up  module  implies  that  something  should 
be  done  in  case  of  a  power  loss.  This  module  is  designed  so  that  all 
the  power-down  module  has  to  do  is  execute  a  halt  command.  Data 
processing  will  automatically  restart  when  power  is  restored. 
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Sumnary 

This  chapter  discussed  the  design  of  the  algorithm  which  when 
Implemented,  will  fulfill  the  requirements  established  in  Chapter  II. 
The  algorithm  was  designed  with  specallsed  modules,  each  having  one 
function.  This  design  method  allows  future  modifications  to  be  made  to 
a  module  without  having  to  modify  the  entire  system.  The  interactive 
aspect  was  designed  to  prompt  the  user  for  single  letter  or  numeral 
responses  thus  making  it  easier  to  operate  by  someone  inexperienced  in 
its  operation.  The  LSI-11  which  was  previously  selected  ns  the  "best" 
computer  to  use  for  this  effort  (Ref  7:25)  will  be  presented,  with  this 
algorithm  in  mind,  in  the  next  chapter. 
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IV.  Hardware  Considerations 


Introduction 

The  algorithm  discussed  in  Chapter  III  describes  the  system  in 
detail,  but  does  not  discuss  its  implementation.  This  chapter  presents 
the  decisions  that  were  made  concerning  the  implementation  of  the 
system.  The  decisions  made  previous  to  this  investigation  (Ref  7:2 3— 
35)  are  discussed  in  the  light  of  the  new  system  requirements. 

LSI  -1 1  Microcomputer  System 

The  LSI-1 1M  Microcomputer  manufactured  by  Norden  Systems,  Inc.  was 
previously  selected  as  the  best  computer  for  this  system  because  of 
its  word  3ize  I/O  capability  instruction  set,  memory  capability,  and 
availability  of  software  development  tools  (Ref  7:24). 

Norden's  LSI-11‘1  is  identical,  in  its  functions  to  the  LSI-11 
produced  by  Digital  Electronics  Corporation  (DEC),  but  uses  a  different 
backplane,  thus  making  it  incompatible  with  standard  DEC  systems. 
Therefore,  a  special  backplane  and  card  cage  would  have  to  be 
constructed  in  order  to  use  the  Norden  LSI-11M  with  DEC's  analog 
converter  and  real-time  clock  that  were  selected  for  the  MVDM  (Ref 
7:29-34).  Although  the  Norden  LSI-  11  is  a  more  durable  nodule  (Ref 
7:25),  DEC's  modules  are  designed  to  operate  in  "Minimally  Controlled 
Environmental  Conditions"  (Ref  6  20)  such  as  the  environment  in  which 
the  MFDM  will  be  operating.  Therefore,  the  LSI-11  was  acquired  and 
used  for  the  development  of  this  system.  The  LSI-1 1  Microcomputer  is 
capable  of  handling  the  additional  requirements  by  using  different 
options  than  initially  determined  (Ref  7* 25-35). 


KD l 1 — J  Processor  Module.  The  KD11-F  processor  module  was 
originally  selected  for  the  system,  but  is  not  adequate  for  the  real¬ 
time  processing  of  data  at  the  higher  frequency  range  of  interest  (6 
K  Hz).  This  is  because  the  KD11-F  processor  module  has  a  A  K  by  16 
bit  semiconductor  RAM  within  the  module.  This  memory  requires  the  on¬ 
board  memory  refresh  to  be  enabled  which  results  in  an  execution 
stoppage  for  up  to  130  micro-seconds  during  a  memory  refresh  cycle  (Ref 
8*76).  This  delay  tine,  when  combined  with  a  AO  micro-second  data 
acquisition  service  routine,  results  in  an  execution  delay  of  up  to  170 
micro-seconds.  This  is  then  the  minimum  interval  between  real-time 
events  for  which  it  can  be  insured  that  no  real-time  event  will  be 
lost.  Thus,  the  maximum  sampling  frequency  would  be  5*8  KHz  (1/170 
micro-seconds) . 

The  KD11-J  processor  module  does  not  have  the  on-board  memory,  so 
the  130  micro-second  delay  does  not  exist  because  the  on-board  memory 
refresh  capability  is  disabled.  This  reduces  the  e  miniraun  interval 
for  real-time  events  to  the  time  required  by  the  data  acquisition 
service  routine  or  the  time  required  for  the  analog-to-digital 
converter  to  digitize  a  data  point,  whichever  is  greater. 

The  KD11  series  processors  can  all  satisfy  the  power-up  auto-3tart 
and  the  manual  start  as  specified  by  the  system  requirements  (Ref  Chap 
II).  These  processors  have  three  power-up  modes  (Ref  8*237).  Mode  0, 
which  will  be  used  for  this  systen,  allows  program  execution  to  start 
at  the  address  contained  in  location  2A  if  the  RUN/HALT  switch  is  in 
the  RUN  position  when  power  is  applied.  Otherwise,  when  power  is 
applied  the  system  is  in  the  halt  node  and  the  user  can  manually 
enter  the  required  starting  address  to  start  program  execution. 
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Thus,  the  KD11-J  processor  module  meets  the  requirements  defined 
in  Chapter  II,  although  the  system  must  include  a  separate  memory  for 
the  program  and  data. 

MRV11  UV-PROM.  The  requirement  that  the  system  be  able  to  start 
processing  data  upon  application  of  system  power  means  that  the  program 
must  be  stored  in  a  non-volatile  memory  such  as  core  memory  or 
programmable  read-only  memory  (PROM).  The  contents  of  core  memory  can 
accidentally  be  changed,  thus  the  use  of  core  memory  to  store  the 
program  would  not  be  the  best  solution.  A  UV-PROM  memory  is  a  read¬ 
only  memory  which  makes  it  difficult  to  change  the  program  accidentally. 
The  program  can  be  executed  directly  from  a  PROM  module  even  though 
data  cannot  be  stored  in  the  PROM  memory.  Therefore,  the  best  memory 
configuration  to  use  is  PROM  to  hold  the  program,  and  core  memory  to 
hold  the  data  arrays. 

A  DEC  MRV-11  module,  which  contains  up  to  4  K  words  of  UV-PROM  and 
256  words  of  RAM  is  available  for  the  LSI-11  (Ref  8:254).  Three  of 
these  nodules  (10  K  words)  will  be  required  to  store  the  conplete 
program.  An  advantage  of  UV-PROM  is  that  if  modifications  to  the 
system  are  required,  the  PROM  can  be  re-programmed  with  the  modifica¬ 
tions-  With  the  program  in  PROM,  the  data  arrays  and  system  variables 
have  to  be  contained  in  a  non-volatile  read/write  memory  such  as  the  MMV- 
11  Core  Memory. 

MMV 1 1 -A  Core  Memory .  A  non-volatile  memory  is  required  to  store 
the  data  array  so  that  it  is  available  to  the  user  after  the  test  run 
has  been  performed.  The  MMV-ll  4  K  by  16  bit  core  memory  meets  the 
requirements  of  the  system  (Ref  7' 29).  The  number  of  core  memory 
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nodules  required  Is  determined  by  the  size  of  the  CDF  array*  The 
maximum  size  of  the  CDF  array  is  1440  elements  ((60000  Hz/50Hz)xl2 
temperature  Intervals),  and  the  requirement  that  any  two  elements  In 
the  array  can  differ  by  as  much  as  10  (Ref  Chap  II)  requires  that  two 
words  (floating  point  numbers)  be  used  to  hold  each  element  in  the 
array  (Ref  8:560).  This  requires  2880  words  of  core  memoy  to  hold  the 
CDF  array.  Thus,  one  4  K  core  memory  module  is  adequate. 

ADV11-A  Analog-to-Dlgital  Converter.  The  A-to-D  converter  must 
be  capable  of  processing  analog  data  at  a  rate  equal  to  or  greater  than 
the  Nyquist  sampling  frequency.  For  a  frequency  range  of  0-6  K  Hz  the 
analog  signal  must  be  sampled  at  a  rate  of  at  least  12  K  Hz,  which  is 
the  minimum  sampling  rate  for  the  analog-to-digital  converter. 

In  order  to  obtain  data  at  a  specific  frequency  interval,  the 
sampling  rate  mu3t  be  controllable  so  as  to  satisfy  equation  5.  The 
ADVll-A  (Analog-t o-Digital  Converter)  and  KWV11-A  (Programmable  Real- 
Time  Clock)  are  the  only  modules  made  for  the  LSI-11  computer  which 
allow  the  sampling  rate  to  be  controlled  so  as  to  obtain  the  required 
frequency  intervals. 

The  ADVll-A  is  capable  of  a  21.5  K  Hz  sampling  rate  (Ref  2:23) 
which  exceeds  the  12  K  Hz  sampling  rate  required  for  this  system.  The 
ADVll-A  converter  can  be  controlled  by  the  program  or  by  an  external 
event.  In  order  to  control  the  real-time  data  acquisition,  the 
conversions  can  be  controlled  by  the  KWV11-A  programmable  real-time 
clock. 
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KWV11-A  Programmable  Real-Time  Clock.  This  option  is  available 
for  the  LSI-11  and  can  be  used  to  control  the  analog-to-digltal 
conversions.  It  can  be  operated  in  any  one  of  four  programmable  modes: 
single  interval,  repeated  interval,  external  event  timing,  and 
external  event  timing  from  zero  base  (Ref  2:3-1). 

This  application  requires  that  it  be  used  in  mode  1  -  repeated 
Interval  with  the  interval  being  set  by  the  program.  The  analog-to- 
digltal  converter  with  the  programmable  real-time  clock  will  handle  the 
data  acquisition  under  program  control.  The  remaining  system  hardware 
that  is  required  is  the  hardware  to  interface  the  system  to  the  user. 

DLV1 1-J  Serial  Line  Interface.  The  DLV11-J  is  the  interface  for 
serial  line  devices  such  as  the  terminal  and  printer,  and  the  LSI-11 
bus.  There  are  four  independent  serial  interfaces  on  the  module,  and 
each  channel  can  independently  be  configured  for  the  following:  baud 
rate  (150-38400),  number  of  data  bits  -  7  or  8,  number  of  stop  bits  - 
1  or  2,  and  parity  -  even  or  odd.  Therefore,  this  module  is  capable  of 
interfacing  with  the  user's  terminal  and  printer. 

PD VI 1-B  Ba ckplane.  The  DDV11-B  is  used  to  hold  and  electrically 
connect  the  different  modules.  It  is  used  with  the  H0341  card  cage 
which  supports  the  nodules.  The  backplane  will  accept  the  KD11-J 
processor,  core  memory,  and  up  to  11  option  modules  (Ref  8:50).  This 
backplane/card  cage  combination  will  hold  all  of  the  required  modules 
for  the  system  in  one  cage,  thus  making  it  a  portable,  self-contained 
system.  The  size  of  the  backplane  with  card  cage  is  17.1  inches  long, 
4.8  Inches  high  and  8.7  inches  deep  (Ref  8:51). 

This  completes  the  description  of  the  required  LSI-11  hardware. 
Digital  Equipment  Corporation  recently  Introduced  two  new  versions  of 
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the  LSI-11  -  the  LSI-11/2,  and  LSI-11/23  (Ref  5:55)  -  cither  of  which 


could  be  used  for  the  MFDM  system.  These  two  processors  are  half  the 
width  of  the  LSI-11  but  are  double  the  thickness  (Ref  5:55-56).  As  of 

December  1979,  DEC  Is  not  yet  making  a  correspondingly  smaller  analog- 

< 

to-digital  converter,  real  time  clock,  or  core  memory  so  there  is?  no 
size  sain  by  using  these  newer  processors.  However,  the  LSI-11/23  is 
four  times  faster  (Ref  5:57)  than  the  LSI-11  which  would  increase  the 
cumulative  damage  matrix  update  rate  if  future  requirements  need  a 
faster  rate. 

In  addition  to  the  microcomputer  hardware  required,  some  hardware 
la  needed  to  acquire  the  analog  signals.  This  hardware  consists  of  two 
sensors  and  associated  linear  amplifiers  which  are  discussed  next. 

IBV1 1-A  LSI-1 1 /Instrument  Bus  Interface.  This  module  is  used  to 
allow  the  computer  to  control  the  auto-ranging  gain  of  the  vibration 
signal  amplifier.  It  provides  the  means  of  inhibiting  the  auto-ranging 
gain,  and  reading  the  gain  setting  of  the  amplifier  before  the 
frequency  signal  is  sampled.  Due  to  a  lack  of  time,  this  nodule  was 
never  Integrated  into  the  bench  model  system. 

KPV1 1-B  Power  Fail/Line  Time  Clock/Termlnator.  This  module 
controls  the  power-up  and  power-down  sequences  and  has  provisions  for 
an  external  clock  and  a  remote  console  capability  (Ref  8:339).  With  a 
DC  power  supply,  no  line  clock  will  be  available.  Thus,  the  external 
clock  input  can  be  used  as  an  input  for  the  system  clock.  The  remote 
console  capability  provides  the  necessary  switching  controls  to  allow 
the  system  to  start  up  in  either  the  automatic  mode  or  manual  mode. 
Although  the  module  is  designed  for  an  AC  input,  the  rectifier  circuit 


can  be  bypassed  to  allow  operation  with  DC  power. 

Analog  Sensors 

In  order  to  gather  the  required  data,  two  analog  sensors  are 
required  -  one  thermocouple  and  one  accelerometer.  The  AFFDL 
laboratory  supplied  the  sensors  and  associated  amplifiers  for  this 
project.  The  thermocouple  was  made  from  type  K  thermocouple  wire  and 
the  accelerometer  was  a  Vibra  Metrics  VM-112  integrated  accelerometer 
(Ref  13).  The  accelerometer  had  a  "built  In"  electronic  pre-amplifier 
allowing  its  signals  to  be  transmitted  over  long  cables  of  conventional 
2-wire,  or  coaxial  construction,  with  no  loss  in  sensitivity  or 
spurious  cable  noise  (Ref  13).  The  output  of  each  sensor  goes  into  an 
amplifier  which  amplifies  the  signal  to  a  level  which  can  be  digitized 
by  the  analog-to-digital  converter. 

The  amplifiers  supplied  with  the  sensors  were  Intech  Model  A2583 
wide  band  data  amplifiers  (Ref  1:1).  The  gain  of  the  amplifiers  can  be 
set  either  automatically  or  manually  to  any  one  of  eight  internally 
programmed  gain  ranges.  The  gain  setting  of  the  amplifier  can  be 
determined  by  either  an  analog  voltage  output  or  a  BCD  output,  and  the 
auto  ranging  can  be  inhibited  under  program  control  (Ref  1:1).  The 
thermocouple  was  connected  to  one  amplifier  with  the  gain  set 
to  +60  db.  The  accelerometer  was  connected  to  another  amplifier  which 
used  a  10  position  thumb  wheel  switch  to  determine  the  gain  setting. 
Figure  9  shows  the  output  of  the  thermocouple  amplifier  for  the 
temperature  range  of  interest,  and  Table  V  shows  the  gain  of  the 
accelerometer  amplifier  for  the  10  possible  thumb  wheel  settings. 

Table  VI  shows  the  specifications  for  the  Intech  amplifiers  (Ref  1:6). 
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Lg  9.  Temperature  Signal  Versus  Temperature 


TABLE  V 


Gain  Settings  For  Frequency  Amplifier 


Switch  Setting 

Gain  (db) 

Digital  Gain  Code 

Analog  Cain  Code 

0 

-10 

111 

-2.40 

1 

0 

110 

-1.70 

2 

10 

101 

-1.00 

3 

20 

100 

-0.35 

4 

30 

oil 

0.35 

5 

40 

010 

1.00 

6 

50 

001 

1.70 

7 

60 

000 

2.40 

8 

Auto- 

9 

Ranging 

- 

TABLE  VI 
SPECIFICATIONS 


A2583  Automatic  Gain  Ranging  Amplifier 


Functions 

Range 

Units 

Min 

Max 

Gain  Range  (lOdb  steps) 

-10.0 

+60.0 

db 

Vernier  Gain/step 

-6  5 

+3.5 

db/step 

Gain  Accuracy 

0.1 

db 

Gain  Stability 

0.1 

db 

Cain  change  threshold 

Down  (adjustable) 

0.5 

2.0 

V  R11S 

UP  (adjustable) 

10  0 

20.0 

db  below 

down 

Temperature  Range 

-25  0 

+85.0 

C 

Power  Supply  Voltage 

14.0 

16.0 

VDC 

Current 

60.0 

ma 
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The  design  elso  requires  two  hardware  modules  which  are  not 
commercially  produced  for  the  LSI-11  system.  The  modules  are:  a  DC  power 
supply  and  a  system  clock.  The  DC  power  supply  will  be  designed  and  built 
by  the  AFFDL  at  Wright-Patterson  AFB  (Ref  11),  and  is  not  covered  in  this 
report . 

System  Clock 

The  system  clock  must  produce  a  TTL  level  compatible  square  wave 
with  a  one-second  period.  This  signal  will  be  the  external  clock  input 
to  the  KPV11-B  module.  This  module  will  produce  an  external  interrupt 
which  will  drive  a  software  system  clock.  A  clock  chip  such  as  the 
National  Semiconductor  MM5316  Digital  Alarm  Clock  can  be  used  to 
generate  the  1  second  pulse  and  will  also  provide  a  digital  readout  of 
the  system  time  if  it  is  desired. 

Summary 

This  section  has  described  the  hardware  available  for  the  MFDM. 

The  hardware  previously  selected  (Ref  7:23-35)  was  re-considered  in 
light  of  the  new  system  requirements.  The  KD11-F  and  DLV11  are  the 
only  previously  selected  modules  which  do  not  meet  the  new  require¬ 
ments.  The  KD11-F  processor  module  has  an  on-board  memory  which 
requires  the  use  of  on-board  memory  refresh  and,  therefore,  could  not 
acquire  the  data  at  12  K  Hz  as  required  to  obtain  the  desired  frequency 
range.  Thus,  the  KD11-J  processor  module  was  selected.  It  is  the 
same  as  the  KD11-F  processor  except  there  is  no  on-board  memory,  and 
therefore  the  on-board  memory  refresh  is  not  required.  The  DLV11 
nodule  is  a  single  port  serial  interface  module.  In  order  to  use  a 
terminal  and  printer,  two  modules  would  be  required;  thus,  a  DLV11-J 
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four-port  serial  interface  nodule  was  selected .  This  module  can  be 
used  for  both  the  printer  and  the  terminal. 

Two  additional  modules,  not  previously  considered,  are  also 
required.  The  first  one  is  an  1IRV11  UV-PROM  (Ref  8-137)  which  is 
used  to  store  the  program,  and  the  second  one  is  an  IBV11  LSI- 
11/Instrunent  Bus  Interface  (Ref  8:180)  which  is  used  to  control  the 
auto-ranging  gain  of  the  vibration  signal  amplifier.  These 
additional  modules,  along  with  the  modules  previously  selected, 
establish  the  hardware  requirements  for  the  MFDM.  The  next  chapter 
uses  this  information  in  designing  the  software  to  implement  the 
system  algorithm. 
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V.  Software  Design 


Introduction 

The  next  step  in  the  development  of  the  MFDM  is  to  design  an 
algorithm  for  each  of  the  software  modules.  In  this  chapter,  the 
possible  program  languages  available  to  implement  the  algorithm  are 
discussed,  and  a  decision  is  made  as  to  which  programming  language  is 
the  best  to  use.  Each  of  the  required  software  modules  is  then 
developed,  insuring  that  each  nodule  will  interface  properly  with  the 
other  modules. 

Programming  Languages 

The  language  in  which  a  module  is  coded  affects  the  complexity  of 
the  program,  execution  speed,  and  the  ease  with  which  the  module  can  be 
maintained.  The  programming  languages  which  were  available  on  the  LSI- 
11  development  system  were  RT-11  assembly  language,  RT-11  FORTRAN,  and 
UCSD  Pascal.  The  Pascal  operating  system  required  24  K  words  of  memory 
to  operate  (Ref  12:iii),  and  the  development  system  only  had  20  K  words 
of  memory  which  could  be  used  without  memory  refresh  when  executing  the 
program  (Ref  Chap  IV).  Therefore,  the  Pascal  operating  system  could 
not  be  used  when  executing  the  program.  In  addition,  no  debugging 
routine  was  included  in  the  Pascal  operating  system.  Without  such  a 
routine,  program  debugging  would  be  very  difficult.  Because  of  these 
two  restrictions,  it  was  decided  not  to  use  Pascal  to  develop  the 
software. 

Generally,  a  structured  FORTRAN  program  is  easier  to  understand 
than  a  structured  assembly  language  program  which  performs  the  same 
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task-  A  FORTRAN  program,  being  easier  to  understand,  will  generally  be 
easier  to  maintain  during  its  life;  and  therefore,  it  will  have  a  lower 
life  cycle  cost.  The  disadvantage  of  a  FORTRAN  program  as  compared  to 
an  assembly  language  program  is  that  the  execution  time  is  generally 
longer.  Therefore  it  was  decided  to  code  the  modules  which  require 
the  most  computational  time  or  which  require  real-time  data  processing 
in  assembly  language,  while  the  other  modules  were  coded  in  FORTRAN. 

It  was  decided  to  code  the  EXECUTIVE  and  FOURIER  modules  in  assembly 
language  because  EXECUTIVE  controls  the  real-time  data  aquisition  and 
FOURIER,  which  computes  a  fast  Fourier  transform,  requires  the  greatest 
percentage  of  execution  time  (Ref  7  115).  The  remaining  modules  were 
all  coded  in  FORTRAN  for  ease  in  future  maintenance. 

The  system  algorithm  W3S  constructed  using  top-down  structured 
design,  which  was  considered  by  the  author  to  be  the  best  design 
methodology  available.  The  process  started  with  the  development  of  a 
module  which  controlled  the  entire  process.  This  module  was  then 
expanded  into  subordinate  nodules,  each  of  which  performed  one  major 
function.  These  subordinate  modules  were  then  expanded  into 
subordinate  nodules  if  the  major  function  could  be  reduced  to  a  series 
of  independent  functions.  The  result  of  this  type  of  design  was  a 
structured  system  where  each  module  accomplished  a  specific  function, 
and  the  coupling  between  modules  was  minimized.  A  structure  chart  for 
the  complete  system  is  9hown  in  Figure  10,  and  the  module  interfaces 
are  shown  in  Table  VII.  The  following  paragraphs  describe  each  of  the 
modules  which  make  up  the  system  in  detail,  and  appendices  A  and  D 
contain  flow  charts  and  the  source  code  for  each  of  the  modules. 


Arguments  of  the  MFDM  Modules 
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Configuration  Modules 

The  configuration  modules  are  the  modules  which  initialize  the 
parameters  used  by  the  data  processing  modules.  These  modules  are 
Interactive  modules  which  allow  the  user  to  configure  the  system  for  a 
particular  test.  In  order  to  make  the  configuration  process  easy  for 
the  user  to  accomplish,  the  modules  were  designed  to  display  the 
system's  current  configuration  and  then  prompt  the  user  to  input  the 
first  letter  of  the  function  to  be  executed.  Thus,  by  using  prompting 
messages,  the  system  instructs  the  user  what  to  do  to  configure  the 
system.  The  following  sections  discuss  the  design  of  the  modules  which 
are  used  to  configure  the  system  to  process  data  according  to  the 
requirements  of  the  user. 

MFDM.  Program  MFDM  is  the  controlling  nodule  for  the  entire 
system.  It  is  an  interactive  module  which  is  initally  executed  when 
the  system  is  started  in  the  manual  mode.  It  prints  the  program's 
name  at  the  user's  terminal  and  then  requests  that  the  user  enter  the 
first  letter  of  the  function  to  be  performed.  The  format  of  these  two 
lines  is  as  follows: 

"MICROPROCESSOR  FLIGHT  DAMAGE  MONITOR" 

"ENTER  S(ETUP,  E(XECUTE,  0(UTPUT,  Q(UIT  -  <CR>" 

where  <CR>  means  the  user's  response  followed  by  a  carriage  return.  If 
a  "Q"  is  entered  by  the  user,  program  execution  halts.  If  an  "S",  "E", 
or  "0"  is  entered,  the  appropriate  subordinate  module  is  executed.  If 
an  illegal  entry  is  made,  the  system's  ERROR  module  is  called  with  an 
error  message  (Table  VII,  line  4)  being  passed  to  it.  After  the 
subordinate  module  completes  execution,  control  is  passed  back  to  MFDM. 


i 

\  . 


The  modules  which  are  subordinate  to  MFDM  are  discussed  next 


SETUP.  Subroutine  SETUP  Is  called  when  the  user  wants  to 
configure  the  system  before  a  test.  The  module  first  displays  the  last 
used  system  configuration  and  asks  what  Is  to  be  changed.  If  the  last 
configuration  was  0-3  K  Hz  and  25  Hr  the  following  lines  would  be 
printed  at  the  user's  terminal: 

0-3000  Hz 
25  Hz 

(Last  Time)  Hours 
(Last  Comments) 

(Last  Date) 

CHANCE  F(REQ,  R(UN  DUR,  COMMENTS,  D(ATE,  Q(UIT,  or  V(ERIFY  -  ?<CR> 

The  user's  input  is  then  read  and,  based  upon  validity,  the 
corresponding  module  is  called.  An  "F"  input  calls  the  GETFREQUENCY 
module,  an  "R"  calls  GETTIME,  a  "D"  calls  GETDATE,  and  a  "C"  calls 
GETCOMMENT.  The  arguments  passed  to  the  subordinate  modules  are  shown 
in  Table  VII  lines  8-18.  If  a  "Q"  is  entered,  program  control  returns 
to  MFDM;  if  a  "V"  is  entered,  the  current  system  configuration  is 
printed  at  the  terminal.  If  an  invalid  input  is  made,  the  error 
module  is  called  -  passing  the  error  message  shown  in  Table  VII  line 
il.  After  the  error  message  is  printed,  the  current  system  config¬ 
uration  is  printed  at  the  terminal,  followed  by  the  prompt  line 
requesting  another  input. 

GETFREQUENCY.  Subroutine  GETFREQUENCY  is  called  by  SETUP  when  the 
frequency  range  and/or  interval  is  to  be  changed.  Tills  subroutine 
prints  the  possible  configuration  at  the  terminal  and  requests  which 
configuration  is  desired.  This  is  accomplished  by  printing  the 


FREQUENCY  RANGE: 
FREQUENCY  INTERVAL: 
RUN  DURATION: 
COMMENTS: 

DATE: 
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following  lines: 


FREQ.  RANGE  <K  Hz)  DELTAF  (Hz)  ENTER 


0  - 

6.0 

100 

1 

0  - 

6.0 

50 

2 

0  - 

3.0 

100 

3 

0  - 

3.0 

50 

4 

0  - 

3.0 

25 

5 

0  - 

1.5 

100 

6 

0  - 

1.5 

50 

7 

0  - 

1.5 

25 

8 

ENTER  A  NUMBER  BETWEEN  1  &  8  CORRESPONDING  TO  THE  DESIRED  SETUP  ?<CR> 


An  Integer  Is  then  read  from  the  terminal  and  used  In  an 
arithmetic  go-to  statement  which  passes  control  to  an  appropriate 
initialization  section  corresponding  to  the  desired  configuration.  If 
the  first  character  entered  is  not  an  integer  between  I  and  8,  program 
control  falls  through  the  go-to  statement,  and  the  error  routine  is 
called  with  the  message  shown  in  Table  VII  line  39.  After  this,  the 
configuration  table  is  re-printed  at  the  terminal  and  the  user  is 
prompted  to  enter  an  Integer.  When  a  valid  entry  is  read,  control 
passes  to  the  corresponding  label  where  the  following  variables  are 
initialized: 


COMMON  LOCAL 

B  EXP 

INDX  INDEX 


DEFINITION 

The  power  of  2  which  gives  the  number 
of  data  (N)  points  to  be  acquired 
(N-2**EXP)  -  used  in  EXECUTIVE. 

An  integer  index  used  with  the  sine- 
look-up  table  in  the  subroutine  FFT. 
It  is  initialized  to  the  offset 
required  for  the  look-up  table. 


INTVL 


INTERVAL 


DELTAP 


MAXFREQ 


RTCMODE  CLKHODE 


RTCRATE  CLKRATE 


The  number  of  frequency  Intervals 
In  the  selected  frequency  range 
(INTERVAL-MAXFREQ/DELTAF  +1).  It 
Is  used  In  the  DAMAGEFACTOR  and 
CDFOUTPUT  subroutine. 

The  frequency  interval  for  which  the 
CDF  is  to  be  computed.  It  is  used 
in  the  CDFOUTPUT  subroutine. 

The  maximum  frequency  for  which  the 
CDF  is  to  be  computed.  It  is  used 
in  the  CDFOUTPUT  subroutine. 

CLKMODE  is  a  local  variable 
initialized  with  a  data  statement 
to  11  (base  10).  This  value  sets 
the  real-time  clock  to  1  MHz 
(Ref  x7:3-12).  It  is  used  in  EXECUTIVE. 

A  variable  which  controls  the 
frequency  with  which  the  A/D 
converter  samples  the  vibration 
signal.  The  CLKMODE  (initialized  to 
1  MHz  with  a  data  statement  in  this 
routine)  divided  by  the  CLKRATE 
determines  the  sampling  rate. It  is 
used  in  EXECUTIVE. 


After  initializing  these  variables,  thus  re-configuring  the  system, 
program  control  returns  to  MFDM. 

GETTIME.  Subroutine  GETTIME  is  called  when  the  user  wants  to 
change  the  length  of  time  for  which  the  system  is  to  process  data.  The 
arguments  passed  to  and  from  this  routine  are  shown  in  Table  VII  line 
9.  The  argument  "MINUTES"  is  a  local  variable  corresponding  to  the 
"PTIME"  (Planned  Time)  common  block.  "PTIME"  is  the  global  assembly 
language  variable  used  in  the  assembly  language  routine  EXEC.  Thus, 
when  the  value  of  "MINUTES"  is  changed,  the  global  variable  is  also 
changed. 

This  routine  first  prompts  the  user  to  enter  the  desired  run 
duration.  In  hours,  by  printing  the  following  line  at  the  terminal: 


“ENTER  DURATION  OF  TEST  (HOURS)  EXAMPLE  15.25  <CR> 

The  user  then  enters  the  desired  run  duration  and  it  is  stored  in  the 
variable  "DURATION".  This  value  is  tested,  and  if  it  is  zero  or  less, 
the  error  module  is  called,  passing  to  it  the  error  message  shown  in 
Table  VII  line  32,  and  a  new  time  is  requested.  When  a  positive  time  is 
entered  it  is  multiplied  by  60  to  form  the  number  of  minutes  the  test 
will  run,  which  is  stored  in  the  global  variable  PTIME.  Program 
control  then  goes  back  to  the  calling  module,  returning  the  new 
DURATION. 

GETDATE.  Subroutine  GETDATE  is  called  when  the  user  desires  to 
change  the  test  date.  The  purpose  of  the  test  date  is  to  allow  the 
user  to  have  the  date  of  the  test  printed  in  the  heading  of  the  CDF 
table. 

This  subroutine  prompts  the  user  to  enter  the  desired  date  with 
the  following  message: 

ENTER  THE  NEU  DATE  Ex.  14-DEC-79  <CR> 

The  user  then  responds  by  entering  the  date  -  which  can  be  up  to  13 
characters  long.  The  date  is  stored  in  the  array  DATE  (13)  of  type 
BITE.  No  error  checking  is  performed,  and  if  the  user  enters  more  than 
13  characters,  the  input  is  truncated  after  13  characters. 

GETCOMMENT .  Subroutine  GETCOMMENT  is  the  same  as  subroutine 
GETDATE  except  that  an  array  COMMENTS  (40)  is  first  cleared  before  it 
Is  loaded  with  the  user's  supplied  comments.  As  with  GETDATE,  the  only 
purpose  of  the  COMMENTS  is  to  allow  the  user  to  have  Identifying 
comments  printed  in  the  heading  of  the  CDF  table.  The  user  is  prompted 
with  the  following  line: 
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ENTER  COMMENTS  -  40  CHARACTERS  MAX  <CR> 

and  the  user  then  enters  up  to  40  characters .  If  more  than  40 
characters  are  entered,  the  input  string  is  truncated  to  the  first  40 
characters. 

ERROR.  Subroutine  ERROR  is  called  by  MFDM,  SETUP,  GETFREQUENCY, 
GETTIME,  and  XECUTE  when  an  illegal  input  is  entered.  The  argument 
passed  to  this  routine  consists  of  an  error  message  up  to  40  characters 
long  and  in  single  quotation  marks,  printed  at  the  user's  terminal. 
ERROR  concatenates  "ERROR"  to  the  beginning  of  the  message,  prints  it 
at  the  terminal,  and  returns  program  control  to  the  calling  module.  If 
a  message  sent  to  this  routine  contains  more  than  40  characters,  the 
message  is  truncated  to  the  first  40  characters. 

CDFOUTPUT.  This  subroutine  is  called  when  the  cumulative  damage 
matrix  is  to  be  printed.  It  is  designed  to  operate  with  a  Heathkit 
H14  printer,  and  uses  its  software  selectable  font  size  feature  to 
control  the  size  of  the  print.  The  table  heading  is  printed  at  10 
characters  per  inch,  while  the  table  entries  are  printed  at  12  char¬ 
acters  per  inch.  The  table  is  printed  as  two  sections  ♦  the  first  for 
a  temperature  range  from  -50  to  100  degrees  Fahrenheit,  and  the  s 
second  for  a  75  to  225  degree  Fahrenbbtt  range.  The  75  and  100  degree 
temperature  intervals  were  repeated  on  both  sections  to  make  interpo¬ 
lation  in  this  temperature  range  easier. 

In  addition,  each  section  consists  of  up  to  six  pages,  depending 
on  the  frequency  range  and  Interval  selected,  with  41  lines  per  page. 
This  allows  one  page  of  the  table  to  be  printed  on  an  8  1/2  by  11 
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inch  piece  of  paper.  The  table  heading  is  printed  on  the  top  of  every 
page  and  consists  of : 


CUMULATIVE  DAMAGE  MATRIX  PAGE  (page) 

DATE  -  test  date 

DURATION  -  test  duration 

COMMENTS  -  user's  comments 

Following  the  table  heading,  one  page  of  the  table  is  printed.  An 
example  of  the  output  is  shown  in  Table  IX,  and  is  explained  in 
Chapter  VII. 

The  variables  used  in  this  routine  are  defined  as  follows: 


CDF  (120,12) 
CONTENTS  (40) 

DATE  (13) 

DELTAF 

FREQPAGE 

I,J,K,L 

LARGE 

LOW,  HIGH 


RICHFREQ 


The  cumulative  damage  factor 
array 

A  string  for  any  comments 
the  user  wants  on  the  output 

A  string  to  hold  the  date  which 
will  be  printed  on  the  output 

The  frequency  between  elements 
of  the  CDF  array 

The  frequency  range  per  page  of 
output 

General  purpose  indexes 

A  three  character  alpha  constant 
which,  when  sent  to  the  printer, 
causes  SO  characters  per  line  to 
be  printed 

Two  indexing  limits  used  to  con¬ 
trol  the  printing  of  the  temper¬ 
ature  scale.  Their  values  are: 

LOW-1-50  deg  F;  1IICH-7-100  deg  F 
LOW-6-75  deg  F;  HIGH-12-225  deg  F 

Index  limits  to  control  the  print¬ 
ing  of  the  frequency  in  the  table 
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H1CHTEMP 


LOWTEMP, 

LTSYM 

MAXFREQ 

NFREQ 

NPAGE 

RMINUTES 


Index  llraitfi  to  control  the  print¬ 
ing  of  the  temperature  in  the  table 

The  symbol  to  the  left  of  the  first 
temperature  in  the  table  heading. 

It  is  either  a  space  or  a  <  (less  than) 

The  maximum  frequency  in  the  6DF  array 

The  number  of  frequency  components 
in  the  CDF  array 

The  number  of  pages  required  to 
print  the  table  for  one  temperature 
range 

The  actual  duration  of  the  test  run 


RTSYM  The  symbol  printed  to  the  left  of 

the  last  temperature  on  the  table 
heading.  It  is  either  a  space  or  a 
>  (greater  than) 


After  printing  the  complete  table  program  control  returns  to 
the  MFDM  module,  where,  if  the  user  desires  another  copy  of  the 
output  table  can  be  obtained  by  entering  another  "0". 

This  completes  the  discussion  of  the  configuration  modules.  The 
next  section  discusses  the  execution  modules. 

Execution  Modules 

The  requirements  defined  in  Chapter  II  state  that  there  be  two 
methods  to  initiate  data  processing.  One  method  required  that  the  user 
be  capable  of  starting  data  processing  from  a  terminal,  and  the  second 
method  required  that  data  processing  commence  with  the  application  of 
power  to  the  system.  The  XECUTE  module  is  used  to  initiate  data 
processing  from  a  terninal.  The  PURUPDN  module  is  used  to  Initiate 
data  processing  when  power  is  applied  if  the  processor's  RUN/HALT 
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switch  is  in  the  run  position  (Ref  Chap  IV).  An  explanation  of  the  two 
modules  follows: 

XECUTE .  Subroutine  XECUTE  is  an  interactive  module  for  user- 
controlled  data  processing.  The  variables  used  in  this  subroutine  are 
defined  below: 


Global  variable: 


PTIME 


Local  variables: 


ANSWER 


MINUTES 


DURATION 


A  label  common  block  with  length  of  one 
word  used  to  store  the  planned  time 
(in  minutes)  for  which  data  processing 
is  to  take  place. 


A  variable  which  contains  the  user's 
response.  It  should  contain  be  either 
a  "Y"  or  an  "N".  Anything  else 
results  in  an  error. 

A  local  variable  equal  to  PTIME 

The  planned  processing  tine  in  hours. 


When  called  by  MFDM,  this  procedure  first  checks  that  PTIME  is 
greater  than  zero.  If  it  isn't,  GETTIHE  is  called  to  get  a  new  run 
time.  The  new  run  time  Is  then  printed  at  the  terminal  followed  by  a 
prompt  asking  if  it  needs  to  be  changed.  The  user's  input  is  then 
read,  and  if  it  is  an  "N"  (NO),  then  the  following  message  is  printed 
at  the  terminal: 


"DATA  BEING  PROCESSED" 


and  program  control  passes  to  the  EXECUTIVE  module  where  data 
processing  is  initiated.  If  a  "Y"  (YES)  is  entered,  then  GETTIME  is 
called  so  a  new  time  can  be  entered.  Afer  the  new  time  is  entered. 
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program  control  returns  to  this  module  and  the  new  time  Is  printed  nt 
the  terminal*  The  user  is  then  asked  if  it  needs  to  be  changed.  This 
second  printing  of  the  time  Is  to  provide  a  "Inst  chance"  recovery 
capability  in  case  a  wrong  time  is  entered.  When  the  user  responds 
with  a  "Y"  or  an  "N’’,  the  process  is  repeated.  If  anything  other  than 
a  "Y"  or  an  "N"  is  entered,  the  error  module  is  called  with  the  error 
message  given  In  Table  VII  line  27.  After  data  processing  Is  complete, 
control  returnn  to  this  module  and  a  "DATA  PROCESS  INC  COMPLETE"  message 
is  printed  at  the  terminal.  Program  control  then  returns  t:o  M!'I)M. 

PWRUPDN.  Subroutine  PUR UPON  controls  the  start  of  data  processing 
when  the  system  starts  in  the  "automatic  node".  When  the  processor's 
RUN/HALT  switch  Is  in  the  "RUN"  position  and  power  in  applied  to  the 
system,  this  routine  is  entered  via  the  interrupt  vector  at  location  24. 

The  global  varlnbles/lahela  which  are  used  in  this  subroutine  are 
defined  below: 

GLOBAL  VARIABLES /LAB ELS  DEFINITION 

BEGIN  The  label  of  the  second  instruction 

in  subroutine  EXEC.  If  power  fails 
and  comes  up  in  the  middle  or  data 
processing,  control  returns  to  this 
point  * 

EXEC  The  label  of  the  first  Instruction 

In  subroutine  EXEC.  Initial  data 
processing  starts  at  this  location. 

PTIME  A  global  variable  containing  the 

planned  duration  of  the  test. 

PWRFLG  A  global  variable  used  to  indicate 

if  it  is  the  first  power-up  or  not. 

PTIME  A  global  variable  used  to  store 

the  total  time  the  system  has  been 
processing  data. 
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This  routine  first  sets  location  24  to  the  address  of  the  power¬ 
down  sequence  (PWRDN)  in  case  of  a  power  failure,  and  then  checks  to 
see  if  it  is  the  first  power-up  sequence*  If  it  is,  PWRFLG  is  set 
equal  to  one,  and  program  control  passes  to  subroutine  EXECUTIVE.  If 
it  is  not  (PWRFLG  -  1),  then  PTIME  and  RTIME  are  compared.  If  RTIME  is 
greater  than  or  equal  to  PTIME,  program  execution  halts.  Otherwise, 
program  control  jumps  to  label  "BEGIN"  in  EXECUTIVE. 

After  the  first  power-up  sequence,  location  24  contains  the 
address  of  the  power-down  (PWRDN)  sequence.  When  power  fails,  the 
power-down  sequence  loads  location  35  with  the  address  of  the  power-up 
sequence  (PWRUP)  and  executes  a  halt  instruction.  This  allows  for 
multiple  power  up/down  sequences  during  the  data  processing  period. 

Data  Processing  Modules 

The  modules  discussed  in  this  section  -  EXECUTIVE,  FFT , 
POWERSPECRUM,  and  DAMACEFACTOR  -  are  the  data  acquisition  and 
processing  modules.  The  functions  which  these  modules  must  perform 
are:  acquiring  and  digitizing  the  analog  vibration  and  temperature 
signals,  computing  the  fast  Fourier  transform  of  the  vibration  signal 
producing  the  frequency  spectrum  which  the  test  item  is  "seeing", 
computing  the  power  in  each  of  the  frequency  intervals  which  make  up 
the  frequency  spectrum,  computing  the  damage  factor  for  each  frequecy 
Interval,  and  summing  the  damage  factors  into  a  cumulative  damage 
matrix.  The  FFT  module  was  taken  from  the  previous  work  on  this  system 
(Ref  7:75-85)  and  modified  to  meet  the  additional  requirements  of  this 
system  (Ref  Chap  II).  Two  modules,  POWERSPECTRUM  and  DAMACEFACTOR,  were 
vrltten  in  FORTRAN  so  as  to  make  future  modifications  easier  -  should 
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they  be  required  -  while  EXECUTIVE  was  written  in  assembly  language  for 
speed  of  execution  for  the  data  acquisition.  The  design  of  the  four 
modules  will  now  be  explained. 

EXECUTIVE.  The  subroutine  EXECUTIVE  is  the  controlling  module  for 
the  data  processing.  When  called  from  XECUTE,  or  the  first  time  called 
from  PWRUPDN,  a  macro-setup  routine  (INITAL)  is  called  to  clear  the 
power  spectrum  array  (PWRSPEC)  used  in  POWERSPRC,  clear  the  CDF  array 
(DATA)  used  in  DAMACEFACTOR,  and  to  initialize  the  global  variables 
used  in  the  data  processing  modules.  INITAL  uses  the  global  variable  B 
to  initialize  the  following  macro-global  variables: 


BITCOUNT  A  counter  used  in  the  FFT  module's  "bit  inversion 
routine". It  is  initialized  to  half  the  number  of 
binary  bits  required  to  write  the  value  of  N. 

COUNT  A  counter  used  in  FFT  and  initialized  to  half  the 

number  of  "complex"  data  points. 

INTCELNUM  An  initial  value  for  the  local  variable  "CELNUM" 

which  is  used  in  FFT's  "2nd  PASS  AND  ONWARD"  routine. 


N  The  number  of  real  data  points  to  be  processed.  It 

is  used  in  EXECUTIVE  and  FFT  and  is  equal  to  the 
value  2N  discussed  in  Chapter  III.  (A  macro  variable 
cannot  begin  with  a  number). 

N2  The  number  of  "complex"  data  points.  It  is  equal 

to  N/2  and  is  used  in  FFT,  POWERS PECTRUM, 
DAMACEFACTOR. 

NM  An  index  used  in  FFT's  "POST  PROCESSING"  routine. 

It  is  Initialized  to  the  relative  location  of  the 
N*lst  "complex"  data  point. 

PWRFLG  A  flag  set  by  PWRUPDN  upon  the  first  system  power 

up  in  the  auto-start  mode.  It  is  initialized  to 

sero. 

RTIME  A  variable  used  in  EXECUTIVE  to  hold  the  actual 

time  data  has  been  processed.  It  is  initialized  to 
ierq. 


TIME  A  variable  used  in  the  EXECUTIVE'S  time  Interrupt 

<  routine*  It  is  used  as  a  "60's"  counter  to  count  60 

seconds.  It  is  initialized  to  zero. 

TSUM  A  variable  used  in  EXECUTIVE  to  hold  the  sum  of  ten 

temperature  values,  which  is  used  to  set  the  global 
variable  TFLG.  It  is  initialized  to  zero. 

In  addition  to  initializing  these  variables,  INITAL  also  loads  the 
time  and  power  interrupt  vectors  with  their  proper  values.  The  next 
function  of  EXECUTIVE  is  to  control  the  acquisition  of  the  temperature 
data. 

The  temperature  data  is  acquired  using  channel  05  of  the  analog-to 
digital  converter.  This  is  accomplished  by  writing  2401  (octal)  into 
the  analog-to-digital  (A/D)  command-status  register  (ADCSR).  This 
selects  channel  05  as  the  temperature  channel  and  starts  the  A/D 
conversion.  After  the  conversion  is  completed,  the  4000  (octal)  bias 
Is  subtracted,  and  the  converted  data  is  summed  into  the  variable  TSUM, 
which  is  used  to  determine  the  temperature  interval  in  which  to  store 
the  damage  factor.  After  acquiring  the  temperature  data,  EXECUTIVE 
must  control  the  acquisition  of  the  vibration  data.  The  amplifier 
used  to  amplify  the  vibration  signal  has  an  auto-ranging  feature  which 
adjusts  the  gain  of  the  amplifier  depending  on  the  magnitude  of  the 
vibration  signal.  Before  any  vibration  data  can  be  digitized,  the  auto 
ranging  feature  must  be  disabled  -  freezing  it  at  one  gain  for  the 
sampling  period  -  and  then  enabling  it  after  the  sampling  period. 

This  ability  to  Inhibit  the  auto-ranging  gain  has  not  yet  been 
Implemented,  due  to  a  lack  of  time.  It  will  require  the  use  of  the 
IBV11-A  LSI-11  Instrument  Bus  Interface  module.  First,  a  logical  ”0" 
will  have  to  be  written  to  a  bit  in  the  instrument  bus  data  register, 
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which  is  connected  to  the  inhibit  line  of  the  vibration  signal 
amplifier.  This  will  inhibit  the  auto-ranging  feature  of  the 
amplifier,  allowing  one  set  of  data  points  to  be  acquired  at  the  same 
gain  setting. 

After  inhibiting  the  auto-ranging  gain,  the  BCD  gain  setting  can 
be  read  from  the  instrument  bus  data  register  and  converted  into  an 
integer  gain  multiplication  factor.  The  multiplication  factor  is 
stored  as  the  variable  "ADGAIN"  which  is  used  in  the  POWERSPECTRUM 
module.  After  acquiring  one  set  of  data  po^fs,  the  auto-ranging 
feature  must  be  released  by  writing  a  logical  "1"  to  the  inhibit  line. 
This  will  allow  the  gain  of  the  amplifier  to  vary  with  the  intensity  of 
the  input  signal  until  the  next  data  set  is  required. 

Since  this  function  has  not  yet  been  implemented  in  EXECUTIVE,  the 
ADGAIN  has  been  set  to  1  In  EXECUTIVE.  After  obtaining  the 
multiplication  factor,  EXECUTIVE  loads  the  ADCSR  with  7040  (octal), 
which  selects  channel  16  as  the  channel  to  be  used  for  the  vibration 
signal,  and  disables  A/D  interrupts.  The  real-time  clock  (RTC)  i3  then 
initialized  aa  required  by  the  GETFREQUENCY  module  and  "N"  vibration 
data  samples  are  acquired  and  stored  in  the  array  (ADDATA)  under 
program  wait.  The  reason  for  using  program  wait  rather  than  interrupt- 
driven  data  acquisition  is  that  there  is  a  44  micro-second  delay  after 
an  interrupt  is  requested  before  the  first  interrupt  service  routine 
instruction  is  executed  (Ref  8:B-5),  during  which  time  nothing  is 
done.  This,  combined  with  a  40  micro-second  service  routine  (the 
length  of  the  routine  used  in  EXECUTIVE),  would  mean  the  fastest 
sampling  rate  would  be  11.9  K  Hz  which  is  not  fast  enough  (Ref  Chap 
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XII)*  With  program  wait,  the  maximum  sampling  rate  (22.2  K  Hz)  is 
determined  by  the  execution  time  of  the  service  routine  (approximately 
45  micro-seconds),  or  the  analog-to-digital  conversion  time  (37  micro¬ 
seconds),  whichever  is  longer.  After  the  sampling  of  "N"  data  points, 
the  real-time  clock  is  stopped  and  the  damage  factor  is  computed. 

EXECUTIVE  then  calls  FFT  followed  by  POWERSPECTRUM.  The  sequence  - 
sample  data,  call  FFT,  call  POWERSPECTRUM  -  is  repeated  ten  times  in 
order  to  compute  an  average  power  spectrum  (Ref  Chap  III).  After  this, 
the  value  in  TSUM  is  compared  to  a  series  of  pre-set  values  which 
correspond  to  the  temperature  intervals.  These  values,  located  at 
label  TINTVL,  are  obtained  by  taking  ten  times  the  theoretical 
digitized  temperature  signal  for  the  temperature  intervals  of 
Interest.  After  each  comparison  where  TSUM  was  greater  than  the  value 
in  TINTVL,  TFLG  is  incremented  by  one.  Thus,  TFLG  will  take  on  a  value 
between  1  and  12  corresponding  to  the  average  temperature  of  the  test 
item*  The  DAMAGE  module  is  then  called  to  compute  the  cumulative 
damage  factor  and  update  the  CDF  array.  After  this,  the  variables 
TFLG,  TSUM,  and  the  power  spectrum  array  PWRSPC  are  cleared,  the  run 
time  is  checked  and  if  RTIME  (actual  run-time)  does  not  equal  PTIME 
(planned  run-time)  the  procedure  is  repeated.  When  RTIME  equals  PTIME, 
program  control  is  passed  back  to  the  module  which  called  EXECUTIVE  - 
either  XECUTE  or  PWRUPDN. 

FFT.  This  module  was  taken  from  the  previous  work  performed  on  a 
cumulative  damage  acquisition  system  (Ref  7:75-85)  and  modified  to 
allow  for  different  values  of  "N"  as  defined  in  the  EXECUTIVE  section 
just  discussed*  Reference  7  should  be  consulted  to  obtain  the  theory 
and  the  actual  equations  implemented  in  this  routine*  This  section 
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will  cover  the  modifications  that  were  made  to  the  module  to  meet  the 
requirements  defined  in  Chapter  II* 

Originally,  the  module  was  designed  to  compute  a  fast  Fourier 
transform  (FFT)  of  32  data  points.  In  the  present  system,  an  FFT  of 
32,64,128,  or  256  data  points  may  be  required.  To  accomplish  this, 
constants  which  were  originally  initialization  constants  were  changed 
to  variables  initialized  by  GETFREQUENCY  and  by  INITAL  in  EXECUTIVE. 
These  global  variables  are:  1NDX  -  initialized  by  GETFREQUENCY;  and 
BITCOUNT,  COUNT,  INTCELHUM,  N,  N2,  and  NM  -  initialized  by  INITAL  in 
EXECUTIVE.  Each  of  the  variables  is  defined  in  the  EXECUTIVE  or 
GETFREQUENCY  section  of  this  chapter,  and  Table  VIII  shows  the  values 
to  which  these  variables  are  initialized. 

The  largest  modification  necessary  was  the  modification 
of  the  sine-cosine  look-up  table  which  is  used  in  the  ’’2nd  pass  and 
onward"  routine  and  in  the  "post  processing"  routine.  The  "2nd  pass 
and  onward"  routine  requires  the  sine  and  cosine  values  of 
(2  *  PHI  *  K)/N2  for  K  -  0,1,2,...,(N  2/2)  -  1,  and  the  "post 
processing"  routine  requires  the  sine  and  cosine  values  of 
(2  *  PHI  *  K)/N  for  K  -  0,1,2,...,  (N/2)  -  1.  One  look-up  table  was 
constructed  containing  the  sine  and  cosine  values  of  (PHI  *  K)/12S  for 
K  ■  0, 1, 2, . . . , 127  -  which  contains  all  of  the  required  values  for  N  - 
32,64,128,256.  The  appropriate  value  can  then  be  extracted  from  the 
table  by  adjusting  the  arguments  passed  to  the  look-up  routine  according 
to  N.  This  adjustment  is  accomplished  with  the  global  constant  INDX. 

The  exact  details  of  how  this  is  done  can  be  found  in  the  source  code 
end  in  Reference  7. 
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TABLE  VIII 


Initialized  Constants  Used  in  FFT 


N 

N2 

INTCELNUM 

NM 

INDX 

BITC0UNT 

256 

128 

64 

32 

508 

2 

7 

128 

64 

32 

16 

236 

4 

6 

64 

32 

16 

8 

124 

8 

5 

32 

16 

8 

4 

60 

16 

4 

This  modification,  along  with  the  variables  initialized  by  the 
other  modules,  allows  the  FFT  nodule  to  function  with  values  of  N  = 
32,64,128,  or  256.  Thus,  the  FFT  module  meets  the  requirements  defined 
in  Chapter  II  in  addition  to  allowing  frequency  ranges /intervals  of  0- 
10  K  Hz  at  100  Hz,  0-6  K  Hz  at  50  Hz,  and  0-3  K  Hz  at  25  Hz  (N*=256) 
which,  if  desired,  can  be  implemented  in  the  future.  The  output  of 
FFT,  the  frequency  spectrum  of  the  vibrational  signal,  is  stored  as 
complex  fixed-point  values  in  the  input  array  ADDATA.  The  binary  point 
is  located  between  bits  14  and  15,  and  the  scale  factor  for  the  array  is 
stored  in  "SFACTR".  After  completing  the  FFT,  program  control  is 
passed  back  to  EXECUTIVE  which  then  calls  POWERSPECTRUM  to  compute  the 
power  spectrum. 

POWERSPECTRUM .  After  the  FFT  of  the  input  data  is  calculated,  the 
instantaneous  power  spectrum  for  each  frequency  interval  in  the  FFT 
spectrum  must  be  computed,  except  for  the  zero  frequency  interval.  The 
power  spectrum  for  the  zero  frequency  interval  is  not  computed  because 
it  is  undefined  (computation  would  involve  division  by  zero).  The 
computed  Instantaneous  power  spectrum  is  then  summed  into  the  array 
PWRSPEC.  The  FFT  and  instantaneous  power  spectrum  calculations  are 


repeated  for  ten  vibration  data  sets,  thus  forming  ten  flues  the 
average  power  spectrum  in  the  global  array  PWRSPEC  (Ref  7:73.85) 


This  module  was  coded  in  FORTRAN  because  of  the  case  of  coding  and 
future  modification.  A  FORTRAN  macro  was  used  to  convert  the  fixed 
point  frequency  spectrum  values  to  floating  point  values.  The  macro 
assumes  the  fixed  point  value  is  an  integer,  converts  it  to  a  floating 
point  number,  and  multiplies  it  by  the  scale  factor  and  the  amplifier's 
gain.  The  macro  is  encoded  as  follows: 

FPDATA(l)  -  DATA(I)  x  2(scale)x  GAIN  ,  (11) 

where  DATA(I)  Is  the  fixed  point  value  (treated  as  an  integer),  SCALE 
is  the  scale  factor  from  the  FFT  module,  and  GAIN  is  the  gain  of  the 
amplifier  ADGAIN  which  is  determined  in  EXECUTIVE.  The  square  of  the 
instantaneous  power  spectrum  is  computed  by  squaring  the  real  and 
complex  floating  point  values  of  the  frequency  spectrum  values  and 
adding  them  together.  This  value  is  then  summed  into  the  local  power 
spectrum  array  -  PTOSPECTRUM  (global  array  PWRSPC).  This  is  represented 
mathematically  by  Eq  12: 


POWERS PECTRUM ( I )  -  POWERS PECTRUM ( I )  +  (FPDATA(K))2  +  (FPDATA(K+1))2  , 

1  -  1  2,3,..-NFREQ,  K-  3,5,7,...  (12) 

where  FPOATA(K)  and  FPDATA(K+1)  are  the  real  and  iraaginery  values  of 
the  complex  frequency  spectrum  data,  and  Imax  is  the  maximum  number  of 
frequency  Intervals  in  the  desired  frequency  range. 


The  output  of  this  module  Is  the  sura  of  the  squares  of  the 
Instantaneous  power  spectrums.  After  ten  summations  the  array 
PWRSPECTRUM  contains  ten  times  the  square  of  the  power  spectrum  - 
which  is  the  data  passed  to  the  damage  module* 

DAMAGEFACTOR .  This  module  conputes  the  damage  factor  for  each 
frequency  interval,  of  interest,  in  the  FFT  spectrum  and  suns  the  values 
into  a  cumulative  danage  matrix.  The  input  arguments  to  this  nodule  are 
the  power  spectrum  array  from  POUERSPECTRUM  and  the  temperature  flag 
TFLG  from  EXECUTIVE.  Since  the  power  spectrum  array  actually  contains 
ten  times  the  square  of  the  average  of  the  power  spectrum,  Eq  1  is  re¬ 
written  as: 


DAMAGEFACTOR (n)  -  <K(A' (n))3' 3)/f (n)5*/'  n  -  1,2,3,.*  . NFREQ  (13) 

where  A' (n)  -  (PONRESPECTRUM(n) )°* 5/10  ,  (14) 

n  +  0  because  DAMAGEFACTOR (0)  is  undefined,  and  NFREQ  is  the 
maximum  number  of  frequency  intervals  in  the  frequency  range  of 
Interest.  Since  K  is  a  constant  for  all  n,  and  the  danage  factor 
is  a  relative  value  it  can  be  left  out.  Re-writing  Eq  14  without  the 
"K"  tern  and  substituting  for  A'(n)  yields: 

DF(n)  -  (PURSPEC(n)1,65)/f(n)5*4,  n  -  1,2,3,..- , NFREQ  (15) 

3.3 

where  the  10  *  term  was  left  out  because  it  is  also  a  constant  for  all 
«•  Eq  15  was  coded  as  a  FORTRAN  macro  DAMAGE  (n,f(n))  and  was  used 
to  sum  the  damage  factor  into  the  CDF  matrix.  Using  n  local  variable 
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TEMP  as  the  global  variable  TFLG  the  equation  for  the  CDF  matrix  Is: 

CDF(n.TEMP)  -  CDF (n, TEMP)  +  DF(n)  ,  n-  1,2,3,** . ,NFREQ  (16) 

Eqs  15  and  16  were  implemented,  and  the  resulting  cumulative  damage 
matrix  is  an  nmax  by  12  array.  The  F^T  routine  is  coded  for  a  maximum 
nmax  of  120  (6000  Hz  range  with  a  frequency  interval  of  50  Hz),  which 
means  the  largest  CDF  array  is  120  by  12,  or  1440  elements,  which 
requires  2880  16  bit  words.  In  order  to  operate  with  values  of  N  between 
32  and  256,  the  CDF  array  was  defined  as  CDF  (120,12).  For  values  of  N 
less  than  256,  only  part  of  the  array  is  used. 

Summary 

A  top-down  structured  design  was  used  to  design  each  of  the 
software  nodules  needed  to  implement  this  system.  The  structure  chart, 
shown  in  Figure  10,  and  the  module  interface  table  shovm  in  Table  VII, 
describe  the  system  as  a  whole  while  appendices  A  and  8  contain  the 
flow  charts  and  source  code  for  each  of  the  modules  The  next  step  in 
implementing  the  system  is  the  individual  testing  of  the  software 
modules,  followed  by  the  assembly  and  testing  of  the  complete  system. 
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VI.  Testing  and  Experimental  Results 


Introduction 

This  chapter  discusses  both  the  individual  testing  of  the  software 
modules  and  the  construction  and  testing  of  a  bench  model  system.  The 
experimental  results  obtained  from  the  bench  model  testing  are  then 
discussed,  and  an  explanation  of  the  CDF  output  table  is  given. 

Testing  of  Software  Modules 

Each  of  the  14  software  modules  was  individually  tested  for  proper 
operation  by  writing  a  controlling  program  which  executed  the  nodule 
being  tested,  and  simulated  the  other  modules.  After  testing  each 
module,  the  "configuration"  modules  were  tested  as  a  unit  by  using  a 
dummy  execution  routine  which  simulated  the  processing  of  data.  This 
"configuration"  system  was  tested  using  the  RT-11  operating  system's  On- 
Line-Debugging-Technique  (ODT)  (Ref  8:645)  to  control  program  execution. 
Using  ODT,  the  program  execution  was  halted  before  and  after  each 
subroutine  call  in  order  to  verify  the  proper  operation  of  the 
subroutines.  After  verifying  the  proper  operation  of  each  nodule  in  the 
system,  the  interactive  functions  were  checked  for  proper  operation  and 
variable  initialization.  When  an  error  was  discovered,  the  module  where 
the  error  occurred  was  checked  and  corrected. 

After  completing  the  tests  on  the  "configuration"  system,  the 
"data  processing"  modules  were  tested  as  a  unit.  First,  proper  array 
indexing  was  verified  and  then  the  calculation  of  the  damage  factor  was 
checked  by  a  manual  calculation.  This  completed  the  software  testing 
of  the  MFDM.  The  next  step  was  the  Integration  of  the  hardware  and  the 
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software  into  a  complete  system. 

Integration  of  the  Hardware  and  the  Software 

After  completing  the  software  testing,  the  hardware  and  the 
software  was  configured  Into  a  complete  system  which  was  then  tested 
for  functional  accuracy. 

Bench  Model  Configuration.  In  order  to  assemble  a  bench  model 
system,  the  RD11-J  processor  board  had  to  have  the  4  K  of  on-board 
memory  disabled.  This  was  accomplished  by  installing  wire  wrap  posts 
for  jumper  W9  which  must  be  installed  to  disable  a  reply  from  the 
resident  memory  (Ref  8:234).  The  processor  jumpers  were  then 
configured  as  a  KD11-J  processor  (Ref  8:234).  The  4  K  core  memory  was 
set  to  be  the  first  bank  of  memory  and  16  K  of  semiconductor  memory  was 
set  to  memory  banks  one  through  four.  A  DDVll-B  backplane  was  not 
available,  so  two  H9270  backplanes  were  used  -  connecting  them  together 
with  a  BCV1B  option  nodule  (Ref  8:380).  Figure  11  shows  the  system 
configuration  used  for  the  bench  model. 

Bench  Model  Testing.  To  test  the  computational  accuracy  of  the 
bench  model,  it  was  taken  to  an  AFFDL  laboratory  where  both  a  shaker 
table  and  spectrum  analyzing  equipment  were  available.  The  lab's 
analog  data  processing  system,  which  computes  the  cumulative  damage 
factor,  was  not  working,  so  a  B  and  K  Instrument's  graphic  spectrum 
equalizer  was  used  to  simulate  the  vibrational  input  signal.  A 
Spectral  Dynamics  Corporation's  one-third  octave  analyzer  was  connected 
to  the  input  jignal  in  order  to  determine  the  input  spectrum.  There 

t 

was  no  equipment  available  to  test  the  temperature  signal  acquisition 
function,  so  the  temperature  input  leads  were  shorted  together 
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simulating  a  0  voltage  Input.  After  the  MFDM  processed  the  simulated 
data.  It  was  determined  that  the  resolution  of  the  spectrum  analyzer 
was  so  much  lesj  than  that  of  the  MFDM  that  no  estimations  of  the 
accuracy  of  the  MFDM  were  possible. 

The  testing  of  the  system  was  then  continued  at  the  Air  Force 
Institute  of  Technology  where  the  bench  model  was  designed  and 
constructed.  The  input  signals  were  simulated  by  using  four  sine-wave 
generators  to  simulate  the  vibration  signal,  and  a  DC  power  supply  was 
used  to  simulate  the  temperature  signal.  The  accuracy  of  the 
temperature  scale  was  determined  by  using  a  test  routine  which  printed 
the  computed  temperature  interval  corresponding  to  the  DC  voltage  being 
used  as  the  simulated  temperature  input  signal.  By  varying  the  DC 
voltage,  the  voltage  at  which  the  computed  temperature  interval  changed 
was  determined.  From  this  voltage,  the  corresponding  temperature  at 
which  the  computed  temperature  interval  changed  was  computed.  This 
test  determined  that  the  temperature  scale  for  the  CDF  output  table  is 
accurate  within  one  degree. 

The  vibration  data  was  checked  by  usinj,  four  sine-wave  generators 
as  simulated  input  signals.  Sine  wave  generators  were  used  in  order  to 
have  a  known  frequency  spectrum  for  an  input  to  the  system.  This 
allowed  the  resolution  of  the  system  to  be  determined  by  varying  the 
frequencies  of  the  input  signals  until  the  frequencies  were  so  close  as 
to  make  them  indistinguishable.  The  smallest  difference  in  frequencies 
fot'  which  the  two  signals  were  distinguishable  was  the  resolution. 

Using  this  method,  it  was  determined  that  the  resolution  of  the  CDF 
output  was  twice  the  frequency  interval  of  the  output  table.  The 
roasoa  for  this  large  resolution  is  the  following:  if  two  input 


frequencies  are  within  two  frequency  Intervals  of  each  other,  the 
output  table  will  show  a  peak  at  both  frequencies.  Since  the  two 
frequency  peaks  are  next  to  each  other,  it  is  impossible  to  determine 
if  there  was  one  frequency  input  between  the  two  intervals,  or  one 
frequency  input  in  each  interval.  By  using  only  one  frequency  as  an 
input  signal,  it  was  deteimined  that  the  system's  accuracy  is  within  15 
Hz  when  using  the  25  Hz  frequency  interval,  within  25  Hz  when  using  the 
50  Hz  frequency  interval,  and  within  50  Hz  when  using  the  100  Hz 
frequency  interval. 

Experimental  Results 

Without  an  analog  data  processing  system  to  compute  the  damage 
factor,  no  accurate  method  of  checking  the  computational  accuracy  of 
the  MFDM  was  available.  Also,  without  an  audio  frequency  spectrum 
analyzer,  no  accurate  method  of  checking  the  output  spectrum  was 
available.  This  left  only  the  testing  that  was  done  with  the  simulated 
input  signals.  Table  IX  shows  the  results  of  one  of  the  tests  using  a 
system  configuration  of  0-3  K  Hz  and  an  interval  of  25  Hz.  The  input 
vibration  signal  consisted  of  500  Hz,  835  Hz,  875  Hz,  and  2685  Hz  sine 
waves.  The  amplitudes  were  3,5,4,  and  1  volt  respectively.  By 
varying  the  input  frequency,  it  was  determined  that  the  frequency 
resolution,  for  one  input  signal,  is  +/—  20Hz  up  to  6000  Hz.  A 
temperature  resolution  was  not  obtained  because  the  thermocouple's 
linear  amplifier  did  not  work. 

Explanation  of  CPF  Table.  The  CDF  Table  shown  in  Table  IX  is  an 


example  of  the  output  of  the  MFDM.  The  date  and  comments  in  the  table 


TABLE  IX 


Cumulative  Damage  Table  Output 


CUMULAT I UE  DAMAGE  MATRIX  PAGE  1 

DATE  12-NOU-79 

DURATION  14.72  HOURS 

COMMENTS  5G0  HZ  <*3U,  335  *S5U,  375  «4U,  2635  SlU 
FREQ  <H2>  TEMPERATURE  <DEG  F> 


<-50 

-25 

3 

25 

50 

75 

100 

25. 

0.000E+00 

0. 0005*00 

0.0005+00 

O. 156E+06 

0. 0005+00 

0.6605+00 

0.0005+00 

5(3. 

0. 0O0E+00’ 

0. 000E-00 

O.OOOE-OO 

0. 370E+0'4 

0. OOOE+OO 

0. 0005+06 

6. OOOE+OO 

75. 

6. 600E+06 

0. 000c -00 

0. 0000+00 

O. 4245+03 

01.0005+00 

0.000E+O0 

0. 000E+00 

ieo. 

6.  OijijE+OO 

0. 000c +00 

0. 000E+00 

0. 834E+02 

0. OOOE+OO 

O. 000E+60 

0. 0005+60 

125. 

0. 600E+06 

0. 000c -00 

0.0065+60 

6. 2865+02 

0.0605+00 

0. OOOE+OO 

0. 0OOE+0O 

156. 

0.O60E+CO 

6. 0665-00 

0 . 0005+601 

0.1135+02 

0. OOOE+OO 

0. 0005+00 

0. OOOE+OO 

175. 

0. OOOE+OO 

0.000c -00 

0.0005+00 

0. 5275+01 

0.0005+00 

0 . 000E+00 

0. 000E+00 

206. 

O. 6OOE+60 

0. 000E+00 

0.0005+00 

6.272E+61 

0 . OOOE+OO 

0 . OOOE+OO 

6 . 600E+00 

225. 

6 . 0G6E+00 

0. 000E+00 

0.0005+00 

O. 1565+01 

0.0005+00 

O. 0005+00 

0. 606E+00 

256. 

O. 60OE+60 

0.0005+00 

0.0005+00 

0. 863E+00 

0.0005+00 

0. OOOE+00 

0.000E+O0 

275. 

6. 660E+00 

0. 000c- 00 

0.0005+00 

0. 644E+00 

0. OOOE+OO 

0.6005+00 

0. 600E+00 

360. 

0. 0O0E+00 

0. OOOE-OO 

0 . OOOE+OO 

0.4515+00 

0. 6605+60 

0.0065+00 

0. 000E+00 

325. 

6.60OE+O0 

0. OOOE+OO 

6 . 0005*60 

0. 33-E+06 

0. OOOE+OO 

0.0005+60 

0.00OE+O0 

350. 

6. 660E+00 

0.0005+00 

0.000E+00 

0. 257E+00 

0. OOOE+OO 

6.0005+00 

0. 600E+00 

375. 

0. OOOE+00 

0 . OOOE+OO 

0.6005+66 

0.214E+00 

0 .  OOOE+OO* 

0 . OOOE+OO 

0.006E+O0 

-106. 

0. 00OE+0O 

0. 0005*00 

0. 0065+06 

O. 1 035+00 

0, OOOE+OO 

0. 0005+00 

0. 000E+00 

425. 

0. 600E+00 

0. OOOE+OO 

0.0005+00 

6. 170E+O0 

0 , OOOE+OO 

0 . 0065+00 

0 . 006E+6'0 

450. 

6.  OO0E+O0 

0.0005*00 

0. 0005+00 

0.204E+00 

0.060E+60 

0 . OOOE+OO 

O.  O0OE+00 

475. 

0. O00E+60 

0. 0005+00 

0.0005+00 

0.431E+00 

6.6605+06 

0.6005+00 

O. 600E+O0 

506. 

0.60OE+00 

0.  0005*00 

6.0605+00 

0. 4005+04 

6. 0605+00 

0.0005+60 

0.600E+60 

525. 

0. 006E+G0 

■0.0065+00 

0.0005+010 

0 . 257E+O0 

0. OOOE+OO 

0.6006+00 

0. 060E+O0 

556. 

6 . 660E+O0 

0.0005+00 

0 . 0005+00 

0.  l6'SE+00 

0.0005+0O 

0 . 0005+00 

0. OOOE+OO 

575. 

6. 0005+06 

0. 0005+00 

0. 0065+00 

0.0465-0! 

01.  0005+60 

0 . 0006+00 

0. OOOE+OO 

600. 

0. 0O0E+O0 

0.0005*00 

0.0605+06' 

6.0835-01 

0 . OOOE+OO 

0.0005+00 

O .  OOOE+OO' 

625. 

0. 000E+06 

0.0O05-00 

0. 0005+00 

0. 1 l 15+00 

0.  OOOE+OO 

0. 0005+00 

0. 000E+00 

656. 

6. 066E+00 

0.000E+00 

0. 0005+00 

6. 1355+60 

0. OOOE+OO 

0. 0005+00 

0. 006E+60 

675. 

0. 000c +00 

0. 000c -00 

0 . 000E+00 

0. 177E+00 

0. OOOE+00 

0.0005+00 

0. 000E+00 

700. 

0. 0O0E+00 

0. 0005+00 

0. 6065+00 

0. 253E+00 

0.O0OE+O0 

0.0006+00 

0.000E+00 

725. 

0 . 600’E+OO 

0. 0005*00 

0. 0005+00 

0.412E+0O 

0. OOOE+OO 

0. 0005+00 

0. 060E+00 

756. 

0. 0005*00 

0.  0005+00 

0. 0005+06 

0. 864E+O0 

0 . OOOE+OO 

0.0005+00 

0. 000E+00 

775. 

6.00OE+00 

0.0065*60 

6.0005+00 

6.2155+01 

0.000E+00 

0.0605+00 

0.  OOOE+OO 

800. 

0. 0005+00 

0.  0005*00 

0. 00O5+O0 

0. 108E+02 

0. OOOE+OO 

0. 000E+00 

6. 000E+00 

825. 

0. 000E+00 

0.  0005+00 

•j.  0005+00 

0. 67 iE+03 

0 . OOOE+OO 

0 . OOOE+06 

0 .  OC'OE+60 

850. 

0. OOOE+00 

0. 0605*00 

0.0005+00 

6.  ! 085+03 

0. 000E+OO 

0, OOOE+OO 

O.  0605+00 

875. 

0. 0OOE+0O 

0.0005*00 

0. 0005+00 

6. 50SE+03 

0. OOOE+OO 

0.O6OE+00 

0. 6005+00 

860. 

6. 00OE+00 

0 . 000- *<-00 

0.000E+00 

0.7015+60 

O. OOOE+OO 

0.0005+00 

0.606E+00 

025. 

0.00OE+O0 

0. 0005*00 

0. 000E+00 

6. 268E+00 

0. 000E+0O 

0.  0005+00 

0.  060E+00 

050. 

0 . O00E+O6' 

0.60O5+0O 

6.  C’06’5+00 

6.8HE-01 

6. 6605+00 

0.006E+0O 

O. OOOE+OO 

075. 

0.  000E+00 

0.O005+O0 

0. 0005+00 

0. 37*4£-0l 

O. OOOE+OO 

0 . OOOE+OO 

0. OOOE+OO 

1000. 

0.  OOOE+OO 

0. 000 £-00 

0.0005+00 

6. 1035-01 

0. OOOE+00 

0 . 0005+00 

0. 0'X'E+00 

1025. 

0. 600E+00 

0. 0005*00 

0.0005+00 

O’!.  1065-01 

0.  OOOE+OO 

0.6005+00 

O.O0OE+00 

FREQ 

<-50 

<HZ> 

-25 

C  25 

TEMPERATURE  <DEG 

50 

F> 

75 

100 

73 


TABLE  IX 


Cumulative  Damage  Table  Output  (cont) 


CUHULATIUE  DAMAGE  MATRIX 

PAGE 

2 

DATE 

12--N0U 

-79 

DURATION  M.72 

HOURS 

COMMENTS  500  H2 

«3U,  S35 

■S5U,  STS 

MU,  26S5 

«1U 

FREQ 

•H2> 

TEMPERATURE  <DEG 

F> 

<-50 

-25 

0 

25 

50 

75 

100 

1025. 

0. OOOE+OO 

0. OOOE+OO 

0 .  O'.’OE+OO 

0. lOOE-Ol 

0 . OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO  • 

1050. 

0. OOOE+OO 

0.  OOOE+OO 

0 . OOOE+OO 

0. 020E-02 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  I 

1075. 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. 3SSE-02 

0 . OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  i 

1100. 

0. OOOE+OO 

0. OOOE+OO 

0 . OOOE+OO 

0.24SE-02 

O.OOO'E+OO 

0. OOOE+OO 

0. OOOE+OO  I 

1125. 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0.  lt.GE-02 

ft .  ftftftE +00 

O.GOOE+OO 

0. OOOE+OO  i 

1150. 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0. 1 13E-02 

0. OOOE+OO 

0.  OftOE+ftft 

0. OOOE+OO  ■ 

1175. 

0. OOOE+OO 

0 . 000 E*  W 

0. OOOE+OO 

0 . 7,X'E~03 

0. OOOE+OO 

0 . OOOE+OO 

0. OOOE+OO  ! 

1200. 

0.  OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. 5bME-ft3 

0 . OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  I 

1225. 

0;  OOOE+OO 

0. OOOE+OO 

0 . OOOE+OO 

0.41 5E-03 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  i 

1250. 

0, OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0. 307E-03 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  i 

1275. 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0.231E-03 

0. OOOE+OO 

0.  OOOE+O’O 

O. OOOE+OO  i 

1300. 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. 175E-03 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  i 

1325. 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0. 135E-03 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO  ! 

1350. 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0, 1GSE-03 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  ! 

1375. 

0. OOOE+OO 

o*occs^oo 

0. OOOE+OO 

0.03 lE-04 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  ! 

1400. 

0. OOOE+OO 

l>.  OOOE+OO 

0. OOOE+OO 

0.07 IE-04 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  i 

1425.  • 

0. OOOE+OO 

0. OOOE-OO 

0. OOOE+OO 

0. 535E-04 

0 . OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  ! 

1450. 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0.42SE-04 

0.  OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  ! 

1475. 

0.  OOOE+OO 

0.  WOE  ♦CO 

0. OOOE+OO 

0. “40E-04 

0 . 000 E”*  00 

0.  OOOE+OO 

0. OOOE+OO  i 

1500. 

0. OOOE+OO 

0.  OOOE+OO 

0 . OOOE+OO 

0. 2*7E-04 

0. OOOE+OO 

O. OOOE+OO 

O. OOOE+OO  i 

1525. 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0.251E-02 

0.  OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  i 

1550. 

0. OOOE+OO 

Tj  +  COOE-OO 

0. OOOE+OO 

0. 207E-02 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  I 

1575. 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0. 171E-02 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  i 

1000. 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0. 144E-02 

O.OOOE-OO 

0. OOOE+OO 

0. OOOE+OO  ! 

1025. 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. 122E-02 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  i 

1050. 

0. OOOE+OO 

0. OOOE+OO 

0 . OOOE+OO 

0. l0cE-02 

0. OOOE+OO 

0 . OOOE+OO 

0. OOOE+OO  1 

1075. 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0. lOlE-02 

0. OOOE+OO 

0 . OOOE+OO 

0. OOOE+OO  i 

1700. 

0. OOOE+OO 

0. 000 E* 00 

0. OOOE+OO 

0. 753E-03 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  i 

1725. 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0.047E-G3 

0. OOOE+OO 

0. OOOE+OO 

O. OOOE+OO  i 

1750. 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0.  t>»XE“03 

0. OOOE+OO 

0. OOOE+OO 

O. OOOE+OO  I 

1775. 

0. OOOE+OO 

0. OOCE+OO 

0. OOOE+OO 

0. +02E-03 

0. OOOE+OO 

0 . OOOE +00 

0. OOOE+OO  ! 

1000. 

0. OOOE+OO 

ft, OOOE+OO 

0. OOOE+OO 

0.417E-03 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  ! 

1025. 

0. OOOE+OO 

0. OOOE+OO 

0. 3G0E-03 

O.OOOE-OO 

0.  OOOE+OO 

0. OOOE+OO  i 

1050. 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. 310E-O3 

0.  OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  1 

1075. 

0.  OOOE+OO 

0. OOOE+OO 

0. 204E-03 

0.  OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  ! 

1*00. 

0. OOOE+OO 

0. 0*X*E*00 

O'.  OOOE+OO 

0.240E-03 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  ! 

1*25. 

0.  OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO 

0.2i*E-03 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO  i 

1*50. 

0.  OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0.  P>JE“03 

0. OOOE+OO 

0. OOOE+OO 

O. OOOE+OO  ! 

1*75. 

0.  OOOE+OO 

0. OOOE+OO 

0 . 00»>£  *00 

0. 170E-03 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  ! 

2000. 

0.  OOOE+OO 

0. OOOE+OO 

0.  iX'OE +00 

0. I55E-03 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO  i 

2025. 

0. OOOE+OO 

•  0. OOOE+OO 

0. OOOE+OO 

0. l+OE-03 

0.  OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO  1 

2050. 

0. OOOE+OO 

0. OOOE+OO 

0. OOOE+OO 

0. l2tE-03 

0. OOOE+OO 

0.  OOOE+OO 

0. OOOE+OO  1 

<-50 

-25 

0 

25 

50 

75 

too 

FREQ 

<HZ> 

TEMPERATURE  <DEG 

F> 

TABLE  IX 


Cumulative  Damage  Table  Output  (cont) 


CUMULAT I UE  DAMAGE  MATRIX  PAGE  3 

DATE  12-NOU-79 

DURATION  M.  72  HOURS 


COMMENTS 

500  H2  «3U,  £35  «5U,  875  «4U> 

2685 

eiu 

FREQ  <HZ> 
< 

TEMPERATURE 

-50  -25  G 

<DEG 

25 

F> 

50 

75 

100 

2050. 

0. OOOE+OO. 

O.OOOE+OO 

O . OOOE+OO 

0. 126E-03 

U.OOOE+UO 

0.000E+00 

o.oooe+oo  : 

2075. 

0. 000E+00 

O.OOOE-OO 

O.OOOiE+OO 

6. i 15E-03 

0. 

0.00OE— 00 

O.OOOE+OO  ! 

2100. 

O.OOOE+OO 

0.  OOOE+OO 

O.OOOE+OO 

0. lO+E-03 

O.OOOE+OO 

O.OOOE+OO  ! 

2125. 

0. OOOE+OO 

0.  OOOE-O’O 

0*.  OOOE+OO 

O. 858E-04 

0. OOOE+OO 

O.OOOE+OO 

O.OOOE+OO  i 

2150. 

0. 000E+00 

0.  OOOE-OO 

0.  OOOE-O'O 

0.857E-C’4 

O.OOOE+OO 

0.000E-0O 

O.OOOE+OO  i 

2175. 

O.OOOE+OO 

O.OOOE-OO 

0*.  OOOE+OO* 

0. 785E-04 

O.OOOE+OO 

0. OOOE+OO 

0. OOOE+OO  ! 

2200. 

0.000E+00 

O.  OOOE-OO 

O.O'OOE+OO 

6.737E-04 

O.OOOE+OO 

0.000E+00 

O.OOOE+00  ! 

2225. 

0. 000E+00 

O.  OOOE+OO 

O.OO’OE+OO 

O.607E-04 

0 . OOOE+OO 

O.OOOE+OO 

O.OOOE+OO  ! 

2250. 

0. OOOE+OO 

O.OOOE-OO* 

O.OOOE+OO 

6. 635E-04 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO  1 

2275. 

O.OOOE+OO 

0.  OOOE-OO 

O.OOOE+OO 

0. 602E-04 

0 . OOOE+OO 

0.000E+O0 

0.00OE+00  i 

2300. 

O. OOOE+OO 

O.OOOE-OO 

G.O0OE+O0 

0.564E-04 

0. OOOE+OO 

O.OOOE-OO 

O. OOOE+OO  i 

2325. 

O.OOOE+OO 

0.  OOOE-OO 

0. OOOE+OO 

0.544E-04 

O.OOOE+OO 

0.  OOOE+OO 

O.OOOE+OO  i 

2350. 

O.OOOE+OO 

O.OOOE-OO 

O.OOOE+OO 

6.526E-0’4 

0 . OOOE+OO 

0 . 0O0E- GO 

O.OOOE+OO  ! 

2375. 

0. OOOE+OO 

0.  OO'OE-O0 

O.O’OOE+OO 

0.515E-04 

0 .  OOOE-rOO 

0.0O0E+00 

O.OOOE+OO  1 

2-100. 

O.OOOE+OO 

O.OOOE-OO 

0 .  OOOE+OO* 

0.5l4E-0a 

O.OOOE+OO 

0’ .  OOOE-OO 

O.OOOE+OO  : 

2425. 

0. OOOE+OO 

0.0O0E+O0 

O.OOOE+OO 

0.53 IE-04 

O.OOOE+OO 

O.OOOE-OO 

O.OOOE+OO  i 

2450. 

0. OOOE+OO 

0.  OOOE-OO 
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heading  are  user-supplied  and  are  printed  In  the  same  format  as  Input 
by  the  user*  The  duration  is  the  length  of  time  the  system  actually 
processed  data*  If,  for  some  reason,  data  processing  was  stopped 
before  the  expiration  of  the  planned  duration,  the  duration  printed  on 
the  table  would  be  the  actual  time  -  not  the  planned  time. 

Each  column  of  the  table  corresponds  to  a  25  degree  Fahrenheit 
temperature  range,  the  center  of  the  range  being  the  temperature 
printed  at  the  top  of  the  column.  If  the  temperature  of  the  test  item 
is  below  -37*5  degrees  the  cumulative  damage  factor  will  appear  in  the 
column  labelled  <-50.  Also,  if  the  temperature  is  greater  than  217*5 
degrees  .  the  cumulative  damage  factor  will  appear  in  the  column 
labelled  >225.  The  temperature  range  of  -50  to  +225  Is  split  into  two 
ranges:  -50  to  +100  and  +75  to  +225  in  order  to  fit  the  table  on  an  8 
1/2  x  11”  sheet  of  paper.  The  75  and  100  degree  columns  have  been 
reproduced  in  both  ranges  to  make  analysis  of  the  data  at  the  page 
boundary  easier.  Each  row  of  the  table  contains  the  CDF  at  the 
frequency  printed  in  the  left-hand  column.  After  printing  41  lines, 
the  temperature  scale  is  re-printed  and  a  new  page  is  started.  Once 
the  complete  frequency  spectrum  is  printed  for  the  lower  temperature 
range,  the  spectrum  is  printed  for  the  high  temperature  range,  thus 
completing  the  CDF  Table. 

Summary 

As  stated  earlier,  due  to  a  lack  of  time  and  of  the  proper  teat 
equipment,  it  was  impossible  to  give  the  MFDM  system  a  complete 
testing.  The  testing  with  the  simulated  input  signals  shows  that  the 
system  is  working  properly,  although  it  was  not  possible  to  test  the 
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accuracy  of  Its  output*  The  system  still  needs  to  be  tested  with  an 
audio  spectrum  analyzer,  and  a  test  with  an  analog  data  processing 
system  should  be  performed* 
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VII.  Proposed  Flight -Worthy  System 


Introduction 

The  previous  chapters  have  discussed  the  design  and  development  of 
a  bench  model  MFDM.  This  chapter  presents  a  proposed  flight-worthy 
MFDM  system.  Proposed  hardware  will  be  discussed,  Including  that 
required  to  control  the  auto-ranging  gain  function  of  the  vibration 
signal  amplifier,  which  was  not  implemented  in  the  bench  model.  The 
program  modifications  required  to  execute  the  program  from  PROM  are 
also  discussed.  In  conclusion,  the  power  requirements,  as  well  as  a 
mechanical  description  of  the  proposed  MFDM  system  will  be  presented. 

Hardware 

The  proposed  MFDM  system  configuration,  shown  in  Figure  12,  uses  a 
DEC  DDV11-B  Backplane  and  an  H0341  card  cage  to  hold  the  hardware 
modules.  The  modules  are  located  in  the  card  cage  according  to  their 
priority  -  the  processor  being  the  highest  priority,  and  the  serial 
Interface  being  the  second  highest.  The  position  of  the  remaining 
modules  does  not  matter  since  they  do  not  require  any  specific 
priority  In  addition  to  holding  the  required  hardware  modules,  the 
H0341  has  enough  space  to  hold  the  hardware  required  for  the  signal 
amplifiers  and  for  the  system  clock.  This  allows  all  of  the  required 
hardware  to  be  in  one  chassis,  making  it  easier  to  package  (Ref  11). 

The  hardware  not  implemented  in  the  bench  model  was  the  controller 
for  the  auto-ranging  gain  of  the  vibration  signal  amplifier.  The  gain 
must  be  constant  during  the  data  sampling  period,  and  allowed  to  auto- 
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range  when  data  Is  not  being  sampled  (Ref  11).  This  allows  the 
Intensity  of  the  signals  to  vary  without  a  loss  of  data.  If  the  gain  of 
the  amplifier  is  held  constant  for  the  duration  of  the  test,  then  the 
Intensity  of  some  signals  may  be  too  small  to  be  digitized,  while  the 
intensity  of  some  other  signals  may  be  too  large  for  the  analog-to- 
digital  converter. 

The  auto-ranging  function  adjusts  the  amplifier  gain,  within 
limits,  to  keep  the  input  signal  to  the  A/D  converter  less  than  the  5 
volts  maximum,  thus  keeping  it  within  the  allowable  limits  of  the  A/D 
converter  (Ref  2:2-1).  The  IBVll-A  Instrument  Bus  Interface  Module  is 
designed  to  control  external  devices  using  the  "Instrument  Bus"  (Ref 
8:314). 

In  this  case,  the  module  allows  the  program  to  control  the 
vibration  signal  amplifier.  In  addition.  If  future  uses  of  the  system 
require  the  control  of  multiple  signal  amplifiers,  this  module  is 
capable  of  controlling  them  (Ref  8:315).  In  order  to  control  the 
vibration  signal  amplifier's  auto-ranging  gain  the  Inhibit  line  of  the 
amplifier  can  be  controlled  by  one  bit  of  the  data  register,  and  the  3 
bit  gain  code  can  be  read  from  the  data  register  (Ref  8:318).  The  code 
necessary  to  use  this  module  will  be  discussed  in  the  next  section. 

Software 

In  order  to  control  the  IBVll-A  Bus  Interface  Module,  some 
assembly  language  code  must  be  added  to  the  EXECUTIVE  module.  Before 
the  "N"  data  points  are  sampled,  the  auto-ranging  feature  must  be 
inhibited,  and  the  gain  must  be  stored  as  the  multiplication  factor 
"ADGAIN".  This  value,  which  has  been  set  to  one  in  EXECUTIVE  until 


80 


the  actual  code  la  Implemented,  is  the  value  by  which  the  data  is 
multiplied  in  the  POWERSPECTRUM  module  to  scale  the  data* 


The  valid  range  for  "ADGAIN"  was  determined  by  computing  the 

damage  factor  for  the  largest  and  smallest  possible  values  of  the 

power  spectrum  at  the  highest  and  lowest  frequency,  respectively,  and 

then  determining  what  multiplication  factor  would  cause  data  overflow 

or  underflow.  The  largest  value  for  any  one  element  of  the  power 

spectrum  array  is  107  (Ref  Chap  V),  which  when  used  in  Eq  15  with  a 

frequency  of  25  Hz,  gives  a  damage  factor  on  the  order  of  107.  Since 

28 

the  largest  floating  point  number  the  computer  can  store  is  10  (Ref 
8:560),  the  largest  multiplication  factor  ("ADGAIN")  is  102* 

(1028/107). 

The  smallest  "ADGAIN"  was  determined  by  using  the  smallest 

possible  value  (1)  for  an  element  of  the  power  spectrum  array  to 

compute  the  damage  factor  at  a  frequency  of  6000  Hz.  The  damage  factor 

-21 

obtained,  using  Eq  15,  was  10  ,  and  the  smallest  number  which  can  be 

—28 

repesented  by  a  floating  point  number  is  on  the  order  of  10  (Ref 
8:560).  Therefore,  the  smallest  "ADGAIN"  is  I0~7  ( 10~28/ 1 0~2 1 ) .  If  a 
smaller  value  is  used,  the  resulting  damage  factor  will  be  zero  due  to 
the  16  bit  word  size. 

After  the  "N"  data  points  are  sampled,  the  gain  must  be  released 
to  auto-range  until  the  next  data  sample  is  required.  A  flow  chart 
shoving  the  control  of  the  auto-ranging  gain  of  the  signal  amplifier  is 
shown  in  figure  13. 

In  addition  to  the  code  necessary  to  control  the  amplifier,  the 
program  must  be  stored  in  MRVll-BA  UV-PROM  modules.  Because  a  PROM  is 
s  read-only  memory,  the  program  must  be  linked  so  that  all  variables 
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Flow  Chart  for  Auto-Ranging  Gain 


have  addresses  in  the  core  memory*  Also,  there  are  five  instructions 
which  cannot  be  executed  with  operands  located  in  PROM.  These 
Instructions  are  MTPS,  MUL,  DIV,  ASH,  and  ASHC  (Ref  8:403).  Only  two 
modules  use  these  instructions.  EXECUTIVE  uses  MTPS;  and  FFT  uses  MUL, 
ASH,  and  ASHC.  EXECUTIVE  is  coded  so  that  it  will  be  executable  from 
PROM,  but  the  FFT  module  must  be  modified  to  store  the  values  for  "SINE" 
and  "COS"  in  the  core  memory.  This  can  be  accomplished  by  defining  the 
"SINE"  and  "COS"  variables  in  the  PSECT  section  of  the  EXECUTIVE  module 
where  the  CDF  nrray  is  defined  (Ref  Appendix  B).  With  this  change,  all 
of  the  nodules  included  in  the  MFDM  system  will  be  executable  from  the 
PROM  memory. 

Figure  12  shows  3  MRV11-  BA  (4  K  by  16  bit  PROMS),  although  the 
program  can  be  contained  in  2  modules  •  8  K  words.  The  third  PROM  is 
Included  for  future  program  expansion.  Although  not  a  requirement,  it 
is  desired  to  have  the  CDF  matrix  plotted  in  "3  dimensions"  (Ref  10). 

A  program  to  accomplish  this  can  be  stored  in  the  third  PROM. 

Power  Requirements 

The  power  requirements  for  each  hardware  nodule  are  given  in  Table 
X  (Ref  8:43-46).  The  total  power  required  is  the  sum  of  the  values 
given  in  Table  X  plus  the  power  required  by  the  amplifiers  and  clock. 

The  total  power  for  the  modules  is  15.4  amps  at  5  volts  DC  and  2.7  amps 
at  12  volts  DC,  while  the  amplifier  requires  .06  amps  at  15  volts  DC 
and  the  clock  requires  .003  amps  between  8  and  29  volts  DC.  This  power 
is  to  be  supplied  by  rechargeable  batteries,  and  they  should  be  capable 
of  operating  the  system  for  up  to  20  hours  (Ref  Chap  2). 
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TABLE  X 


Proposed  Module  Specifications  (Ref  5:43-46) 


Module 

Module  No. 

Description 

Power  Requirements 
+5V  +12V 

KT-ll-J 

M7264 

LSI-11  Processor 

2.0  A 

0.9  A 

MMV11-A 

H223 

4  K  x  16  Core  Memory 

7.0  A 

0.6  A 

DLV11-J 

M8043 

4-Channel  Asynchronous 

1.3  A 

0.2  A 

Serial  Line  Unit 

Interface 

MRV11-BA 

M8021 

UV-PROM-RAM 

0.7  A 

0.5  A 

KWV11-A 

M7952 

Programmable  Real-Time 

Clock 

1.8  A 

0.1  A 

ADV11-A 

AO  12 

16  Channel  12-bit  A/D 

Converter 

2.0  A 

0.5  A 

IBV11-A 

M7954 

Instrument  Bus  Interface 

0.1  A 

KPV11-B 

M8016YB 

Power  Fail/Line  Time 

Clock/Terminator 

0.6  A 

Dnvil-B 

Backplane 

H0341 

Card  Cage 

RrAD-AO0O  417 
UNCLASSIFIED 


AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OH  SCHOO— ETC  F/G  1/3 
DESIGN  ANO  IMPLEMENTATION  OF  A  MICROPROCESSOR  FLIGHT  DAMAGE  MON— -ETC (U) 
DEC  79  B  C  JOHNSON 
AFIT/GCS/EE/79-6 
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Packaging 

Figure  14  shows  the  author's  concept  of  the  MFDM's  final  shape, 
and  the  layout  of  the  subunits  which  make  up  the  system  is  shown  in 
figure  15.  Some  things  which  must  be  considered  when  designing  the 
package  are  the  cooling  and  damping  of  the  card  cage,  the  strength  of 
the  case  and  the  method  by  which  it  will  be  fastened  to  the  airplane. 
The  two  switches  on  the  front  panel  are  keylock  switches  The  reason 
for  this  type  of  switch  is  that  the  "RUN/HALT"  switch  can  be  locked  in 
the  run  position  when  automatic  start-up  is  required.  Then,  when  the 
DC  power  switch  is  turned  on,  data  processing  will  start  automatically. 

The  use  of  keylock  switches,  although  not  required,  will  prevent 
the  inadvertent  application  or  removal  of  the  system's  power,  or  turn¬ 
ing  the  system  on  with  the  "RUII/HALT"  switch  in  the  "HALT"  position 
when  it  should  be  in  the  "RUN"  position.  Manual  start  up  can  be 
accomplished  by  turning  the  DC  power  switch  on  with  the  "RUII/HALT" 
switch  in  the  "HALT"  position  and  then  switching  it  to  the  "RUN" 
position,  or  by  entering  the  starting  address  of  the  MFDM  from  a 
terminal.  The  package  shown  in  Figures  13  and  14  does  not  include  the 
power  pack  which  may  be  a  separate  box.  A  seperate  power  supply  will 
reduce  the  size  and  allow  it  to  be  used  with  a  DC  power  pack  or  with  an 
AC  power  converter  pack.  The  power  supply  can  be  included  in  the 
system  case  but  this  will  restrict  its  use  in  that  only  one  power 
source  will  be  available  to  use. 
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Fig  14.  Proposed  Concept  for  the  MFDM 


Fig  15.  Subunit  Layout  for  Proposed  System 


Suaaary 

In  this  chapter,  the  Implementation  of  the  flight-worthy  system 
was  discussed.  The  additional  hardware  needed  to  control  the  auto¬ 
ranging  feature  of  the  vibration  signal  amplifier  was  discussed  along 
with  the  associated  additional  code  in  the  EXECUTIVE  module  required  to 
control  the  gain  of  the  amplifier.  The  modifications  required  to  store 
the  program  in  a  UV-PROM  were  discussed,  and  the  power  requirements  and 
packaging  of  the  MFDM  were  given. 


VIII*  Recommenda t Iona  and  Conclusions 


Introduction 

A  bench  model  of  the  MFDM  currently  exists.  The  Implementation  of 
a  flight -worthy  system  will  require  additional  work.  This  effort  would 
Include  developing  the  auto-ranging  gain  controller,  completing  the 
system  testing,  storing  the  program  in  PROM,  packaging  the  system,  and 
constructing  the  DC  power  supply  required.  This  chapter  presents  some 
recommendations  for  future  development,  and  some  conclusions  drawn  from 
this  effort. 

Recommendations 

Recommendations  for  future  work  include  the  installation  of  the 

IBV11—  A  Instrument  Bus  Interface,  and  the  modification  to  the 

EXECUTIVE  module  necessary  to  control  the  auto-  ranging  amplifier  must 

be  completed.  The  new  code  should  read  the  amplifier's  gain  setting 

from  the  Instrument  Bus  Interface  and  store  a  multiplication  factor 

equivalent  to  the  gain  in  the  global  variable  "ADGAIN".  The  allowable 

••6  22 

range  of  "ADGAIN"  is  from  10  to  10  ,  and  will  insure  proper 

cumulative  damage  factor  computation  (Ref  Chap  VII). 

Further  testing  of  the  system  is  required.  The  testing  should  be 
performed  using  a  random  vibration  signal  as  an  input  to  the  analog 
sensor,  with  the  output  of  the  sensor  going  to  the  MFDM  and  an  audio 
spectrum  analyzer.  A  vibration  table  is  available  at  the  AFFDL  and  can 
be  used  to  generate  the  input  vibration  signals  for  the  system  and  for 
sn  audio  frequency  spectrum  analyzer.  In  addition,  an  analog  data 
acquisition  system  can  be  connected  to  tho  same  vibration  input  signals 
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and  can  be  used  to  check  the  computation  of  the  cumulative  damage 
factor*  This  consists  of  comparing  the  peaks  In  the  outputs.  The 
relative  difference  between  the  two  outputs  should  be  the  sane  for  all 
peaks.  The  spectrum  analyser  can  also  be  used  to  verify  that  the  peaks 
In  the  CDF  table  are  at  the  same  frequency  as  the  Input  spectrum. 

The  system  program  must  be  stored  In  UV-PROM.  This  will  allow  the 
system  to  be  run  without  the  disk  system  that  the  bench  model  used. 

Some  minor  modifications  to  the  code  in  the  FFT  module  will  be  required 
in  order  for  it  to  execute  from  PROM  (Ref  Chap  VII). 

The  system  needs  to  be  assembled  in  an  H0341  card  cage  with  a  DP  I’¬ 
ll  Back  Plane  (Ref  VII).  The  clock  circuit  and  signal  amplifiers  can 
be  included  in  this  card  cage,  resulting  in  a  self-contained  system, 
except  for  batteries 

The  system  casing  must  be  constructed  to  protect  the  modules  from 
damage  in  an  aircraft  environment,  and  to  allow  the  system  to  be 
controlled  with  the  use  of  two  key  switches  (Ref  Chap  VII). 

A  DC  power  supply  must  be  designed.  It  will  have  to  supply  +5,12, 
and  15  volts  DC  at  15.4,  2.7,  and  .007  amps,  respectively,  for  a  period 
of  up  to  20  hours  of  continuous  operation  (Ref  Chap  VII).  The  power 
supply  will  connect  to  the  MFIMl  sy9ten  via  a  power  line  which  plugs 
into  the  back  of  the  MFDM. 

Conclusions 

The  original  design  for  this  system  was  reviewed  with  respect  to  a 
new  set  of  requirements*  The  system  was  then  re-designed  to  meet  the 
new  requirements,  and  a  bench  model  of  the  MFDM  was  constructed* 
Although  time  did  not  allow  testing  of  the  complete  system,  each  of  the 
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■odules  which  stake  up  the  system  was  tested  individually  and  found  to 
be  accurate  within  the  limits  of  a  16  bit  word  size.  The  complete 
computer  system  was  tested  using  sine  wave  generators  to  simulate  the 
output  of  the  signal  amplifiers.  The  sine  wave  generators  produced  a 
known  input  spectrum  which  allowed  the  output  of  the  system  to  be 
compared  to  the  input  signal.  The  frequency  resolution  was  determined 
to  be  half  of  the  frequency  interval  selected.  Due  to  a  lack  of  time, 
a  complete  spectrum  analysis  comparing  the  output  of  the  MFDM  to  a 
known  frequency  spectrum  input  signal  was  not  accomplished. 

The  current  system  demonstrates  the  practicability  of  using  state- 
of-the-art  technology  to  accomplish  the  task  of  real-time  data 
acquisition  and  processing.  In  particular,  the  real-time  data 
acquisition  of  vibration  and  temperature  signals  and  the  real-time 
computation  of  the  cumulative  damage  matrix  from  these  signals  was 
accomplished.  This  system  thus  computes  the  data  required  to  enable 
the  AFFDL  to  produce  a  damping  layer  treatment  for  an  item  being 
tested.  Construction  of  a  flight-worthy  system  can  now  be  accomplished 
following  the  design  discussed  in  Chapter  VIII.  This  investigation  has 
successfully  developed  a  bench  model  system  for  the  real-time 
processing  of  the  cumulative  damage  matrix,  and  has  proposed  a  flight¬ 
worthy  design  which  can  now  be  constructed. 


91 


A  M 


? 

S? 


BIBLIOGRAPHY 


1.  A2S83  Instruction  Manual*  Intech  Corp-  1973. 

2.  ADV11-A,  KWV11-A,  AAVll-A,  DR711  User's  Manual.  Maynard, 
Massachusetts,  Digital  Equipnent  Corporation,  1977. 

3.  Brigham,  E.  Oran.  The  Fast  Fou-ier  Transform.  Englewood  Cliffs. 

New  Jersey:  Prentice-Hall,  Inc.,  1974. 

4.  Cooley,  J.W.,  P.W.  Lewis  and  P.D.  Welch.  "The  Fast  Fourier 
Transform  Algorithm:  Programming  Considerations  in  the  Calculation 
of  Sine,  Cosine,  and  Laplace  Transform."  Journal  of  Sound  and 
Vibration,  12*315-337  (1970). 

5.  Dulaney.  Gerald.  "Designing  the  LSI-11/23  "  Mini-Micro  Systens 
XII  (4): 53-60  (April  1979). 

6.  Girmaldi.  Frank.  Environmental  Standards  For  Computers  and 
Peripherals.  DEC  STD  102  REV  B-  Maynard,  Massachusetts,  Digital 
Equipnent  Corporation,  1976. 

7.  Iftekhar,  Saleem.  Microprocessor  Based  Data  Acquisition 

and  Processing  System.  Master's  Thesis.  Wright-Patterson  Air  Force 
Base,  Ohio:  Air  Force  Institute  of  Technology,  December  1978. 

(AD  A064728) 

8.  Microcomputer  Handbook  (second  edition).  Maynard,  Massachusetts, 
Digital  Equipnent  Corporation,  1977. 

9.  Myers,  J.  Glenford.  Conposite/Structured  Design.  Hew  York:  Van 
Nostrand  Reinhold  Company,  1978. 

10.  Rogers  Lynn,  personal  interviews.  March-November  1979,  Aerospace 
Engineer,  Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson 
Air  Force  Base  Ohio. 

11.  Talmadge,  R.  D.,  personal  interviews  June-November  1979, 

Aerospace  Engineer,  Air  Force  Flight  Dynamics  Laboratory,  Wright- 
Patterson  Air  Force  Base,  Ohio. 

12.  UCSD  (Mini-Micro  Computer)  Pascal.  Version  II. 0  Institute  for 
Information  Systems,  UCSD  Mailcode  C-021,  La  Jolla,  California, 

March  1979. 

13.  VM-112  Integrated  Accelerometer  Short  Form  Brochure.  Vibra 
Metrics  Inc.,  October  1974. 


92 


Appendix  A 

Flow  Chart 8  and  Source  Code 
for  Configuration  Modules 


Fig  16.  Flow  Chart  for  MFDM  (cont) 
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FOPTP.AN  IV  V02.1-11  Hoc'  28-Nov-79  00-00  00 
0001  PROGRAM  MFDM 

C 

C  THIS  MICROPROCESSOR  FLICHT  DAMAGE  MONITOR  (MFDM) 

C  PROCRAM  IS  AN  INTERACTIVE  PROGRAM  WHICH  DOFS  THREE  THINGS 
C  ALLOWS  THE  USER  TO  CONFIGURE  THE  SYSTEM  AS  DESIRED,  ALLOWS 
C  THE  ESER  TO  MANUALLY  CONTROL  DATA  PROCESSING,  ALLOWS  THE  USER 
C  REQUEST  AN  OUTPUT  OF  THE  CDF  MATRIX. 

C 

C  THIS  PROGRAM  IS  STARTED  BY  TYPING  RU  MFDM  IF  UNDER  RT-11, 

C  OR  3Y  TYPING  1000  C  IF  UNDER  SYSTEM  ODT •  DATA  PROCESSING 
C  CAN  BE  STARTED  MANUALLY  BY  TURNING  THE  SYSTEM  HALT /RUN  TO  RUM 
C  AND  THEM  TURNING  THE  DC  POWER  SWITCH  ON. 

C 

C  WHEN  THE  PROCRAM  IS  STARTED  EITHER  OF  THE  FIRST  TWO  WAYS. 

C  THE  USER  WILL  BE  GREETED  BY: 

C 

C  MICROPROCESSOR  FLICHT  DAMAGE  MONITOR 
C  ENTER  S(ET''P,  E(XECUTE,  0(UTPUT,  Q(UIT  ?<CR> 

C 

C  WHEN  THE  USER  ENTERS  THE  FIRST  LETTER  OF  THE  SUBROUTINE  DESIRED 
C  FOLLOWED  BY  A  CARRIAGE  RETURN  (<CR>),,  THAT  SUBROUTINE  WILL  BE 
C  EXECUTED. 

C 


C  SUBROUTINE  SETUP  INITIALIZES  THE  VARIABLES  USED  IN  THE 
C  DIFFERENT  SUBROUTINES  TO  THEIR  INITIAL  VALUES.  THE  DEFAULT 
C  VALUES  ARE  THE  LAST  USED  VALUES.  WITH  THE  VALUES  INITIALIZED  TO 
C  A  FREQUENCY  RANGE  OF  0-3000  HZ  WITH  A  FREQUENCY  INTERVAL  OF 
G  100  HZ  THE  FIRST  TIME  THE  PROGRAM  IS  RUN. 

C 

C  SUBROUTINE  XECUTE  CONTROLS  THE  DATA  ACQUISITION  IN 
C  A  CONTROLLED  ENVIRONMENT  SUCH  AS  A  LABORATORY  WHERE  A  TERMINAL 
C  IS  AVAILABLE- 
C 

C  SUBROUTINE  CDFOUTPUT  PRINTS  OUT  A  CUMULATIVE  DAMACE  MATRIX 
C  FOR  FREQUENCY  VS  TEMPERATURE. 

C 

C  FORTRAN  CO' WON  BLOCKS  (USED  TO  TRANSFER  INFORMATION  TO 
C  AND  FROM  THE  ASSEMBLY  SUBROUTINES). 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


ANSWER 
CDF (1452) 
COMMENTS (40) 

DATE ( 1 3) 
PELTAF 


DURATION 

EXP 


MINUTES 


1USED  TO  HOLD  THE  ANSWER  TO  YES /NO  QUESTIONS 

!THE  CUMULATIVE  DAMAGE  TACTOR  ARRAY 

! AN  ARRAY  OF  ALPHANUMER ICS  WHICH  IS  PRINTED 

!ON  THE  OUTPUT  TABLE- 

!AN  ARRAY  TO  HOLD  THE  DATE 

1THE  SIZE  Or  THE  FREQ  INTERVAL  IN  THE 

!CDF  ARRAY 

!THE  DURATION  OF  A  TEST  IN  HOURS 
ITHE  EXPONENT  FOR  N=2**EXP  WHERE 
! EXP  IS  THE  LOCAL  VARIABLE  CORRESPONDING 
J TO  TUF  GLOBAL  VARIABLE  R  USED  IN  EXEC 
1AND  FOUR  I F, 

ITHE  DURATION  OF  A  TEST  IN  MINUTES 
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C  NFREQ  JTHE  NUMBER  OF  FREQ  INTERVALS  IN  THE 

C  !CDF  ARRAY 

C  POWERFLAG  -  !A  FLAC  USED  IN  THE  POWER-UP  ROUTINE  TO 

C  1INDICATE  THE  FIRST  "AUTOMATIC"  POWER-UP. 

C 

C  VARIABLE  TYPE  DEFINITION 
C 

0002  L0GICAL*1  ANSWER, DATE (13)  COMMENTS (40)  POWERFLAG 

0003  INTEGER* 2  EXP, NFREQ 

0004  INTEGER* 2  MINUTES . INDEX 

0005  REAL *4  DURATION, DELTAF, CDF ( 120, 12) 

C 

C  GLOBAL  VARIABLES  DEFINED 
C 

0006  COMMON  /CDFDATA/CDF. /PTIME /MINUTES , /DELTA/DELTAF 

0007  COMMON  /B/EXP, /INDX/INDEX  /NINTVL /NFREQ, /PWRFLG/POUEP.FLAG 

C 

C  SUSTEM  DEFAULT  ARGUMENTS  INITIALIZED 
C 

0008  DATA  EXP  INDEX, NFREQ  DELTAF/6, 8, 30, ’00. 2/ 

0009  DATA  DATE/'  1'  4'.'-'.'D'  F,','C','-','7','9',4*'  '/ 

0010  DATA  MINUTES/30/ 

0011  DATA  POWERFLAG/ 0/  ’POWERFLAG  IS  INITIALIZED  TO  1 

C  ’FOR  USE  IN  THE  POWER-UP  ROUTINE. 

C 

C  PROGRAM  EXECUTION  STARTS  HERE-  TITLE  AND  A  REQUEST  FOR  THE  USERS 
C  RESPONSE  IS  PRINTED  AT  THE  TERMINAL  FIRST- 
C 


0012 

TYPE  100 

0013 

100 

FOP.’IAT  (2X,  MICROPROCESSOR 

FLIGHT  DA!1AGE  MONITOR') 

0014 

200 

TYPE  300 

0015 

300 

FORMAT  (/2X,  ENTER  S(ETUP, 

E(XECUTE,  OOJTPUT,  0(UIT  -  ?<CR>'$) 

0016 

ACCEPT  400, ANSWER 

0017 

400 

FORMAT  (Al) 

C 

C  THE  USERS  ANSWER  IS  CHECKED  AND  IF  IT  IS  VALID  PROGRAM  CONTROL 
C  PASSES  TO  THE  REQUESTED  ROUTINE- 
C 

0018  IF  ((ANSWER- NE.  S' ) .AND. (ANSWER. NE. 'E' ). AND. (ANSWER. NE. '0' ) .AND. 

1  (ANSWER. NE. 'Q'))  CALL  ERROR  ('INPUT  MUST  BE  AN  S  E  .O.OP,  Q') 

C 

0020  IF  (ANSWER. EQ.  S') 

1CALL  SETUP (DELTAF , NFREQ , DATE . COMMENTS , MINUTES /60 . ) 

C 

0022  IF  (ANSWER. EQ.'E')  CALL  XEC’JTE 

C 

0024  IF  (ANSWER. EQ  '0') 

1CALL  COFOUTPUT (CDF  NFREQ  DELTAF , DATE , CO’PtENTS ) 

C 

0026  IF  (ANSWER. EQ.  Q')  STOP  '  ALL  DONE' 

C 

0028  GO  TO  200 

0029  END 
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FORTRAN  IV  V02.1-11  Frl  23-Nov-79  14:55:18  PAGE  001 

0001  SUBROUTINE  ERROR  (MESSAGE) 

C 

C  ERROR  ROUTINE  WHICH  TYPES  THE  ERROR  MESSAGE  "MESSAGE". 

C  "MESSAGE"  IS  A  STRING  OF  TYPE  LOGICAL* 1  OF  LENGTH  UP  TO  40. 

C 

0002  LOGICAL*!  MESSAGE (40) 

0003  TYPE  100, MESSAGE 

0004  100  FORMAT  (/2X  'ERROR  -  40A1) 

0005  RETURN 

0006  END 
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Fig  18.  Flow  Chart  for  SETUP  (cont) 
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FORTRAN 

0001 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0002 

0003 

000 A 


IV  V02- l-l 1  Tue  27-NOV-79  00-00:00  PACK  001 

SUBROUTINE  SETUP (DELTAF , NFRF.Q,  DATE , COMMENTS , DURATION) 

SETUP  CONTROLS  THE  INITIALIZATION  OF  THE  PROGRAM  VARIABLES- 
IT  IS  AN  INTERACTIVE  ROUTINE  BY  WHICH  THE  USER  CAN  CHANGE  THE 
FREQUENCY  RANGE  AND  INTERVAL  FOR  WHICH  DATA  IS  TO  BE  PROCESSED. 

IN  ADDITION  THE  USER  CAN  CHANCE  THE  RUN  DURATION, DATE  AND  CO'fflENTS. 


THREE  SUBROUTINES  ARE  CALLED  BY  THIS  SUBROUTINE.  THEY  ARE: 

GETFREQUENCY(DELTAF)  -  A  ROUTINE  WHICH  INITIALIZES  THE 
GLOBAL  VARIABLES  USED  IN  THE  ASSEMBLY 
LANGUAGE  ROUTINE  EXEC-  ALSO  VARIABLES 
MAXPREQ  AND  DELTAF  ARE  SET  TO  THE 
MAXIMUM  FREQUENCY  IN  THE  CDF  ARRAY, 

AND  DELTAF  IS  THE  FREQUENCY  INTERVAL 
IN  THE  CDF  ARRAY. 

GETTIME (DURATION)  -  A  ROUTINE  WHICH  READS  THE  DESIRED  RUN 
TIME  FROM  THE  TERMINAL  AND  STORES  IT  IN 
DURATION  (IN  HOURS)  ,A!ID  IN  MINUTES  (IN 
MINUTES). 


VARIABLES  USED  IN  THERE  ROUTINE  ARE 
ANSWER 


COMMENTS (40) 

DATE (1 3) 

DELTAF 

DURATION 

NFKEQ 


-  LOCATION  WHERE  THE  RESULTS  OF  A  Y/N  QUESTION 
IS  STORED. 

-  A  40  BYTE  STRING  TO  CONTAIN  ANY  MESSAGE  THE 
USER  WANTS  TO  BE  PUT  IN  THE  OUTPUT  TABLE- 

-  A  9  BYTE  STRING  TO  HOLD  THE  DATE  WHICH  IS 
PRINTED  ON  THE  OUTPUT  TABLE 

-  THE  FREQUENCY  INTERVAL  IN  THE  CDF  ARRAY. 

-  THE  PLANNED  RUN  TIME  IN  HOURS 

-  THE  NUMBER  OF  FREQUENCY  INTERVALS  IN  THE  CDF 
ARRAY. 


LOGICAL1'  1  ANSWER. DATE (13)  COMMENTS (40) 
INTEGER* 2  NFREQ 
RF.AL*4  DURATION, DELTAF 
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C  THIS  ROUTINE,  BEING  INTERACTIVE.  FIRST  DISPLAYS  THE  CURRENT 
C  VALUE  OF  A  VARIABLE  AND  THEN  ASKS  IF  IT  SHOULD  BE  CHANGED.  THE 
C  USER  RESPONDS  WITH  A  YES  (Y)  OR  NO  (N).  IF  IT  IS  TO  BE  CHANGED 
C  THE  USER  IS  PROMPTED  TO  ENTER  THE  NEW  VALUE. 

C  FOR  THE  DATE  AND  COMMENTS  NO  ERROR  CHECKING  IS  DONE  SO 

C  WHAT  EVER  IS  ENTERED  (PRINTABLE  CHARACTERS  ONLY)  WILL  BE  PRINTED 
C  ON  THE  OUTPUT  TABLE. 


TYPE  20,NFREQ*DELTAF.  DF.LTAF 

FORMAT(///2X, 'FREQUENCY  RANGE:  0  -  '  F6.1.'  HZ'/ 

12X  'DELTAF:  'F5.1,'  HZ') 

TYPE  30, DURATION 

FORMAT (2X  'RUN  DURATION:  '.F5-2,'  HOURS') 

TYPE  40  COMMENTS 

FORMAT (2X, 'COMMENTS :', 10X . 40A 1 ) 

TYPE  50, DATE 

FORMAT (2X, 'DATE: ', 14X, 13A1) 

TYPE  70 

FORMAT (///2X, 'CHANGE  F(RF.Q,  R(UN  DUR  C  ( OMMENTS ' , 

1'  D(ATE,  Q(UIT  ,OR  V(ERIFY  -  ?<CR>'$) 

ACCEPT  80, ANSWER 
FORMAT(Al) 

IF  ((ANSWER. NE- 'F'). AND. (ANSWER. NE. ' R ' ) . AND . (ANSWER. NE. 'C').AND. 
1( ANSWER. NE. ' D ' ) . AN!) . ( ANSWER. NE. 'Q' ) .AND. (ANSWER. NE. 'V'))  CALL 
2  ERROR('  INPUT  MUST  BE  A:  F,R  C  P  Q,OR  V') 


0019 

IF 

(ANSWER  .EQ 

0021 

IF 

(ANSWER  .EO 

0023 

IF 

(ANSWER  .EQ 

0025 

IF 

(ANSWER  .EQ 

0027 

IF 

(ANSWER  .EO 

0029 

IF 

(ANSWER  .EQ 

0031 

GO 

TO  10 

C 

0032 

END 

! CHANGED  FROM  60  TO  10  AFTER  THESIS  WAS 
PRINTED  -  VERIFY  IS  NOW  AUTOMATIC. 
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FORTRAN  IV  V02.1-11  Fri  23-Nov-79  14  54  13 
0001  SUBROUTINE  GETFREQUENCY 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CHANG EF  IS  USED  TO  INITIALIZE  THE  FOLLOWING  VARIABLES 


COMMON 

LOCAL 

DEFINITION 

/B/ 

EXP 

-  THE  INTEGER  WHICH  SATISFIES 

THE  EQ.  N-2**EXP.  WHERE  "N" 

IS  THE  HUMBER  OF  COMPLEX" 

DATA  POINTS  PROCESSED. 

/DELTA/ 

DELTAF 

-  THE  FREQUENCY  INTERVAL  IN 

THE  CDF  ARRAY. 

/INDX/ 

INDEX 

-  AN  INTEGER  INDEX  USED  WITH 

THE  SINE  LOOK-UP  TABLE  IN  THE 

FFT  ROUTINE. 

/NINTVL/ 

NFREQ 

-  THE  NUMBER  OF  FREQUENCY  INTERVALS 
IN  THE  CDF  ARRAY. 

/RTCMODE/  CLKMODE  -  THE  MODE  OF  OPERATION  OF  THE 
PROGRAMMABLE  REAL  TIME  CLOCK, 
POSSIBLE  VALUES  ARE 


CLKMODE  FREQUENCY 
(DECfflAL) 

1 l  1  MHZ 
19  100  KHZ 
27  10  KHZ 


/RTCRATE/  CLKRATE  -  THE  FREQUENCY  DIVISOR  FOR  THE  CLKMODE 


VARIABLE  TYPE  DEFINITION 


0002  INTEGER*2  EXP  INDEX. INTERVAL, ENTRY. NFREQ 

0003  INTEGER* 2  CLKMODE  CLKRATE 

0004  REALM  DELTAF 


C 

C  GLOBAL  VARIABLES 
C 


0005  COMMON  /B/EXP. /INDX/INDEX. /NtNTVL/NFREQ. /DELTA/DELTAF 

0006  COMMON  /RTCMODE/CLKMODE  /RTCRATE /CLKRATE 

C 

C  DEFAULT  REAL-TIME  CLOCK  SETTINGS 
C 

0007  DATA  CLKMODE  / 11/, CLKRATE  /-1 56/ 


0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 


0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 


C 

C 

C  THE  POSSIBLE  COMBINATIONS  OF  THE  FREQUENCY  RANGE  ANT)  FREQUENCY 
C  INTERVAL  ARE  DISPLAYED  FIRST.  THE  USER  IS  THEN  PROMPTED  TO  ENTER 
C  THE  INTEGER  (1-6)  CORRESPONDING  TO  THE  COMBINATION  DESIRED.  THE 
C  INTEGER  IS  USED  IN  A  COMPUTED  GO  TO  STATEMENT  TO  BRANCH  TO  THE 
C  APPROPRIATE  LOCATOIN  FOR  VARIABLE  INITIALIZATION. 

C 

C 

100  TYPE  200 

200  FORMAT (//2X, 'FREQ.  RANGE  (KHZ)  DELTAF  (HZ) 9X, 'ENTER'/) 

TYPE  300, 6. .100. 2. 1,6. ,50.1. 2, 3. ,100.2, 3. 3.  50*1.4  3. .25. 0.5, 
91.5.100.2,6  1.5.50.1,7, 1.5,25.0,8 

300  FORMAT ( 8 ( 8X  F4-1.12X  F5. 1, 15X. II/)) 

TYPE  400 

400  F0RMAT(//1X. 'ENTER  A  NUMBER  BETWEEN  1  &  8  CORRESPONDING  TO  THE', 
1'  DESIRED  SETUP  ?<CR>'$) 

ACCEPT  500, ENTRY 

500  FORMAT (II) 

GO  TO  (1  2, 3, 4, 5  6,7.8)  ENTRY  ! IF  ENRTY>8  OR  <1  CONTROL 

C  ! FALLS  THROUGH. 

CALL  ERROR( 'INPUT  MUST  BE  AN  INTEGER  BETWEEN  1  AND  8') 

GO  TO  100 
C 

C  VARIABLE  INITIALIZATION  STARTS  HERE 
C 

1  EXP-7 
INDEX-4 
NFREQ-60 
DELTAF- 100. 2 
CLKRATE— 78 
GO  TO  600 

2  EXP-S 
INDEX-2 
NFREQ-120 
DELTAF-50. 1 
CLKRATE— 78 
GO  TO  600 

3  EXP-6 
INDEX-8 
NFREQ-30 
DELTAF- 100. 2 
CLKRATE— 156 
GO  TO  600 

4  EXP-7 
INDEX-4 
NFREQ-60 
DELTAF-50. 1 
CLKRATE— 156 
GO  TO  600 
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0043 

5 

EXP-8 

0044 

INDEX-2 

0045 

NFREQ-120 

0046 

DELTAF-25. 

0047 

CLKRATE— 156 

0048 

GO  TO  600 

0049 

6 

EXP-5 

0050 

INDEX- 15 

0051 

NFREQ-1 5 

0052 

DELTAF-100.2 

0053 

CLKRATE— 312 

0054 

GO  TO  600 

0055 

7 

EXP-6 

0056 

INDEX-8 

0057 

NFREQ-30 

0058 

DELTAF-50. 1 

0059 

CLKRATE—  312 

0060 

GO  TO  600 

0061 

8 

EXP-7 

0062 

INDEX-4 

0063 

NFREQ-60 

0064 

DELTAF-25. 

0065 

CLKRATE— 312 

0066 

600 

RETURN 

0067 

END 

107 


Fig  20.  Flov  Chart  for  GETTIME 
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FORTRAN  IV  V02.1-U  Frl  23-Nov-79  14:  56' 19 
0001  SUBROUTINE  GETTIME  (  DURATION  ) 

C 

C  TIME  GETS  THE  DESIRED  DATA  PROCESSING  DURATION  FROM  THE  USER 

C  AND  STORES  IT  IN  THE  GLOBAL  VARIABLE  RTIME  (USING  CO!EtON). 

C  RTIME  IS  AN  INTEGER  (LOCAL  VARIABLE  MINUTES),  AND  DURATION  IS  THE 
C  REAL  VARIABLE  EQUIVALENT  TO  MINUTES. 

C 

C  VARIABLES  DEFINED  HERE 

V 

0002  INTEGER*2  MINUTES 

0003  REAL*4  DURATION 

0004  COMMON  /PTIME/MINUTES 

C 

C  REQUEST  USER  TO  ENTER  A  NF.U  TIME 
C 


0005 

100 

TYPE  200 

0006 

200 

F0RMAT(/2X, 'ENTER  DURATION  OF  TEST  (HOURS)  EXAMPLE  15  25  <CR>'$) 

0007 

ACCEPT  300, DURATION 

0008 

300 

FORMAT (F5.0) 

0009 

IF  (DURATION. LE  0- ) 

1CALL  ERROR ('DURATION  MUST  BE  IN  HOURS  AND  >  .01666') 

0011 

IF  (ERR)  GO  TO  100 

0013 

MINUTES “DURATION *60. 

0014 

RETURN 

0015 

END 
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21.  Flow  Chart  for  GETDATE 
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FORTRAN  IV  V02. 1-11  Tue  27-Hov-79  00:00t00 
0001  SUBROUTINE  GETDATE (  DATE  ) 

C 

C  DATE  READS  IN  THE  USER-SUPPLIED  DATE  AND  STORES  IT  IN  THE  STRING 
C  DATE.  UP  TO  13  CHARACTERS  ARE  ALLOWED •  DATE  IS  ONLY  USED  IN 
C  SUBROUTINE  GRAPH  TO  PRINT  OUT  THE  DATE. 

C 

C 


0002 

0003 

C 

0004 

0005 

100 

0006 

0007 

200 

0008 

0009 

LOG I CAL *1  DATE (13) 

INTEGER* 2  I 

TYPE  100 

FORMAT (/2X  'ENTER  THE  NEW  DATE  EX. 
ACCEPT  200, (DATE(I) , 1=1 . 13) 

FORMAT (13A1) 

RETURN 

END 


13-AUG-79<CR>'$) 


A 
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Fig  22.  Flow  Chart  for  GETCOMMENT 
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FORTRAN  IV  V02.1-11  Tue  27-Nov-79  00:00  00 
0001  SUBROUTINE  GETCOMMENT (  COMMENTS  ) 

C 

C  COMM  READS  IN  USERS  COMMENTS,  UP  TO  40  CHARACTERS 
C  AND  STORES  THEM  IN  THE  ARRAY  COMMENTS.  NO  ERROR  CHECKING 
C  IS  DONE.  WHAT  EVER  IS  ENTERED  WILL  BE  PRINTED  ON  THE 
C  OUTPUT  TABLE. 

C 


0002 

LOGICAL* 1  COMMENTS (40) 

0003 

C 

INTEGER* 2  I 

0004 

C 

DO  100, 1=1, 40 

0005 

COMMENTS (I )='  ' 

0006 

100 

CONTINUE 

0007 

TYPE  200 

0008 

200 

FORMAT (/2X, 'ENTER  COMMENTS  -  40  CHARACTERS  MAX  <CR>'$) 

0009 

ACCEPT  300, (COMMENTS (I) .1=1.40) 

0010 

300 

FORMAT (:40A1) 

0011 

RETURN 

0012 

END 

START 

N - 


Fig  23.  Flow  Chart  for  CDFOUTPUT 
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0001  SUBROUTINE  CDFOUTPUT (CDF . NFREQ  DELTAF, DATE, COMMENTS) 


OUTPUT  IS  USED  TO  PRINT  OUT  THE  CDF  ARRAY.  THE  STRUCTURE  OF 
THE  CDF  ARRAY  IS: 


CDF(l, 1)- 


{NFREQ  NUMBER  OF  FREQUENCIES  DELTAF  APART 


CDF(NFREO, 1) 
CDF(1, 2) 
CDF(2. 2)- 


: (0-MAXFREQ) 

; START  OF  TEMP  INTERVAL  #2 


CDF (NFREQ. 2) 


CDF ( 1 , 12) 
CDF(2, 12) 


CDF (120  12) 


THE  VARIABLES  USED  IN  THIS  SUBROUTINE  ARE: 


CDF (120  12) 
CONTENTS (40) 

DATE (13) 

DELTAF 

FREQPAGE 
I  J  K,L 
LARGE 


LOW, HIGH 


HICHFREQ 


;THE  CUMULATIVE  DAMAGE  FACTOR  ARRAY 
A  STRING  FOR  USERS  COTP1EUTS 
THE  USER  WANTS  ON  THE  OUTPUT 
A  STRING  TO  HOLD  THE  DATE  WHICH  WILL 
BE  PRINTED  ON  THE  OUTPUT- 
THE  FREQUENCY  NFREQ  BETWEEN  ELEMENTS 
OF  THE  CDF  ARRAY 

THE  FREQUENCY  RANGE  PER  PACE  OF  OUTPUT 
GENERAL  PURPOSE  INDEXES 
A  THREE  CHARACTER  ALPHA  CONSTANT  WHICH 
WHEN  SENT  TO  THE  PRINTER,  CAUSES  80 
CHARACTERS  PER  LINE  TO  BE  PRINTED 
TWO  INDEXING  LIMITS  USED  TO  CONTROL 
THE  TEMPERATURE  SCALE  PRINTING 
LOW-1—50  F;HIGll-7-100  F 
LOW-6-75  F’HIGH-12-225  DEG  F. 

INDEX  LIMITS  TO  CONTROL  THE  PRINTING 
OF  THE  FREQUENCY  IN  THE  TABLE 


•-in-'  'i1* 


..vmxiff'*-1';  :->**srV *****»*'. */*"  nr. » -at-  r.>  ...»  ., 


0002 

0003 

0004 

0005 

0006 


0007 


0008 

0009 


0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


L0WTF.MP ,  HIGHTEMP 
LTSYM 


MAXFREQ 

NFREQ 

NPAGE 

/RTIME/RUINUTES 

RTSYM 


•.INDEX  LIMITS  TO  CONTROL  THE  PRINTING 
;0F  THE  TEMPERATURE  IN  THE  TABLE 
;THE  SYMBOL  TO  THE  LEFT  OF  THE  FIRST 
; TEMPERATURE,  IN  THE  TABLE  HEADING.  IT  IS 
•.EITHER  A  SPACE  OR  <  (LESS  THAN). 

;THF.  MAXIMUM  FREQUENCY  IN  THE  CDF  ARRAY 
;THE  NUMBER  OF  FREQUENCY  COMPONENTS  IN 
;THE  CDF  ARRAY 

;THE  NUMBER  OF  PAGES  REQUIRED  TO  PRINT 
;THE  TABLE  FOR  ONE  TEMPERATURE  RANCF 
‘.THE  ACTUAL  DURATION  OF  THE  TEST  RUN 
;THE  SYMBOL  PRINTED  TO  THE  LEFT  OF  THE 
;LAST  TEMPERTURF,  ON  THE  TABLE  HEADING. 
•IT  IS  EITHER  A  SPACE  OR  A  >  (GREATER 
;THAN) 


TYPE  DEFINITION 

LOG I CAL* 1  DATE ( 1 3 ) , COMMENTS ( 4  0 ) . LARC  E ( 4 ) , SMALL ( 4 ) , LT5  YM , RTS  YM 
INTEGER* 2  NFREQ  PAGE, LOWTEMP. HICHTEMP, I, J.K  L.M  N 
INTEGER*2  LOW  HIGH  NPAGE  MINUTES 
REAL*4  DURATION. DELTAF, FREQ, CDF(120  12) 

REAL*4  MAXFREQ  HIGHFRF.Q  FREQPAGF. 

GLOBAL  VARIABLE 

COMMON  /RTIME/MINUTES 

SMALL  INITIALIZED  FOR  H14  PRINTER  -  12  CHAR/INCH  PRINTING. 

LARGE  -  80  CHAR/INCH. 

DATA  SMALL/27. 117.24,0/ 

DATA  LAP.GE/27  117  4,0/ 

PROGRAM  EXECUTATION  STARTS  HERE. 

DURATION-MINUTES / 60 . 

NPAGE- (NFREQ /4 1 ) +1 

LOW  -  1 

HIGH-7 

FRF.QPAGE-4 1 .  *DF,LTAF 
LOWTEMP— 50 
HIGHTEMP- 100 
PAGE-0 
LTSYM- '<' 

RTSYM- '  ' 

MAXFREQ-NFREQ*DEI,TAF 


»  • 


c 

C  THE  CDF  TABLE  FOR  HALF  OF  THE  TEMPERATURE  RANGE  IS  HOW  PRINTED  OUT 
C 


0021 

DO  150,1-1  2 

0022 

FREQ-DELTAF 

0023 

HIGHFREQ-O. 

0024 

K-0 

0025 

TYPE  80  LARGE 

C 

C  PRINTING  OF  EACH  PAGE  STARTS  HERE. 


0026 

C 

DO  140.J-1  MPAGE 

0027 

HIGHFREQ-HIG'irREQJ-FREQPAGE 

0028 

IF(HIGHFP.EQ.CT.MAXFKEQ)  HIGHFREQ-HAXFREQ 

0030 

TYPE  10 

0031 

10 

FORMAT (IH1////) 

0032 

PACE-PAGK+l 

0033 

TYPE  20, PAGE 

0034 

20 

FORMAT (24X  'CUMULATIVE  DAMAGE  MATRIX' 13X, 'PAGE  ',12/) 

0035 

TYPE  30, DATE 

0036 

30 

FORMAT (IX  'DATE'  6X  13A1) 

0037 

TYPE  35, DURATION 

0038 

35 

FORMAT (IX  DUR AT ION ' , 2X , F 5 . 2  '  HOURS') 

0039 

TYPE  40, COMMENTS 

0040 

40 

FORMAT (IX, 'COMMENTS'  2X  40A1/) 

0041 

TYPE  50 

0042 

50 

FORMAT ( IX .  FREQ  (HZ)'  18X  'TEMPERATURE  (DEG  F)') 

0043 

TYPE  70.LTSYM,  (L,L-T,OWTEMP  HIGHTEMP-25  25) ,  RTSYM .  HTGHTEMP 

0044 

70 

FORMAT (9X. A  1,5(13. 7X),I3.6X  A1,I3) 

0045 

TYPE  80  SMALL 

0046 

80 

FORMAT ('+' ,3A1) 

0047 

TYPE  90 

0048 

90 

r* 

FORMAT ( 7X, 86 ('-' )) 

Vi 

C  CDF  DATA  IS  PRINTED  BY  THE  FOLLOWING  CODE 

n 

0049 

\j 

100 

IF  (FREQ.GT  HICHFREQ)  GO  TO  130 

0051 

K-K+l 

0052 

TYPE  110, FREQ, (CDF(K. N) , N-LOW  HIGH) 

0053 

110 

FORMAT( IX, F5. 0 . IX. ' | '  7 ( IX  E10.3  IX) ,  ' | ' ) 

0054 

FREQ-FREQ+DELTAF 

0055 

GO  TO  100 

0056 

130 

CONTINUE 

C  END  OF  PRINTING  TABLE;  NOW  PRINT  TEMP  AXIS. 

O 

0057 

v» 

TYPE  90 

0058 

TYPE  80  LARGE 

0059 

TYPE  70 . LTSYM, (L  L-LOWTEMP  HIGHTEMP-25, 25) .RTSYM  1IIGHTEMP 

0060 

TYPE  50 

0061 

K-K-l 

0062 

FREQ-FREQ-DELTAF 

0063 

140 

CONTINUE 

117 


0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 


C 

C  SET  VARIABLES  TO  PRINT  HIGH  TEMPERATURE  RANGE. 
C 

LOW-6 
HIGH-12 
LOWTEMP-75 
HIGHTEMP-225 
LTSYM-*  ' 

RTSYM-'>' 

150  CONTINUE 
RETURN 
END 
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Appendix  B 

Flow  Charts  and  Source  Code  for 
Data  Processing  and  Execution  Modules 
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;  24.  Flow  Chart  for  EXECUTIVE 
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MACRO  V03.01  27-NOV-79  PAGE  1 


.TITLE  EXECUTIVE 

. SBTTL  INTRODUCTION  TO  CONTROL  EXECUTIVE 


1.  INTRODUCTION  TO  CONTROL  EXECUTIVE: 


THIS  PART  OF  THE  CONTROL  EXECUTIVE  GETS  CONTROL  OF 
THE  SYSTEM  ONCr  POWER  HAS  BEEN  SNITCHED  ON  AND  THE 
RUN/STANDBY  SNITCH  HAS  BEEN  PLACED  IN  RUN  POSITION.  THE 
EXECUTIVE  WILL  RELINQUISH  CONTROL  IF  THE  POWER 
HAS  BEEN  SWITCHED  OFF,  THE  SYSTEM  HAS  BEEN  SWITCHED 
TO  STANDBY  MODE.  OR  THE  RUN  TIME  HAS  EXPIRED. 


NOTE:  THE  GLOBL  DATA  STORAGE  AREAS  ARE  DEFINED  BELOW 


. GLOBL  B  BITCOUNT. COUNT. INTCELNUM 

.GLOBL  N.N2.NM.NINTVL.TFLG 

.GLOBL  EXEC  BEGIN, RTCRATK . RTCMODE  TIMEINT 

•  GLOBL  ADDATA  PWRSPC ,  CDFDATA 

.GLOBL  INDX.RTIME  PTIME. ER, ADCAIN 


;NOTE:  THE  GLOBAL  SUBROUTINES  ARE  DEFINED  BELOW- 

t 

.GLOBL  FFT  POWERSPECTRUM.DAMACEFACTOR 

•.NOTE-  THE  REGISTERS  ARE  DEFINED  BELOW 

.MCALL  .REGDEF 
•REGDEF 
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DEFINITION  OF  DEVICES 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 


•SBTTL  DEFINITION  OF  DEVICES 
; 2.  DEFINITION  OF  DEVICES: 


THE  DEVICES  THAT  ARE  CALLED  IN  BY  THE  EXECUTIVE 
ARE  DEFINED  HERE. 

ANALOG  TO  DIGITAL  CONVERTOR  (ADC) 


ADCSR-1 70400  -ADC  CONTROL  AND  STATUS 

; REGISTER. 

ADBFR-1 70402  ;ADC  DATA  BUFFER  REGISTER 


.PSECT  ADGAIN , RN , D ,  GBL .  REL . OVR 

ADGAIN :  0  ;ADGAIN  IS  FOR  FUTURE 

•  PSECT  '.IMPLEMENTATION.  IT  WILL  CONTAIN 

;THE  GAIN  OF  THE  VIBRATION  SIGNAL 
^AMPLIFIER. 

;  REAL  TIME  CLOCK  (RTC) 


RTCSR-1 70420  ;RTC  CONTROL  AND  STATUS  REG 

RTBFR=1 70422  ;RTC  DATA  BUFFER  REGISTER 

> 

-.NOTE: 


;  BEFORE  THESE  POINTS  ARE  PROCESSED,  THE  PROCESS 

LOOP  COUNTER  (PRCLOP)  IS  DEFINED 
;  AND  INITIALIZED  BELOW • 

y 

PRCLOP'  .WORD  II  ;THIS  COUNTER  IS  USED  TO 

.-DETERMINE  HOW  MANY  POWER 
:SPECTRUH  SETS  HAVE  TO  BE 
ADDED  TOGETHER. 

; 12 (OCTAL)  INTERVALS  0-11. 
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INTERRUPT  SERVICE  ROUTINE- 

1  .SBTTl  INTERRUPT  SERVICE  ROUTINE. 

2  ; 3.  INTERRUPT  SERVICE  ROUTINE; 

3  ; - - 

4  TIME  INTERRUPT  ROUTINE 

5 


6 

;AN  EXTERNAL 

INTERRUPT  (A  TIME 

SIGNAL  EVERY  MINUTE) 

7 

; CAUSES  A  BRANCH  TO  THIS  ROUTINE  WHICH  INCREMENTS 

8 

;THE  MINUTES 

COUNTER  -  RTIME. 

y 

10 

TIME IN?:  INC 

TIME 

•.INCREMENT  TIME  BY  ONE  COUNT 

11 

CMP 

#74. TIME 

;WHEN  TIME  EQUALS  60  (BASE  10) 

12 

; INCREMENT  RTIME  1 . IF  THE  SYSTEM 

13 

CLOCK  IS  SET  TO  6011Z  RTIME 

14 

WILL  BE  IN  HOURS.  IF  IT  IS 

15 

SET  TO  1HZ  RTIME  WILL  BE  IN  MIN 

16 

BNE 

TIM 

17 

INC 

RTIME 

18 

CLR 

TIME 

•RESET  TIME  TO  0. 

19 

TIM  RTI 

20 

21 

22 

\ 

23 

» 

24 

;  4.  REAL  TIME 

ERROR  ROUTINE- 

25 

; - 

26 

» 

27 

;THIS  ERROR  ROUTINE  IS  OES IONS') 

TO  HANDLE  REAL-TIME  ERRORS. 

28 

;TF  THE  ANALOG  DATA  IS  BEING  SAMPLED  TOO  FAST,  THIS  ROUTINE 

29 

; IS  CALLED. 

A  LOCAL  VARIABLE  ” 

ER"  IS  A  COUNTER  WHICH 

30 

INDICATES  THE  NUMBER  OF  TIMES 

THIS  ROUTINE  HAS  BEEN  CALLED. 

31 

;IT  IS  FOR  MAINTENCE  PURPOSES  ONLY.  "ER"  CAN  BE  EXAMINED 

32 

; USING  ODT . 

33 

J 

34 

ERROR;  MOV 

#0  <3#RTCSR 

;STOP  RTC 

35 

MOV 

0#/,^BFR  (R3) 

; CL EAR  A/D  DONE  FLAG. 

36 

BIG 

#100200, P#ADCSR 

;CLF,AR  A/D  ERROR  FLAG. 

37 

INC 

ER 

INCREMENT  ERROR  COUNTER. 

38 

BR 

ERRTN 

;  RETURN 

39 

TMPERR:  MOV 

(R 1 ) , R3  REMOVE 

BAD  DATA  &  CLEAR  DONE  FLAG. 

40 

•DONE  FLAG. 

41 

BIC 

#100000, (RO) 

-.CLEAR  ERROR  FLAG. 

42 

INC 

(RO) 

; RESTART  CONVERSION. 

43 

BR 

LOOP 
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FPASS  —FIRST  PASS  THROUGH  THE  EXECUTIVE 

1  .SBTTL  FPASS  —FIRST  PASS  THROUGH  THE  EXECUTIVE 

2  .SBTTL  PROCESSING  CONTROL 


3 

4 

C 

•5.  PROCESSING 

CONTROL: 

6 

THE  TEMPERATURE  SIGNAL  IS  SAMPLED  FIRST.  THE  DIGITIZED 

7 

,VALUE  IS  SUMMED  IN  A  STORAGE 

VARIABLE  "TSUM". 

8 

THEN  A  SET 

OF  N  FREQUENCY  DATA  POINTS  ARE  DIGITIZED. 

9 

THE  FOURIE  TRANSFORM,  POWER  SPECTRUM,  ARE  THEN  COMPUTED.  THIS 

10 

PROCESS  -  ACQUIRE  TEMPERATURE 

DATA  FREQUENCY  DATA.  COMPUTE 

11 

THE  FFT  AND  THE  POWER  SPECTRUM 

-  IS  PERFORM'D)  TEN  TIMES  FORMING 

12 

THE  AVERAGE  POWER  SPECTRUM  AND 

THE  AVERAGE  TEMPERATURE. 

13 

THESE  AVERAGE  VALUES  ARE  THEN  USED  TO  COMPUTE  THE  CUMULATIVE 

14 

DAMAGE  FACTOR 

15 

16  EXEC:  JSR 

PC . IHITAL 

•.INITIALIZE  VARIABLES 

17  BEGIN:  MOV 

@#ADRFR  RO 

;CLF.AR  AD  BUFFER  AND  DONE  FLG. 

18 

MOV 

0ADCSR.RO 

M10VE  ADDRESS  OF  ADC SR  AMD  ADBFR 

19 

MOV 

02400, (RO) 

;SET  TEMP  CONVERSION  CH=05 

20 

AND  START  TEMP  CONVERSION. 

21 

MOV 

0ADBFR  R1 

;INTO  REGISTERS  FOR  FASTER 

22 

; INSTRUCTION  EXECUTION  IN  THE 

23 

;TIME  CRITICAL  LOOP  (LOOP2) 

24 

; WHICH  FOLLOWS. 

25 

26 

INC 

(RO) 

; START  TEMP  ACQUISITION. 

27  LOOP  TSTB 

(RO) 

{WAIT  FOR  TEMPERATURE  CONVERSION 

28 

BPL 

LOOP 

TO  BE  COMPLETED 

29 

BIT 

#100000, (RO) 

;CK  FOR  A/D  ERROR 

30 

BHE 

TMPERR 

31 

MOV 

(Rl) ,R3 

{GET  TEMP. 

32 

SUB 

#4000, R3 

{REMOVE  BIAS 

33 

ADD 

R3.TSUM 

{ADD  TEMP  TO  FORM  A  AVERAGE  TIMES 

34 

12 (OCTAL) 

35 

MOV 

#7040, (RO) 

{LOAD  ADCSR-  CH- 16. INTERRUPT 

36 

{DISABLED 

37 

NOTE: 

38 

***************************************************************** 

39 

40 

THE  ADGAIN  OF 

THE  VIBRATION  SIGNAL  AMPLIFIER  WILL  BE  DETERMINED 

41 

HERE.  THIS  HAS  NOT  YET  BEEN  IMPLEMENTED  SO  THE  VARIABLE 

42  ; 

"ADGAIN"  IS  SET  EQUAL  TO  1,  SO 

THAT  IT  HAS  NO  EFFECT  ON  THE 

43 

POWER  SPECTRUM 

CALCULATION  MADE  IN  POWERS PECTRUM 

44 

45 

MOV 

#1 , ADCAIN 

46 

47 

48 

49 

MOV 

ADDATA  R3 

{LOAR  R3  WITH  LOC  FOR  DATA 

50 

MOV 

N,R2 

{LOAD  R2  WITH  THE  #  OF  DATA  POINTS 

51 

MOV 

RTCRATE ,  @#RTRFP, 

•SET  RT  COUNTER  TO  -RTRATE 
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PROCESSING  CONTROL 


1 

MTPS 

#340 

;DISABLE  ALL  INTERRUPTS 

2 

{(TIME  INTERRUPTS) 

3 

A 

ERRTN  MOV 

RTCMODF.,@#RTCSR  ; START  RTC  MODE  1 .  1  OK  HZ. 

5 

f. 

* 

;NOTE: 

7 

8 

THE  LOOP  STARTING  AT  "LOOP:" 

AND  GOING  THROUGH  SOB  R 2  LOOP  IS 

9 

;A  TIME  CRITICAL  LOOP.  IT  IS 

THE  A/D  CONVERSION  ROUTINE  WHICH  IS 

10 

•FASTER  THAN 

AN  INTERRUPT  ROUTINE. 

11 

;  THE  EXECUTION  TIME  OF  THE 

LOOP  MUST  BE  LESS  THAN  THE  AD  CONVERSION 

12 

;TIME  OR  THE 

SAMPLING  PERIOD, 

WHICHEVER  IS  SMALLER. 

13 

• 

14 

15 

LOOP2:  TSTB 

(RO)  {WAIT 

FOR  AD  CONVERSION 

16 

BPL 

LOOP  2 

;TO  FINISH. 

17 

BIT 

#100000, (RO) 

{DONE:  CK  FOR  A  CONVERSION  ERROR 

18 

BNE 

ERROR 

•.ERROR  GO  TO  ERROR  ROUTINE 

19 

MOV 

(Rl) , (R3) 

;MOVE  DATA  INTO  DATA  ARRAY 

20 

SUB 

#4000, (R3)+ 

; REMOVE  4000  BIAS  AND  INC  R3 

21 

SOB 

R2.LOOP2 

;GO  GET  MORE  DATE  UNTILL  ALL 

22 

{DATA  POINTS  HAVE  BEEN  PROCESSED 

23 

J 

24 

1 

25 

END  OF  TIME 

CRITICAL  SECTION 

26 

• 

27 

CLR 

@#RTCSR 

{STOP  RTC- THE  RTC  MUST  BE  STOPPED 

28 

SO  THE  PROPER  START-UP  SEQUENCE 

29 

{STARTINC  AT  LOOP2  IS  EXECUTED. 

30 

MTPS 

#0 

ENABLE  INTERRUPTS 

31 

JSR 

PC.FFT 

32 

;NOTE: 

33 

• 

34 

;  AT  THIS 

POINT,  THE  VALUE  OF  SFACTR  WILL  SHOW 

35 

THE  NUMBER  OF  BITS  THE 

FOURIER  TRANSFORMED  ARRAY 

36 

HAS  BEEN  SHIFTED  1-RIT 

TO  THE  RIGHT. 

37 

i 

38 

» 

39 

JSR 

PC, POWERS PECTRUM 

40 

*1 


V 

t 
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PROCESSING  CONTROL 

* 1  l  ; 

2  ; 

3  ; 


4 

TST 

PRCLOP  ;IS  PRCLOP-O 

5 

BEQ 

TEMP  ;IF  NOT  GO  BACK 

6 

DEC 

PRCLOP  ; COMPUTE  ANOTHER  SET  OF 

7 

JMP 

BEGIN  ; POWER  SPECTRUM. 

8 

9 

10 

MOTH: 

11 

12 

BEFORE 

DF  IS  CALCULATED,  THE 

13 

,THE  TEMPERATURE  INTERVAL  IS 

14 

DETERMINED,  AN!)  "TFLG"  IS 

15 

LOADED 

WITH  THE  APPROPRIATE 

16 

VALUE. 

17 

18 

19 

TEMP: 

MOV 

#T INTVL, Rl 

20 

MOV 

#1 3, R2 

21 

MOV 

*1  TFLG 

22 

IN3: 

CMP 

TSIL1 ,  (R1H 

23 

BLE 

IN5 

24 

INC 

TFLG 

25 

SOB 

R2.IN3 

26 

IN5: 

3R 

DAMCE 

27 

28  ;  INTVL 

1 

2  3  4 

29  TINTVL 

164413, 

170452,174561  000737 

30 

TEMP: 

-37.5 

-12.5  12.5  37.5 

31 

INTVL 

5 

6  7  8 

32 

005166, 

011466,016034,022453 

33 

TEMp: 

62-  5 

87.5  112.5  137  5 

34 

INTVL: 

9 

10  11  12 

35 

027022, 

033371  040010.044427 

36 

TEMP: 

162.  5 

187.5  212.5  237.5 

PROCESSING  CONTROL 


1  : 

2  ; PROCESSING  CONTROL  CONTINUED . 

3  ; 

4 

5  DAMGE :  JSR  PC  DAMAGEFACTOR  COMPUTE  DF 

6 

7  ,-NOTE:  ONCE  THE  DF  IS  COMPUTED,  THE  WHOLE  LOOP 

8  • -  IS  REPEATED  AGAIN.  THE  FOLLOWING  CLEARS 

9  ;  THE  REQUISITE  VARIABLES  AND  ARRAYS. 


10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 

21  CLEAR:  CLR  (R0)+ 

22  SOB  Rl, CLEAR 

23  ; 

24  ;PTIME  IS  A  GLOBAL  VARIABLE  INITIALIZED  IN 

25  .-SUBROUTINE  SETUP  IN  MFDM.  IT  IS  THE  TIME  IN  MINUTES. 

26  ;FOR  WHICH  THE  TEST  SHOULD  RUN.  P.TIUE  IS  A  GLOBAL 

27  ; VARIABLE  WHICH  IS  THE  ACTUAL  RUN  TIME  RETURNED 

28  TO  MFDM  . 

29  ;RTIME  IS  NOW  CHECKED  AGAINST  PTIME  AND 

30  ;IF  RTIME  IS  LESS  THAN  PTIME  PROCESSING  CONTINUES 

31  ;EI.SE  CONTROL  IS  RETURNED  TO  MFDM  OR  THE  SYSTEM 

32  ; HALTS  (IF  RUNNING  WITH  OUT  A  TERMINAL). 

33  ; 

34  CMP  PTIME. RTIME  ;IF  PLANNED  RUN  TIME 

35  •  BCT  BEGIN  ;>  RUN  TIME  CONTINUE 

36  RTS  PC  ;ELSE  RETURN 

37 

38 

39  END  OF  PROCESSING  CONTROL 


CLR  TSUM 

MOV  (?]  1  PRCLOP 

MOV  #PWRSPC  RO  ;THE  POWER  SPECTRUM 

IS  NOW  CLEARED. 

MOV  NINTVL.Rl  ; LOAD  RO  WITH  THE 

STARTING  ADDRESS  AND 
:R1  WITH  ITS  SIZE- 
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DATA  STRUCTURE  AND  VARIABLES  USED 


1  .SBTTL  DATA  STRUCTURE  AND  VARIABLES  USED 


6.  DATA  STRUCTURE  AND  VARIABLES  USED- 


,  THE  DATA  STRUCTURE  WHICH  IS  USED  TO  STORE  ACQUIRED 
DATA  AND  DIFFERENT  VARIABLES  USED  ARE  GIVEN  BELOW: 


.PSECT  NINTVL  RW,D,CBL,REL  OVR 
NINTVL  .WORD  0 

.PSECT  TFLC,RW,D,GBL,REL,OVR 
TFLG :  .WORD  0 

.PSECT  ADDATA  RW  D,GBL,REL,OVR 
DATALOC  .REPT  AOO 
•WORD  0 
.  ENDM 

ADDATA-  DATALOC 

.PSECT  PTIME  RW ,  D , GBL , REL , OVR 
PTIKE  .WORD  0 


•  -  - 


INITIALIZE 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


.SBTTL  INITIALIZE 
.PSECT 


1.  DESCRIPTION 


THIS  SUBROUTINE  IS  USED  TO  INITIALIZE  VARIOUS 
; PARAMETERS  USED  IN  EXEC  AND  FOURIE.  TWO  PARAMETERS 
r(B  AND  1NDX)  ARE  INITIALIZED  BY  SUBROUTINE  SETUP  WHICH 
;IS  CALLED  BY  MFDM.  THESE  TWO  VARIABLES  ARE  USED  TO 
; INITIALIZE  THE  REST  OF  THE  VARIABLES  IN  THIS  ROUTINE. 

;THE  VARIABLES  WHICH  ARE  INITIALIZED  IN  THIS 
{SUBROUTINE  ARE: 

{BITCOUNT  A  COUNTER  USED  IN  THE  FOURIE 

BIT  INVERSION  ROUTINE  IT  IS 
{  EQUAL  TO  THE  NUMBER  OF  BITS 

IN  THE  NUMBER  -  (2**B)-1. 

OR  BITCOUNT-B-1. 

COUNT  A  COUNTER  USED  IN  FOURIE  AND 

INITIALIZED  TO  HALF  THE  NUMBER 
OF  COMPLEX  DATA  POINTS. 

POSSIBLE  VALUES  ARE  64  32.16,8 

ER  A  COUNTER  USED  IN  EXEC  WHICH 

CONTAINS  THE  NUMBER  OF  TIMES 
THE  ERROR  ROUTINE  WAS  ENTERED 

INTCELNUM  AN  INITIAL  VALUE  FOR  "CELNUM" 

WHICH  IS  USED  IN  FOURIE- 
"2ND  PASS  AND  ONWARD"  ROUTINE. 

POSSIBLE  VALUES  ARE  32,16,8,4 

N  THE  NUMBER  OF  DATA  POINTS  TO  BE 

PROCESSED.  USED  IN  EXEC  AND 
FOURIE-  POSSIBLE  VALUES  ARE 
256,128,64,32. 

N2  THE  NUMBER  OF  "COMPLEX"  DATA  POINTS 

USED  IN  FOURIE  POSSIBLE  VALUES 
ARE  128  64.32  16. 

NM  AN  INDEX  USED  IN  FOURIE-  POST 

PROCESSING"  ROUTINE.  IT  IS 
INITIALIZED  TO  THE  RELATIVE 
LOCATION  OF  THE  N-l  COMPLEX 
DATA  POINT.  POSSIBLE  VALUES  ARE 
774  374,174  74. 
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INTRODUCTION 


1 

« 

* 

2 

{2.  INITIALIZATION 

3 

• _ _ _ . _ _ — 

• 

4 

t 

5 

INITAL  MOV 

0TIMEINT.001OO 

6 

MOV 

#0, @0102 

7 

MOV 

#1  .RO 

8 

ASH 

B  RO 

9 

MOV 

R0,N 

10 

• 

11 

MOV 

R0,NM 

12 

SUB 

#2  NM 

13 

ASL 

NM 

14 

15 

ASR 

RO 

16 

MOV 

R0.N2 

17 

* 

18 

ASR 

RO 

19 

MOV 

RO. COUNT 

20 

i 

21 

ASR 

RO 

22 

MOV 

RO, INTCELNUM 

23 

• 

» 

24 

CLR 

ER 

25 

CLR 

RTIME 

26 

CLR 

TIKE 

27 

• 

• 

28 

MOV 

B.BITCOUNT 

29 

SUB 

01 , BITCOUNT 

30 

CLR 

TSUM 

31 

• 

> 

32 

• 

> 

33 

MOV 

Ipwrspc.ro 

34 

35 

MOV 

0360, R1 

36 

37 

CLEAR  1:  CLR 

<R0)+ 

38 

SOB 

Rl, CLEAR 1 

39 

MOV 

0CDFDATA.RO 

40 

MOV 

05500, Rl 

41 

CLEAR2:  CLR 

(R0)+ 

42 

SOB 

R1.CLEAR2 

43 

RTS 

PC 
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PUT  INTERRUPT  ADDR.  IN  LOC  100 
{INTERRUPT  PRIORTTY-0 
;N«2**n 


;NM»2(N-2) 


•N2-N/2 


COUNT «»N  2/2 


; INTCELNUM “COUNT / 2 


•ER»0 

{INITIALIZE  RUN  TIME 
-.AND  THE  TIME  COUNTERS . 

; SET  BITCOUNT-B-1 


;THE  POWER  SPECTRUM 
IS  NOW  CLEARED 
{LOAD  RO  WITH  THE  STARTING 
•AND  R1  WITH  THE  SIZE  OF 
{THE  ARRAY  (127  DECIMAL). 
{CLEAR  THE  ARRAY. 

{REPEAT  FOR  THE  CDF  ARRAY. 


i 


x  .  iimuMHiilUMlilhii 


DEFINITION  OF  DEV 


1 

2 

3 

4 

5 

6 
7 
3 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


.SBm  DEFINITION  OF  DF.V 
•PSECT  RTCRATE ,  RW ,  D , C.BL ,  DEL  OVR 

RTCRATF:  .WORD  0  ;REAL  TIMF.  CLOCK  OVERFLOW 

;RATE  CONTROL.  OVERFLOW  SIGNAL 
RATE  CLOCK  FREQIJENCY/RTC-RATE 

.PSECT  RTCMODE  RW.D .GBL.REL, OVR 

RTCMODE: .WORD  0  ;REAL  TIME  CLOCK  SETTING  FOR 

;THE  RTC  CSR. 

TSUM  .WORD  0 

TIME-  .WORD  0 

N  .WORD  0 

.PSECT  N2,RW,D  GBL.REL, OVR 
N2:  .WORD  0 

NM  .WORD  0 

COUNT:  .WORD  0 

INTCELNUM’  .WORD  0 
.PSECT  ER.RW.D, GBL.REL, OVR 
ER  .WORn  o 

.PSECT  B  RW  D. GBL.REL, OVR 
B:  .WORD  0 

BITCOUNT:  .WORD  0 

.PSECT  I ND  X, RW , 0 , GRL , REL , OVR 
INDX:  .WORD  0 

•PSECT  RTIME  RU,D, GBL.REL, OVR 
RTIME  .WORD  0 

.PSECT  PURS PC  RW.D, GBL.REL. OVR 

PWRSPC:  .REPT  360 
•WORD  0 
•  ENDM 

.PSECT  CDFDATA, RW.D. GBL.REL, OVR 

CDFDATA:  .REPT  5500 
.WORD  0 
.ENDM 
.END 
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START 


GET  (P.-0) 
FOR  SCALING 


ADD  2  TO 
SFACTR 


/isx 

BIT  14 
v  SET  , 


SHIFT -ARRAY 
TWO  PLACES 
TO  THE  RIGHT 


/lS\ 
BIT  13 
'  SET  . 


YES  SET  R4  TO 
H  LOGIC  1 


RETURN 


INCREMENT 

R0 


DECREMENT 

R1 


NO  X  IS 
\R1  •  0 


VN  YES 


/IS  ' 
R4  «  1 


Fig  25.  Flow  Chare  for  FFT  (font)  -  Seale  Routine 
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FPT  — -FFT  SUBROUTINE  MACRO  V03.0I  28-NOV-79  PACE  1 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


.TITLE  FFT  —FFT  SUBROUTINE 

•SBTTL  INTRODUCTION  —INTRO  TO  FFT  AND  ITS  USE. 
.CSECT 


l.  DESCRIPTION: 


THIS  SUBROUTINE  CAN  BE  USED  TO  TRANSFORM  "N"  REAL 
DATA  POINTS.  IT  ASSUMES  THE  DATA  EQUIVALENT  TO  "N2" 
COMPLEX  POINTS  (THE  ODD  REAL  POINTS  ARE  CONSIDERED 
REAL  PARTS  AND  THE  EVEN  PARTS  ARE  CONSIDERED  COMPLEX) 
AFTER  THE  TRANSFORMATION,  A  POST  PROCESSING  OF  DATA  GIVES 
"N2"  FOURIER  TRANSFORMED  FREQUENCY  COMPONENTS  (THE  OTHER 
"N2"  COMPONENTS  ARE  THE  COMPLEX  CONJUGATES  OF  THESE 
VALUES). 


2.  GLOBAL  SYMBOLS: 


THE  FOLLOWING  SY’IBOLS  OF  "FFT”  ARE  INITIALIZED 
OR  USED  IN  OTHER  ROUTINES  THAT  ARE  CALLED  BY  EXEC. 

•GLOBL  FFT  ADDATA , W, N2  NM, COUNT, INTCELNUM, INDX 
.  GLORL  RECOUNT,  SFACTR 


3.  REGISTER  DEFINITION 


.MCALL  -REGDEF 
•REGDEF 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 


4.  DEFINITION  OF  PARAMETERS  USED 


THE  FOLLOWING  ARE  THE  DIFFERENT  PARAMETERS  USED 
IN  THE  FFT  SUBROUTINE  "FFT"  WHICH  ARE  DEFINED 
IN  EXEC. 


N  N  IS  THE  ft  OF  REAL  DATA 

POINTS  IN  OCTAL  FORMAT 

N2:  N2  IS  THE  ASSUMED  ft  OF 

COMPLEX  POINTS  IN  OCTAL. 

B  B  IS  THE  EXPONENT  OF  2 

WHICH  SATISFIES  THE  RELA 
-TIONSHIP  N2=2**B. 

ADDATA  ADDATA  IS  A  POINTER  THAT  WILL 

CONTAIN  THE  STARTING  ADDRESS 
OF  THE  INPUT  (REAL)  POINTS  THAT 
ARE  TO  BE  TRANSFORMED. 

COUNT:  A  COUNTER  SET  TO  HALF  THE 

NUMBER  OF  COMPLEX  POINTS. 


5.  DATA  STORAGE  AREA  ALLOCATION: 


THE  CONSTANT  "ADDATA"  IS  DEFINED  IN  EXEC 
AND  IS  THE  STARTING  ADDRESS  OF  AN  "N"  WORD  ARRAY 
WHICH  CONTAINS  THE  DIGITIZED  DATA  AND  WILL 
CONTAIN  THE  TRANSFORMED  OATA 


FFT  —  DIFFERENT  PARTS  AND  Til  I)  I II  CODE. 


1 

. SBTTL  FFT  --DIFFERENT 

PARTS  AND  THEIR  CODE. 

2 

1.  THE  BIT 

INVERSION  ALGORITHM 

J 

4 

THE  INPUT 

DATA  POINTS  ARE  BIT  INVERTED  (OR  SHUFFLED) 

5 

A 

THAT  THE 

FINAL  TRANSFORMED 

ARRAY  IS  IN  CORRECT  ORDER. 

7 

rFT  CLR 

SFACTR 

;SET  SCALE  FACTOR  TO  0 

a 

CLR 

RO 

RO  WILL  CONTAIN  THE 

9 

; STARTING  ADDRESS  OF  DA 

10 

-TA  TO  BE  TRANSFORMED. 

u 

CLR 

R1 

;R1  WILL  CONTAIN  BIT 

12 

; INVERTED  ADDRESS. 

13 

MOV 

N2  P.5 

; COUNTER  TO  KEEP  TRACK 

14 

•OF  If  OF  DATA  POINTS. 

15 

RIV10:  MOV 

RITC0UNT.R4 

;BIT  SHIFTING  COUNTER 

16 

CLR 

R2 

17 

BIV20:  CLC 

•.CARRY  FLAG  IS  CLEARED 

18 

BECAUSE  BITS  WILL  BE 

19 

•.COPIED  INTO  C-FF  BY 

20 

; ROT AT AT IONS .  WHICH  WILL 

21 

{ACCOMPLISH  SHIFT  FROM  RI 

22 

; INTO  R2  OF  LSB. 

23 

ROU 

R 1 

24 

ROL 

R2 

25 

son 

R4.RIV20 

AT  THIS  POINT  4  LSB  OF 

26 

;R1  ARE  COPIED  TNTO  R2 

27 

IN  BIT  REVERSED  ORDER. 

28 

MOV 

RO,  R  1 

{COPY  FOR  SWAP 

29 

CMP 

R2,  R 1 

{COMPARE  TO  PREVENT 

30 

{DOUBLE  SWAP. 

31 

BLT 

BIV30 

32 

ASL 

R1 

{LEFT  SHTFT  Rl  FOR 

33 

{BYTE  ADDRESSING. 

34 

ASL 

R1 

35 

ASL 

R2 

36 

ASL 

R2 

37 

ADD 

ADDATA,R2 

{STARTING  ADDRESS 

38 

ADD 

ADDATA  R1 

{DESTINATION  ADDRESS 

39 

JSR 

PC  SWAP 

40 

CMP 

(Rl)+. (R2)+ 

: ADVANCE  TO  TMAG  PART 

41 

JSR 

PC  SWAP 

{SWAP  DATA 

42 

BIV30:  INC 

RO 

{GO  TO  NEXT  POINT. 

43 

MOV 

RO,  R 1 

{COPY  R0.R1  FOR 

44 

INVER ION. 

45 

SOB 

R5  RIV10 

46 

NOTE: 

47 

48 

AT  THIS  POINT  THE  BIT  INVERSION  IS  COMPLETE 

49 

AND  DATA  IS  READY  TO  BE 

TRANSFORMED. 
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FFT  —DIFFERENT  PARTS  AND  THEIR  CODE. 


2 

3 

DURING  THIS  PASS  THERE  ARE 

NO  MULTIPLICATIONS,  AND 

4 

THE 

FOLLOWING  EQUATIONS  ARE 

IMPLEMENTED 

J 

6 

R(M) ' 

-R(M)+R(N) 

7 

I(M)' 

-I(M)+I(N) 

8 

R(N) ' 

-R(M)-R(N) 

9 

I  (N) ' 

■I (M)-I (N) 

10 

11 

12 

FPASS : 

MOV 

ADUATA  RO 

; LOAD  STARTING  ADDRESS 

13 

;0F  DATA  FOR  SCALING. 

14 

MOV 

N  Rl 

;L0AD  #  OF  POINTS. 

15 

JSR 

PC , SCALE 

.ALWAYS  SCALE  DATA  TO 

16 

‘.PREVENT  OVERFLOW. 

17 

MOV 

ADOATA  RO 

18 

MOV 

RO,  R I 

;FIRST  POINT  R(M) 

19 

CMP 

<R1)+. (Rl)+ 

2ND  POINT  R(N) 

20 

MOV 

COUNT, R 5 

; COUNTER «N 2/2 

21 

22 

PASS  1 : 

MOV 

(RO)  R2 

;SAVE  R(M) 

23 

ADD 

(R1),(R0)+ 

;R(M)'-R(M)+R(N) 

24 

SUB 

(R1).R2 

;R(N)'»R(M)-R(N) 

25 

MOV 

R2, (Rl)+ 

26 

MOV 

(RO) ,R2 

27 

ADD 

(Rl) > (R0)+ 

;I(M)'-I (M)+I(N) 

28 

SUB 

(Rl ) ,R2 

;I(N)'-I(M)-I(N) 

29 

MOV 

R2, (Rl )+ 

30 

CMP 

(R0)+,  (Rl )+ 

31 

CMP 

(R0)+,  (Rl)+ 

GO  TO  NEXT  PAIR 

32 

SOB 

R5, PASS  1 

;G0  BACK  N2/2  TIMES 

33 

BR 

SPASS 

34 

35 

NOTE 

: 

36 

— 

- 

37 

FIRST 

PASS  IS  OVER  NOW 

GO  TO  SECOND  PASS. 
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FFT  —DIFFERENT  PARTS  AND  THEIR  CODE 


1  ;  3.  FAST  FOURIER  TRANSFORM  (2ND  PASS  AND  ONWARD) : 


> 

3 

9 

THE  2ND  3RD  AND  4TH  PASS  USE  THE  SAME  LOOP, 

4 

t 

HOWEVER  IN  EACH  PASS 

THE  PARAMETERS  DEFINED 

5 

C 

• 

BELOW  TAKE  DIFFERENT 

VALUES : 

7 

> 

DELTAY  0 

;SIZE  OF  INCREMENT  FOR  EXPO 

8 

9 

-NENT  OF  "W"  WHICH  SATISF 

9 

• 

-IES  THE  FOLLOWING  RELATIO 

10 

9 

-NSHIP- 

11 

> 

W«EXP(2*PI*K/N) 

12 

U 

K=0. 1,2 . N-l 

13 

9 

N“  #  OF  POINTS 

14 

9 

15 

CELHUM:  0 

•J  OF  CELLS  IN  EACH  PASS. 

16 

9 

17 

PAIRNM:  0 

OF  PAIRS  IN  EACH  PASS. 

18 

9 

19 

CELDIS :  0 

^DISTANCE  (BYTES)  BETWEEN  A  PAIR  OF 

20 

; POINTS  IN  A  CELL. 

21 

• 

9 

22 

POINT-  0 

•.TEMPORARY  POINTER  TO  STARTING 

23 

; 

ADDRESS  OF  DATA  TO  TRANSFORMED. 

24 

• 

9 

25 

CELCNT '  0 

;CELL  COUNTER  TO  KEEP  TRACK 

26 

9 

OF  #  OF  CELLS  IN  EACH  PASS. 

27 

• 

9 

28 

PARCNT:  0 

; COUNTER  TO  KEEP  TRACK  OF 

29 

9 

#  OF  PAIRS  IN  EACH  CELL. 

30 

9 

31 

SINE:  0 

•TEMPORARY  STORAGE  TO  HOLD 

32 

9 

VALUE  FOR  SIN (Y) . 

33 

• 

9 

34 

cos-  0 

•TEMPORARY  STORAGE  TO  HOLD 

35 

T 

VALUE  FOR  COS (Y) . 

36 

j 

37 

Y  0 

THIS  INDICATOR  HOLDS  THE 

38 

• 

VALUE  OF  "Y"  BEFORE  SINE/ 

39 

9 

COSINE  LOOK  UP  TABLES  ARE 

40 

SEARCHED. 

41 

• 

9 

42 

• 

* 

INTCELNUM : 

THE  INITIAL  VALUE  FOR  "CELNUM" 

43 

• 

t 

WHICH  IS  DEFINED  IN  EXEC- 

PFT  —DIFFERENT  PARTS  AND  THEIR  CODE 


V 

4» 


1 

2 

3 

4 

5 

6 


4.  CODE  FOR  2ND  PASS  AND  ONWARD • 

THE  EQUATIONS  THAT  ARE  IMPLEMENTED  IN  2ND  3RD 
AND  4TH  PASS  ARE  GIVEN  BELOW: 

Rl'»Rl+R2*COS(Y)  >I2*SIN(Y) 


7 

1 1 '  -1 1  -R  2*S  IN  ( Y ) «  2*COS ( Y ) 

8 

R2'«R1-R2*C0S(Y)-I2*SIN(Y) 

9 

I2'«I1+R2*SIN(Y)-I2*C0S(Y) 

10 

11 

WHERE  Y-2*PI*K/N2,  N2-#  OF  COMPLEX  POINTS,  AND 

12 

K-0,1,2..  N2/2-1 

13 

14 

15 

NOTE: 

16 

17 

3RD  AND  4TH  PASSES  START 

FROM  LABEL  "NEWPASS" 

18 

WHEREAS  2ND  PASS  STARTS 

AT  LABEL  "SPASS”. 

19 

20 

SPASS:  MOV  #200, DELTA Y 

■INITIALIZE  DELTAY  TO 

21 

;HALF  THE  SIZE  OF  THE 

22 

; LOOK-UP  TABLE 

23 

24 

MOV  INTCELNUM.CELNUM 

25 

MOV  #2  PAIRNM 

26 

MOV  IIO.CELDIS 

27 

J 

28 

NEWPAS:  MOV  ADDATA  RO 

;GF.T  STARTING  ADDRESS 

29 

AND  #  OF  POINTS  BEFO 

30 

■CALLING  SUB: "SCALE" 

31 

MOV  N.Rl 

32 

JSR  PC , SCALE 

33 

1 

34 

MOV  CELNUM . CELCNT 

•SET  UP  CELL  COUNTER. 

35 

MOV  ADDATA  POINT 

36 

BR  NEWCEL 
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—DIFFERENT  PARTS  AND  THEIR  CC' 


1 

;N0TE:  EACH.  NEW 

CELL  (IN  A 

PASS  STARTS  HERE. 

2 

; - SOME  TEMPORARY  STORAGE 

VARIABLES  ARE  DEFIN 

3 

A 

;  -ED  BELOW: 

* 

5 

I2SINY 

0  ;I2*SIN(Y)  IS  STORED  HERE. 

6 

R2SINY 

0  ;R2*COS (Y)  IS  STORED  HERE. 

7 

I2COSY: 

0  ;I2*COS(Y)  IS  STORED  HERE. 

8 

o 

R2COSY: 

0  ;R2*COS(Y)  IS  STORED  HERE. 

y 

10 

NEWCEL  MOV 

PAIRNM  PARCNT 

It 

CLR 

Y 

CLEAR  Y 

12 

NEWC10:  MOV 

POINT.RO 

13 

MOV 

RO  R 1 

14 

ADD 

C ELD IS  Rl 

15 

MOV 

Y.R2 

;LOAD  REG  2  WITH  THE  VALUE 

16 

TO  GET  THE  SIN  &  COS  FOR. 

17 

MOV 

S INTAB (R2) , SINE  -GET  THE  VALUE  OF  THE 

18 

SINE (R2) 

19 

MOV 

COSTAB (R 2)  COS 

;COS(R2) 

20 

MOV 

(R1).R2 

;CET  REAL  PART  R2 

21 

MUL 

COS.R2 

;R2*COS (Y) . 

22 

ASHC 

)»1,R2 

;LEFT  SHIFT  FOR  FRACTIO 

23 

; -NAL  MULTIPLICATION. 

24 

MOV 

R2,R2COSY 

;R2COSY“R2‘vCOS  (Y) . 

25 

MOV 

(Rl) ,R2 

;GET  REAL  PART  R?. 

26 

MUL 

SINE  R2 

;R2*SIN (Y) 

27 

ASHC 

#1  R2 

28 

MOV 

R2.R2SINY 

;R2SINY«R2*SIN(Y). 

29 

MOV 

R1.R4 

;COPY  PRESENT  ADDRESS- 

30 

TST 

(R4)+ 

;GET  ADDRESS  OF  IMAG. 

31 

MOV 

(R4),R2 

;THIS  IS  VALUE  OF  12. 

32 

MUL 

SINE.R2 

;I2*SIN(Y). 

33 

ASHC 

#1  ,R2 

34 

MOV 

R2  I2SINY 

;I2SINY-I2*SIN(Y). 

35 

MOV 

(R4) ,R2 

36 

MUL 

COS  R2 

;I2*COS(Y) . 

37 

ASHC 

<>1.82 

38  MOV  R2  I2COSY 


I2COSY=I2*COS(Y) 


FFT  —DIFFERENT  PARTS  AND  THEIR  CODE 


1 

NOTE: 

3 

— 

NOW  COMBINE  TERMS  TO  CET 

THE  VALUES  FOR 

4 

< 

R1'.R2\I1*  AUDI 2' . 

6 

MOV 

(RO).  (Rl) 

7 

ADD 

R2C0SY  (RO) 

; (RO)-R1+R2*COS(Y) 

8 

ADD 

I2SINY  (RO) 

;R1'-Rl+R2*C0S(Y)+I2*SIN(Y) 

9 

SUB 

I2SINY. (Rl) 

;  (Rl)™Rl-T2*SIN(Y). 

10 

SUB 

R2C0SY. (Rl) 

{R2'=R1-T2*SIN(Y)-R2*C0S(Y) 

11 

MOV 

R4  Rl 

12 

TST 

(R0)+ 

13 

MOV 

(RO).(Rl) 

14 

SUB 

R2SINY  (RO) 

; (R0)»I1-R2*5TN(Y). 

15 

ADD 

I2C0SY  (RO) 

{ 1 1 '  *»1 1-R2*SIN(Y)+I ?*COS (Y) 

16 

ADD 

R2SINY  (Rl) 

; (RI)«I1+R2*SIN(Y). 

17 

SUB 

I2C0SY  (Rl) 

;I2'“1 1+R2*SI1)Y-I 2*C0S  (Y) . 

18 

19 

20 

NOTE: 

21 

22 

ONE  PAIR 

OF  POtNTS  HAS  BEEN  DONE  HOW  DO 

23 

NEXT  PAIR 

24 

25 

> 

26 

ADD 

#4  POINT 

•NEXT  ADDRESS  IN  THE  CELL. 

27 

ADD 

DELTAY . Y 

; INCREMENT  EXPONENT  OF  W. 

28 

DEC. 

PARCNT 

;ONE  PAIR  DONE- 

29 

BNE 

NEWC.10 

;D0  NEXT  PAIR. EXIT  IF  ALL 

30 

PAIRS  WITHIN  A  CELL  ARE 

31 

{DONE. 

32 

33 

34 

NOTE: 

35 

36 

ONE  CELL 

DONE  NOW  RESET 

EXPONENT  OF  W  AND 

37 

PAIR-COUNTER  (PARC.NT). 

38 

39 

40 

ADD 

C ELD IS  POINT 

{ADDRESS  OF  NEXT  CELL. 

41 

DEC 

CELC.NT 

{ALL  CELLS  DONE? 

42 

BNE 

NEWCELL 

;NO,GO  BACK. 

1 A4 


FFT  —DIFFERENT  PARTS  AND  THEIR  CODE 


C 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


n 


e 
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NOTE: 


ONE  PASS  IS  DONE  NOW  GO  ON  TO  NEXT  PASS, 
SET  NEW  VALUES  OF  PARAMETERS  DEFINED 
IN  PARA  3. 


ASR 

CELNUM 

;LESS  CELLS  THIS  TIME 

BEQ 

POST 

;ALL  DONE  UHF.N  VARIABLE 
;CELNUM«0. 

ASL 

PAIRNM 

•THIS  TIME  THERE  ARE 
;MORE  PAIRS  IN  EACH  CELL. 

ASL 

CELDIS 

•PAIRS  ARE  FURTHER  APART. 

ASR 

DELTAY 

LESS  INC  OF  Y  PER  CELL. 

BR 

NEWPAS 

;G0  BACK  AND  DO  NEXT  PASS 

END  OF  FFT  ALGORITHM. 
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POST  PROCESSItJG  ALGORITHM 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 


•SBTTL  POST  PROCESSING  ALGORITHM 


1.  INTRODUCTION  TO  POST  PROCESSING  ALGORITHM 


POST  PROCESSING  OF  THE  TRANSFORMED  ARRAY  IS  DONE 
BECAUSE  THE  INPUT  (REAL)  DATA  WAS  ASSUMED  TO  BE  COMPLEX 
AND  THE  OUTPUT  OF  THE  FFT  IS  NOT  IN  CORRECT  FORM- 
THE  EQUATIONS  THAT  ARE  IMPLEMENTED  IN  THIS  PASS 
ARE  GIVEN  BELOW: 


AR (M ) - (RP+I P*COS (T ) ~RM*S IN (T ) ) /2 
AI(M)-(IM-IP*SIN(?)-RM*C0S(T))/2 
AR(N-M)-(RP-IP*C0S(T)+RM*SIN(T))/2 
AI (N-M)«(-IM-IP*SIN (T)-RM*COS (T) ) /2 


WHERE  T-PI*M/N,  N-16  AND  M-0,1,2  N/2-1,  AND 

RP-I(M)+I(N-H),  RP-R (M)+R (N-M) ,  IM**I  (M) -I (N-M)  AND 
RM-R  (M)  -R  (N-M) . 


SOME  TEMPORARY  VARIABLES  THAT  ARE  USED  IN 


23 

THIS  PASS  ARE 

DEFINED 

BELOW- 

24 

25 

RP: 

0 

26 

IP: 

0 

27 

RM: 

0 

28 

IM: 

0 

29 

ARM 

0 

;ARM-AR(M) 

30 

AR11M 

0 

;ARNM«AR(N-M) 

31 

AIM 

0 

:AIM-AI(M) 

32 

AINM : 

0 

;AINM«AI(N-M) 

33 

INDX: 

INDEX  TO  GET  SIN, COS  VALUES. 

34 

(INITIALIZED  IN  SETUP) 

35 

NM 

A  RELATIVE  OFFSET  EQUAL  TO 

36 

THE  N-l  COMPLEX  DATA  POINT 

37 

NM-2(N-2)  AND  IS  INITIALIZED 

38 

IN  EXEC. 

POST  PROCESSING  ALGORITHM 


1 

2 
3 


2.  CODE  FOR  POST  PROCESSING  : 


4 

5 

NOTE: 

6 

SET  UP 

ADDRESS  POINTERS  AND  OTHER  RELATED 

7 

8 

;  COUNTERS. 

9 

10 

*OST:  MOV 

ADDATA  RO 

:SET  ADDRESS  POINTER 

11 

;AND  {  OF  POINTS  BEFORE 

12 

; SCALING  DATA  ARRAY. 

13 

MOV 

N,R1 

U 

JSR 

PC, SCALE 

15 

MOV 

ADDATA  RO 

16 

MOV 

COUNT  R1 

;SET  COUNTER  TO  DECIMAL  8. 

17 

MOV 

INDX  R5 

; INCREMENT  FOR  SINK  POINTER. 

18 

MOV 

NM  R4 

.•OFFSET  FOR  ADDRESS  OF 

19 

POINT- 

20 

ADD 

ADDATA. R4 

21 

22 

23 

NOTE: 

24 

25 

THE  POST  PROCESSING 

OF  FIRST  POINT  IS  DONE 

26 

SEPARATELY  BECAUSE  WE  DO  NOT  PHYSICALLY  GE 

27 

-NERATE 

THE  (N+1)ST 

POINT. 

28 

29 

30  POSTlj  MOV 

2(R0)  R2 

;GET  IMAG  PART  OF  1ST 

31 

; POINT. 

32 

ADD 

R2  (RO) 

;ADD  IT  TO  REAL  PART. 

33 

ASL 

(R0)+ 

jDOUBLE  THIS  RESULT  TO 

34 

MAINTAIN  SCALE. 

35 

CLR 

(R0)+ 

•.CLEAR  THE  LOCATION  WH 

36 

37 

38 

39 

40 

NOTE: 

-ERE  THE  IMAG  WAS  STORED. 

41 

NOW  THE 

REST  OF  THE 

POINTS  ARE  PROCESSED. SO 

42 

GET  RM 

IM,RP  AND  IP 

FOR  ALL  EXCEPT  THE  FIRST 

43  ; 

POINT 

U7 


POST  PROCESSING  ALGORITHM 


1 

POSTS:  MOV 

S INTAB (R5) 

,  SINE  SIME=SINE(R5) 

2 

MOV 

COSTAB (R 5) 

COS  ;COS=COS(R5) 

3 

MOV 

(R4) ,RP 

;RP(N-M) 

4 

MOV 

(R4),RM 

;RM=R (N-M) 

5 

NEC 

RM 

;RM**-R(M-M) 

6 

MOV 

2(R4),IP 

;IP=I (N-M) 

7 

MOV 

2<R4),IM 

;IM=I(H-M) 

9 

NEC 

Ill 

;IM**-I  (N-M) 

9 

ADD 

(RO)  RP 

;RP=R(M)+R(N-M) 

10 

ADD 

(RO).RM 

;RM“R (M)-R(N-M) 

11 

ADD 

2(R0)  IP 

;IP“I(M)+I(N-M) 

12 

ADD 

2 (RO) ,  Ill 

; IM=I (M)-I(N-M) 

13 

POST 10:  MOV 

RP  ARM 

14 

MOV 

IM  AIM 

;  ARM-AR(ll)  AND  AIM=AI(M). 

15 

MOV 

RP.ARNM 

16 

17 

MOV 

IM.AINM 

! ARNM =AR (N-M)  AMD 
;AINM-AI(N-M). 

18 

NEC 

AINM 

;AINM—  AI(N-M). 

19 

MOV 

IP  R2 

20 

MUL 

COS.R2 

;IP*COS(T) 

21 

ASHC 

01  R2 

22 

POST20:  ADD 

R2,  ARM 

•ARM«A(M)+IP*COS(T) 

23 

SUB 

R 2. ARNM 

•  ARNM  =*AR  (N-M) -IP*COS  (T) 

24 

POST25:  MOV 

IP,R2 

25 

MUL 

SINE.R2 

;IP*SIN(T) 

26 

ASHC 

#1  ,R2 

27 

SUB 

R2, AIM 

;AIM=AI(tl)-IP*SIN(T) 

28 

SUB 

R2.AINM 

:AIN11“AI  (N-M)  -IP*SIN  (T)  . 

29 

MOV 

RM  R2 

30 

MUL 

SINE,R2 

;RM*SIN (T) 

31 

ASHC 

91  R2 

32 

SUB 

R2. ARM 

; ARM-AR (M )  +1  P*COS  (T )  -RM*S  IN  (T ) 

33 

ADD 

R2.ARNM 

; ARNM=AR  (N-M) -I P*COS  (T )  •. RM*S IN  (T ) 

34 

MOV 

RM,R2 

35 

MUL 

COS.R2 

;RM*COS (T) 

36 

ASHC 

01, R2 

37 

POST50:  SUB 

R2  AIM 

; AIM-AI (M) -I P*S IN (T ) -RM*C0S (T9 

38 

SUB 

R2.AINM 

;AINM=AI (N-M)  -IP*SIN(T) -RM*COS(T) 

39 

POST55:  MOV 

ARM, (RO) 

40 

MOV 

ARNM  (R4) 

41 

MOV 

AIM, 2 (RO) 

42 

MOV 

AINM, 2(R4) 

POST  PROCESSING  ALGORITHM 


i* 


14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


NOTE: 


3  i 

;  ONE  PAIR 

OF  POINTS  HAS 

BEEN  DONE,  NOW 

4 

5 

;  GO  ON  TO 

» 

NEXT  PAIR  OF 

POINTS. 

6 

7 

ADD 

INDX  R5 

; INCREMENT  SINE  INDEX 

3 

CMP 

(R0)+, (R0)+ 

;INC  RO  BY  2 

9 

CMP 

-<R4),-(R4) 

;DEC  R4  BY  2 

10 

DEC 

R1 

;R1<— Rl— 1 

11 

BNE 

POST5 

12 

DEC 

SFACTR 

jDIVIDE  DATA  ARRAY  BY 

13 

RTS 

PC 

END  OF  POST  PROCESSING  ALGORITHM. 
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SUBROUTINE  SCALE 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


•SBTTL  SUBROUTINE  SCALE 


1.  INTRODUCTION  TO  SUBROUTINE  "SCALE": 


THIS  SUBROUTINE  CHICKS  BITS  13-14  AND  SHIFTS  THE 
WHOLE  ARRAY  TO  THE  RIGHT  BY  TWO  BITS  IF  BIT-14  IS  SET, 
AND  SHIFTS  THE  ARRAY  ONE  BIT  TO  THE  RIGHT  IF  BIT- 13 
IS  SET. 

BEFORE  CALLING  THIS  SUBROUTINE  IT  SHOULD  BE 
ENSURED  THAT  THE  STARTING  ADDRESS  OF  THE  ARRAY 
IS  LOADED  IN  kO  AMD  THE  LENGTH  OF  THE  ARRAY  IS  LO 
-ADED  IN  Rl.  THE  SUBROUTINE  IS  CALLED  BY  GIVING  THE 
COMMAND  "JSR  PC  SCALE  . 

THE  VARIABLE  "SFACTR",  WHICH  IS  DEFINED  BELOW, 
KEEPS  TRACK  OF  THE  9  OF  TIMES  THE  ARRAY  HAS  BEEN 
SHIFTED  TO  THE  RIGHT  BY  ONE  BIT. 


. PSECT  SFACTR , RW . D . CBL . REL . OVR 

SFACTR:  0  ; SFACTR  IS  DEFINED  HERE. 

.PSECT 
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SUBROUTINE  SCALE 


1 

o 

2.  CODE  FOR 

SUBROUTINE  SCALE 

3 

4 

SCALE:  MOV 

RO  R5 

;GET  STARTING  ADDRESS 

5 

•AND  STORE  IN  REG  R5. 

6 

MOV 

R1,R3 

jGET  ARRAY  LENGTH  AND 

7 

STORE  IN  REG  R3. 

8 

CLR 

R4 

•CLEAR  R4  MUCH  IS  A  FLAG 

9 

TO  KEEP  TRACK  WHETHER  BIT 

10 

• -14  OR  BIT-13  IS  SET. 

11 

SCAL1:  MOV 

(R0)+,R2 

;GET  EACH  POINT 

12 

BPL 

SCAL2 

;SKIP  NEXT  INSTRUCTION  IF 

13 

; VALUE  IS  POSITIVE. 

14 

NEG 

R2 

;GET  A3S  VALUE 

15 

SCAL2:  BIT 

#40000. R2 

;IS  BIT-14  SET? 

16 

BNE 

RDUCE2 

; YES, SCALE  TWICE. 

17 

BIT 

#20000  R2 

;IS  BIT-13  SET? 

18 

BEQ 

SCAL3 

; NO, GET  ANOTHER  POINT. 

19 

MOV 

#1  R4 

; YES, SET  FLAG  AND  CHECK 

20 

;  11  EXT  POINT. 

21 

SCAL3:  SOB 

R1.SCAL1 

;ALL  POINTS  CHECKED? 

22 

TST 

R4 

;TEST  IF  SCALING  IS  REQUIRED. 

23 

BNE 

RDUCE1 

;YES  IT  IS  REQUIRED. 

24 

RTS  PC 

25 

26 

NOTE: 

27 

28 

SHIFT 

ENTIRE  ARRAY  EITHER  ONE  OR  TOO  BITS 

29 

TO  THE 

RIGHT. 

30 

31 

RDUCE2 :  CLR 

R4 

32 

INC 

SFACTR 

33 

RDUCEls  INC 

SFACTR 

34 

TST 

R4 

35 

BNE 

RDUCE4 

36 

RDUCF.3:  ASR 

<R5) 

SHIFT  EACH  POINT  TWO 

37 

;BITS  TO  THE  RIGHT. 

38 

ASR 

(R5)+ 

39 

SOB 

R3,RDUCE3 

;ALL  POINTS  SCALED? 

40 

RTS 

PC 

;YES 

41 

RDUCE4:  ASR 

(R5)+ 

42 

SOB 

R3  RDUCE4 

43 

RTS 

PC 

44 

45 

46 

47 

48 

49 

END  OF  SCALING  SUBROUTINE. 
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SUBROUTINE  "SWAP" 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 


.SBTTL  SUBROUTINE  "SWAP" 

» 

t 


> 

1.  SUBROUTINE  "SWAP"  (INTRODUCTION  AND  CODE): 


;  THIS  SUBROUTINE  IS  USED  TO  SWAP  THE  CONTENTS  OF 
MEMORY  LOCATIONS  POINTED  TO  BY  REGISTERS  R1  AND  R2. 
;  BEFORE  CALLING  THIS  SUBROUTINE  IT  SHOULD  BE  ENSi'P.ED 
;  THAT  THE  ADDRESS  OF  MEMORY  LOCATIONS  THAT  ARE  TO  BE 
;  SWAPPED  ARE  LOADED  IN  R1  AND  R2. 


SWAP:  MOV 

MOV 
MOV 
RTS 


(R1),R3 
(R2)  (Rl) 
R3,  (R2) 

PC 


END  OF  SUBROUTINE  SWAP. 
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SUBROUTINE  "LOOKUP" 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 


•SBTTL  SUBROUTINE  "LOOKUr" 


1.  INTRODUCTION  TO  "LOOKUP": 


THIS  SUBROUTINE  IS  USED  IN  TWO  PLACES.  "PASS 2  AND 
ONWARD  ,  AMD  "POST  PROCESSING".  THE  LOOK-UP  TABLE  CONTAINS 
THE  SINE  AND  COSINE  VALUES  FOR  K*PHI/N2,  K-l  2...  N2/2. 

THE  TABLE  IS  CONSTRUCTED  FOR  THE  LARGEST  ALLOWED  VALUE  OF 
N2  (128)  WITH  THE  VALUES  FOR  THE  SINES  STARTING  AT  SIULOK, 
AND  FOR  COSINES  STARTING  AT  COSLOK.  EACH  FUNCTION  HAS  128 
ENTRIES  (THE  LAST  HALF  OF  THE  SINK  TABLE  IS  THE  FIRST  HALF 
OF  THE  COSINE  TABLE).  THE  "PARS 2  AND  ONWARD"  SECTION  USES 
EVERY  OTHER  VALUE  (K*PHI/(N2/2))  AND  THE  "POST  PROCESSING" 
SECTION  USES  EVERY  VALUE. 
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SUBROUTINE  "LOOKUP" 


*  * 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
lb 

17 

18 

19 

20 
21 
22 
2  3 

24 

25 
2b 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 
43 


2.  LOOK  UP  TABLE 


SINE(K*PHI/128),  K-l  2,  128  STARTS  HERE 

ft 

S INTAB :  000000 ,001444.003110,004532, 00b 2 1 4 . 00  7b S 3 . 0 1 1  3 1 0 
012742,014371  0 1 SO  1 3, 0 1 74  » ' . 02 1 04 4,02 3 4  SO . 02 404 7 
02543/  027021, 0'>0374. 031  737  03  1272  03461S  036177 
037427.04071b  0431  7  1,0434  IS ,  044064. 04b  1 00, 04  7  100 
0  304 6  4 , 0  S 1  b  3  3 , 0  S  2  7 b  b .  0 34 1 0  l .  0  S  2  0  2 .  0  S  i»  2  b 4 . 0  S  7  3  2  7 
060354,  Obi  3b*  062  ISO.  Ob l 31  / ,  064246.  06M  S3,  06MV,/, 
066712, 067S37, 070341. 0  'H  2S.  07  lbbb,07240S,  0731  O  ' 
073334,  074204 , 0746 1 2 , 07S 175,  0 73‘. 35,  0 7bO‘>2 ,0  7b 344 
076612,0/703S,  0772  13, 07  J4 1 2, 07 /S42, 07  '!>4  7, 07  7  7  30 
077766 

I 

;  COSINE  LOOK -UP  TABLE  STARTS  HERE 

• 

C 0 S TAR :  0  7  7  7  7  7 , 0  7  7  7  6  6 , 0  7  7  7  3 0 , 0  7  7  6  4  7 . 0  7  7  S  4  2 , 0  7  7  4 1  2 

0  7  7  2  3  S ,  0  /  7  0  3  S ,  0  7  6 b  1  * ,  0  7  b  1 4  4 . 0  7  b  0  S  2 . 0  7  3  S  3  S ,  0  7  S  l  7  3 
07461  2 .074204 . 07  ISS4 . 073102  07240S .071 bbb, 0 ' 11 23 
070  34  1 . Ob  7317  Obb  712. 0bb044 ,063133, Ob4 2 4 b  Ob  3 3 1 7 
062  3 SO , Ob  t  3b 2  060334,037327, 0 3 b 2 64,0 S S  2 0 7 , 0 3 4 l n  3 
0S2  766,03 1633,  0SO464  7 300,  04b  1 00 . 044bb4  04  74  IS 
042173  040? lb. 0174 '7  Olbl 2 7, 0  146 1 3. 0  13  3  72  031  7  17 
030374  027021 , 02S437, 0  '404 7, 02 2 4 SO, 02 1 044 , 01 74  12 
016013,014  171  0 12  7  4  2 . 0 1 1  3 1 0 . 0 0  7  6  S  3 , 0 0 6  2 1 4  004332 
003110.001444,000000 

» 

;  Kill'  OF  SINK  LOOK-UP  TABLE 

176334,  1  74670, l 7322S, 171564, 1  7012S  166470, 16S0  16 
163407,  161764.  160346.  136734,  l  SS  3  30,  133731  ,  IS.’  141 
1S07S7,  147404,  146041, 144S06,  14316,3,  141631.  140  131 
137062  135605, 134343,  13  11  14, 1  11  700,  130300. 127314 
126 143,12 SO 12,123673,  1 22S7S, 1 2 1 s 1 4  1204  SO.  1  1  7424 
116416, ll  3427,  1  14460,  1  1  3331,1  1  2622, 111  734,  1  1  1066 
110241,1074  IS, 106633,  106112,  103373  I04h’b,  104224 
103573. 103166, 102603, 102243,101726  101434, 101166 
100743,  1 00  S  4  3 , 10036b,  100236  1001  31, 100047.  100012 
100000 
s 

;  KNn  OF  COSINE  LOOK-UP  TABLE 
.END 
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IV  V02.1-11  Frl  23-Nov-79  12:46  06  PAGE  001 

SUBROUTINE  POWERS PECTRUM 

INTRODUCTION 

THIS  ROUTINE  COMPUTES  THE  POWER  WITHIN  EACH  OF  THE  FREQUENCY 
INTERVALS  IN  THE  FREQUENCY  SPECTRUM.  THE  GLOBAL  ARRAY  ADDATA 
CONTAINS  THE  "COMPLEX"  FREQUENCY  SPECTRUM  AND  THE  ARRAY  "POWER" 
CONTAINS  THE  SUM  OF  THE  SQUARE  OF  THE  REAL  AND  COMPLEX 
COMPONENTS  OF  THE  FREQUENCY  SPECTRUM. 

THE  POWER  SPECTRUM  ARRAY  (POWER)  IS  USED  TO  SUM  TEN  POWER 
SPECTRUMS  TO  FORM  THE  AVERAGE  POWER  SPECTRUM. 

AN  ARITHMETIC  STATEMENT  IS  USED  TO  CONVERT  THE  FIXED  POINT 
VALUES  IN  THE  FREQUENCY  ARRAY  (DATA)  TO  FLOATING  POINT  VALUES. 
THE  SCALE  FACTOR  AND  THE  GAIN  OF  THE  TEMPERATURES '  LINEAR 
AMPLIFIER  (GAIN)  IS  USED  TO  CONVERT  THE  DATA  TO  REAL  FORMAT 


DEFINITION  OF  VARIABLES 


LOCAL  GLOBAL 

ADATA  ADDATA 

GAIN  ADGAIN 

NFREQ  NINTVL 

POWER  PWRSPC 


DEFINITION 

THE  OUTPUT  OF  THE  DFT  MODULE 
CORRESPONDING  TO  THE  FREQUENCY 
SPECTRUM. 

THE  GAIN  OF  TRE  VIBRATION  SIGNAL 
AMPLIFIER. 

THE  HUMBER  OF  FREQUENCY  INTERVALS 
IN  THE  DFT.  POWER  SPECTRUM.  AND 
CDF  OUTPUTS. 

THE  ARRAY  WHICH  IS  USED  TO  FORM 
THE  AVERAGE  POWER  SPECTRUM. 


SCALE  SFACTR  THE  SCALE  FACTOR  FROM  THE  DFT 

ROUTINE.  IT  INDICATES  HOW  MANY 
TIMES  THE  DFT  ARRAY  (ADATA)  HAS 
BEEN  SHIFTED  TO  THE  RIGHT. 


VARIABLE  TYPE  DEFINITION 

INTEGER* 2  ADATA (256), SIZE. K,  I , SCALE , GAIN , NFREQ 
REAL*4  POWER (120) 

COMMON  /PWRSPC /POWER . /ADGAIN/GAIN, /NINTVL /NFREQ 
COMMON  /ADDATA/ADATA,  /SFACTR/SCAI.E 


n 


til  «■ f 


0006 


0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 


C 

C  ARITHMETIC  STATEMENT  FUNCTION  TO  COMPUTE  THE  FLOATING  POINT 

C  VARIABLE  EQUAL  TO  THE  VALUE  OF  THE  FIXED  POINT  FREQUENCY 

C  SPECTRUM  COMPUTED  IN  FOURIER  TIMES  THE  GAIN  OF  THE  LINEAR 

C  AMPLIFIER. 

C 

FPDATA(I)-  FL0AT(ADATA(I))*(2.**SCALE)*FPGAIN 

c 

C  PROGRAM  EXECUTION  STARTS  HERE- 
C 

FPGAIN “FLOAT (GAIN ) 

K-3 

DO  100,1*1. NFREQ 

POWER(I)«POWER(I)+FPDATA(K)**2+FPDATA(K+l)**2 

K-K+2 

100  CONTINUE 
RETURN 
END 
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Fig  27. 


Flow  Chart  for  DAMAGF.FACTOR 
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FORTRAN  IV  V02.1-11  Fri  23-Nov-79  12:47  09 
0001  SUBROUTINE  DAMAGEFACTOR 

C  INTRODUCTION 
C 

C  THIS  MODULE  COMPUTES  THE  CUMULATIVE  DAMAGE  FACTOR  (CDF) 
C  BY  SUMMING  THE  DAMAGE  FACTORS.  THE  DAMAGE  FACTORS  ARE 
C  COMPUTED  USING  AN  ARITHMETIC  STATEMENT  FUNCTION  - 
C  DAMACE  (K, FREQUENCY) WHICH  IMPLEMENTS  EQUATION  15. 

C 

C 

C  VARIABLE  DEFINITION 
C 


C  LOCAL 

p 

GLOBAL 

DEFINITION 

C  CDF 

C 

C 

p 

CDFDATA 

DATA  THE  ARRAY  WHICH  CONTAINS  THE  CUMULATIVE 
DAMAGE  FACTORS.  I  RANGES  FROM  1  TO  127 

AND  TEMP  RANGES  FROM  1  TO  12. 

V» 

C  DELTAF 

p 

DELTA 

DELTA  THE  FREQUENCY  INTERVAL 

C  FREQUENCY 

C 

THE  FREQUENCY  AT  WHICH  THE  CACULATION 

IS  TAKING  PLACE. 

C  NFREQ 

C 

p 

NINTVL 

THE  NUMBER  OF  FREQUENCY  INTERVALS  IN 

THE  POWER  AND  CDF  ARRAYS. 

C  POWER 

C 

P 

PWRSPC 

AN  ARRAY  WHICH  CONTAINS  THE  POWER 

SPECTRUM. 

C  TEMP 

C 

TFLG 

THE  TEMPERATURE  ARRAY  IN  CDF  FOR  WHICH 

THE  DATA  TAKEN. 

C 

C  VARIABLE  TYPE  DEFINITION 
C 

0002  INTEGER* 2  TEMP  NFREQ 

0003  REAL*4  POW£R(120),CDF(120. 12)  FREQUENCY , DELTAF 

C 

C  GLOBAL  VARIABLE  DEFINITION 

C 

0004  COMMON  /PWRSPC/POWER, /CDFDATA/CDF  /TFLG/TEMP 

0005  COMMON  /NINTVL /NFREQ, /DELTA/DELTAFV 

C 

C  ARITHMETIC  STATEMENT  FUNCTION  TO  DETERMINE  THE 

C  DAMAGE  FACTOR 

C 

0006  DAMAGE (K, FREQUENCY) - < (POWER (K) **1.65)/ (FREQUENCY**5. 4) ) 

0007  FREQUENCY “DELTAF 

0008  DO  100,1-1. NFREQ 

0009  CDF ( I , TEMP ) -CD  F ( I . TEMP ) +D AM ACE ( I , FREQUENC Y ) 

0010  FREQUENCY-FREQUENCY+DELTAF 

0011  100  CONTINUE 

0012  RETURN 

0013  END 
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START 


Fig  28.  Flow  Chart  for  XECUTE 
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Fig  28,  Flow  Chart  for  XECUTE  (cont) 


PACE  001 


FORTRAN  IV  V02.1-11  Fri  23-Nov-79  14:57:19 
0001  SUBROUTINE  XECUTE 


C 

C 

C 

C 

C 

C 

C 

C 

c 


IF  THE  TIME  IS  LESS  THAN  OR  EQUAL  TO  ZERO,  A  NEW  TIME  IS 
OBTAINED  BY  CALLING  SUBROUTINE  TIME,  AND  THEN 
THE  DATA  PROCESSING  STARTS. 

VARIABLES  DEFINED  HERE 


0002 

0003 

0004 

0005 


0006 

0008 

0009 


C 

C 

C 


LOGICAL* 1  ANSWER, ERR 
INTEGER* 2  MINUTES 
REAL*4  DURATION 
COMMON  /PTIME/MINUTES 

THE  RUN  TIME  IS  CHECKED  AND  IF  IT  IS  GREATER  THAN 
ZERO  DATA  PROCESSING  IS  STARTED.  IF  NOT  A  NEW  TIME 
IS  REQUESTED. 

IF  (MINUTES  .LT-  0)  CALL  GETT IMF, (DURATION) 

100  TYPE  200, MINUTES/60. 

200  FOR?1AT(2X,  'THE  DURATION  IS  '.F5-2  '  HOURS') 

THE  USER  IS  OFFERED  ONE  LAST  CHANCE  TO  CHANCE  THE  RUN  DURATION. 


C 

C 

C 

C 

C 


0010 

TYPE  225 

0011 

225 

FORMAT (/2X, 'DO  YOU  WANT  TO  CHANCE  IT  ?  Y/N<CR>'$) 

0012 

ACCEPT  250, ANSWER 

0013 

250 

FORMAT (A 1) 

0014 

IF  (ANSWER. EQ  Y')  CO  TO  500 

0016 

IF  ((ANSWER. NE. 'Y')-AND. (ANSWER. NE. 'N'))  GO  TO  600 

0018 

TYPE  300 

0019 

300 

FORMAT  (2X, 'DATA  BEINC  PROCESSED') 

0020 

CALL  EXEC 

0021 

399 

TYPE  400 

0022 

400 

FORMAT  (2X. 'DATA  PROCESSING  COMPLETE'/) 

0023 

RETURN 

0024 

500 

CALL  GETTIME  (DURATION) 

0025 

GO  TO  100 

0026 

600 

CALL  ERROR ('INPUT  MUST  BE  A  Y  OR  A  N') 

0027 

GO  TO  100 

0028 

END 
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POWER-UP 


Fig  29.  Flow  Chart  for  POWERUPDN 
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POWERUPDN  PROCEDURE  MACRO  V03.01  23-MOV-79  15  40:26  PACE  1 


1 

2 

3 

4 

5 

6 

7 
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12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


TITLE  POWERUPDN  PROCEDURE 
SBTTL  INTRODUCTION 


1.  INTRODUCTION- 


THIS  PROCEDURE  IS  USED  WHEN  THE  AUTOMATIC  STARTING  MODE 
IS  SELECTED.  WHEN  SYSTEM  POWER  IS  TURNED  ON,  WITH  THE  RUN/HALT 
SWITCH  IN  THE  RUN  POSITION  (LSI-11  OPERATING  IN  MODE  0,)  PROGRAM 
EXECUTION  STARTS  AT  LABEL  "PURUP".  FIRST,  THE  ADDRESS  OF  THE 
POWER  DOWN  LA3EL  "PWRDN"  IS  STORED  AT  LOCATION  24.  THEN  THE 
GLOBAL  VARIABLE  "PWRFLG"  IS  CHECKED.  IF  IT  IS  EQUAL  TO  EERO 
MEANING  IT  IS  THE  FIRST  POWER-UP,  "PURFI.C"  IS  SET  EQUAL  TO 
ONE  AND  PROGRAM  EXECUTION  JUMPS  TO  THE  LABEL  "EXEC"  IN  THE 
EXECUTIVE  MODULE  IF  "PWRFLG"  DOES  NOT  EQUAL  ZERO,  THE  RUN 
TIME  IS  CHECKED  AND  IF  IT  IS  LESS  THAN  OR  EQUAL  TO  THE  PLANNED 
TIME  PROGRAM  EXECUTION  JUMPS  TO  THE  LABEL  "BEGIN"  IN  THE 
EXECUTIVE  MODULE.  IF  THE  RUN  TIME  IS  EQUAL  TO  OR  GREATER  THAN  THE 
PLANNED  RUN  TIME  PROGRAM  EXECUTION  STOPS. 


2.  GLOBAL  VARIABLES: 


25 

GLOBL 

PWRFLG, 

EXEC, BEGIN, PTIME 

,P.TIME 

26 

MCALL 

•REGDEF 

27 

REGDEF 

28 

29 

30 

PSECT 

PWRFLG 

RW,D,GBL  REL.OVR 

31  PWRFLG: 

.WORD 

0  POWER 

FLAG  USED  TO  DETERMINE  IF 

32 

;FIRST 

POWER-UP. 

33 

PSECT 

34 

35 

36 

3.  CODE  FOR  POWER-UP: 

37 

— 

38 

39  PWRUP 

MOV 

# PWRDN  @#24 

;SET  UP  VECTOR  FOR  THE 

40 

MOV 

#200, @#26 

; POWER-DOWN  ROUTINE. 

41 

TST 

PWRFLG 

;CK  IF  IT  IS  THE  FIRST 

42 

: POWER-UP;  IF  NOT. 

43 

BNE 

RESTRT 

; CHECK  THE  RUN  TIME 

44 

MOV 

#1  PWRFLG 

; PWRFLG  INDICATES  THAT 

45 

;THE  FIRST  POWER-UP 

46 

;HAS  OCCURRED. 

47 

JSR 

PC, EXEC 

{START  DATA  PROCESSING 

48 

;AT  LABEL  "EXEC"  IN  THE 

49 

•.EXECUTIVE  MODULE 
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50 


51 

;  WHEN 

EXECUTIVE  RETURNS 

PROGRAM  CONTROL  TO  THIS 

52 

53 

54 

;  MODULE,  RTIME  WILL  BE  EQUAL 
;  HALT. 

TO  PTIME  SO  PROGRAM  WILL 

55 

RESTRT:  CMP 

PT  IMF,.  RTIME 

CK  RUN  TIME 

56 

BCE 

HLT 

■STOP  IF  RTIME=>PTIME 

57 

58 

59 

60 
61 

JSR 

HLT:  HALT 

1 

PC, BEGIN 

;ELSE  FINISH  PROCESSING 
;THE  DATA. 

62 

;4.  CODE  FOR 

POWER  DOWN  MOD 

63 

64 

65 

66 

PWRDN :  MOV 

#PMRUP(?f?24 

REPLACE  THE  POWER-DOWN 
VECTOR  WITH  THE  POWER- 

67 

68 

HALT 

.END 

;UP  VECTOR  AND  STOP. 
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Appendix  C 
Users'  Manual 


► 


Introduction 

Tie  purpose  of  this  manual  is  to  explain  how  to  use  the  existin'”, 
bench  model  of  the  MFDM.  First,  the  physical  configuration  of  the 
system  is  discussed,  and  then  its  operation  is  discussed.  The  pro¬ 
cedure  for  the  automatic  start-up  mode  is  not  discussed,  since  tills 
mode  of  operation  was  not  incorporated  in  the  bench  model. 

Physical  Configuration 

The  bench  model  system  is  designed  to  be  used  with  a  Lear  Rieglcr 
ADM-3  terminal,  a  Heathkit  H27  dual  floppy  disk  drive  and  a  Meathkit 
H14  printer.  A  block  diagram  for  the  system  is  shown  in  figure  30  with 
the  backplane  configuration  shown  in  figure  31.  The  cable  for  the 
floppy  disk  must  he  connected  to  the  disk  Interface  nodule  with  the  red 
dot  on  the  connector  on  the  top  (visible).  The  bus  expansion  cable 
must  be  connected  between  the  two  backplanes  with  the  M9400-YE  module 
in  the  first  backplane  (with  the  processor  board  in  it). 

The  baud  rate  should  be  set  at  4800  at  the  computer,  the  terminal, 
and  the  printer.  A  slower  baud  rate  can  be  used,  if  desired,  as  long 
'is  all  three  baud  rates  are  the  same.  The  accelerometer  is  connected 
to  the  signal  amplifier,  with  the  output  of  the  amplifier  connected  to 
the  channel  16  leads  of  the  A/n  converter.  Since  the  computer  control 
of  the  auto  ranging  gain  has  not  yet  been  implemented,  the  thumb  wheel 
•witch  must  be  manually  sot  to  a  value  between  0  and  7;  thus,  manually 
••(ting  the  gain  of  the  amplifier.  The  value  to  be  set  depends  upon 
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the  intensity  of  the  input  signal.  In  order  to  determine  what  value  to 
use.  observe  the  output  of  the  amplifier  with  an  oscilloscope  and  set 
the  switch  to  a  position  which  keeps  the  magnitude  of  the  strongest 
signal  below  5  volts.  The  gain  of  the  amplifier,  as  a  function  of  the 
switch  setting,  is  shown  below. 


Switch  Setting 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


Gain  (db) 

-10 

0 

10 

20 

30 

40 

50 

60 

DO  NOT 
USE 


For  the  temperature  input,  the  output  of  the  temperature  amplifier 
is  connected  to  the  channel  5  lends  from  the  A/D  converter  (the 
thermocouple  leads  are  permanently  connected  to  the  anpifier).  This 
completes  the  connections  for  the  vibration  and  temperature  signals 
The  system  i9  now  ready  to  he  turned  on  and  to  be  configured  for  data 
acquisition  and  processing  for  the  required  frequency  range  and 
interval. 

System  Start-Up  Procedure 

First,  the  power  should  be  turned  on  for  all  components  (the 
computer,  the  terminal,  the  printer,  and  the  amplifiers).  After 
application  of  system  power,  the  "RUN/HALT"  switch  should  be  placed  in 
the  "RUN"  position,  and  then  the  "DC  Power"  switch  should  be  turned  on. 
A  "$"  should  appear  on  the  terminal.  If  it  does  not  check  that  all 


units  have  power,  and  that  all  baud  rates  are  sot  to  the  same  value. 
Insert  a  system  disk  with  MFDM.SAV  on  it.  in  the  left  disk  drive 
(drive  0)  and  enter  a  "DX",  followed  by  a  carriage  return,  at  the 
terminal.  If  the  system  disk  has  a  STARTF.COM  file  which  contains  the 
Instruction  "R  MFDM  ,  then  the  system  program  will  start  executin'*, 
automatically.  If  not  the  following  commands  must  be  entered  from  the 
terminal : 

SET  TT : FORM 
SET  TT: WIDTH-96 
R  MFDM 

The  first  two  commands  configure  the  system  to  send  form  feeds  to  the 
terminal,  and  to  set  the  length  of  a  line  to  96  characters.  The  last 
command  starts  the  MFDM  system  program. 

Operating  the  MFDM 

The  R  MFDM  command  starts  the  execution  of  the  MFDM  system 
program.  After  starting,  the  user  is  greeted  with  the  following  lines: 

MICROPROCESSOR  FLIGHT  DAMAGE  MONITOR 

ENTER  S (ETUP,  E(XECUTE,  0(UTPUT,  Q(UIT  -  ?<CR> 

The  first  letter  of  the  function  to  be  performed  is  entered,  followed 
by  a  carriage  return  (<CR>),  and  that  procedure  is  executed.  The 
following  sections  describe  each  of  the  functions  and  their  use3. 

SETUP.  If  an  "S"  is  entered,  the  setup  procedure  is  executed. 

This  procedure  allows  the  user  to  configure  the  system.  The  current 
system  configuration  is  printed  at  the  terminal,  along  with  the 
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following  prompt  line: 


CHANCE  D(ATE,  COMMENTS,  F(REQUEMCY,  R(UN,  DURATION, 

V(ERIFY,  Q(UIT  -  ?<CR> 

The  first  letter  of  the  desired  procedure  should  be  entered, 
followed  by  a  carriage  return,  and  t'.at  procedure  will  be  executed. 

If  a  "V"  is  entered,  the  current  system  configuration  will  be  displayed 
again,  thus  allowing  the  configuration  to  be  checked.  If  a  "Q"  is 
entered,  the  set-tip  procedure  will  terminate,  and  the  MFDM  command  line 
will  be  displayed. 

DATE.  If  the  date  (which  will  be  printed  on  the  cumulative 
damage  matrix  table)  is  to  be  changed,  a  "D<CR>"  is  entered,  and  the 
following  message  is  printed: 

ENTER  HEN  DATE  (13  characters  max.)  -  ?<CR> 

The  new  date  is  then  entered  in  any  form  desired.  For  example 
these  dates  would  be  valid  15-DEC-79,  12-15-79,  or  15  DECEMBER  79. 
However,  the  date  must  be  no  more  than  13  characters  long,  and  must  end 
with  a  carriage  return. 

COMMENTS .  By  entering  a  "C<CR>",  the  comments  can  be  changed. 

The  following  line  is  printed: 

ENTER  COMMENTS  <CR> 

Up  to  40  printable  characters  can  then  be  entered  just  as  the  user 
wants  them  to  appear  on  the  output  table  followed  by  a  carriage 
return.  These  comments  can  be  used  for  identification  purposes  such  as 
the  flight  number,  test  object  and  location. 
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FREQUENCY.  When  an  "F<CR>"  is  entered,  the  following  lines  are 
printed: 


FREQ.  RANGE  (K  Hz)  DELTAF  (Hz)  ENTER 


0  -  6.0 
0  -  6.0 
0  -  3  0 
0-3  0 
0  -  30 
0  -  1.5 
0  -  1-5 
0  -  1.5 


100.2  1 

50. 1  2 

100.2  3 

50.1  4 

25.0  5 

100.2  6 

50.1  7 

25  0  8 


ENTER  THE  NUMBER  CORRESPONDING  TO  THE  DESIRED  COMBINATION. 
ENTER  AN  INTEGER  BETWEEN  1  AND  8  -  ?<CR> 


The  frequency  range  for  which  the  cumulative  damage  factor  (CDF) 
will  be  obtained  is  0  through  the  frequency  selected,  with  the 
corresponding  frequency  interval  -  DELTAF.  The  frequency  range  and 
interval  are  initialized  in  the  system  by  entering  the  integer  (1  to  8) 
corresponding  to  the  desired  frequency  range  and  intervals.  Only  the 
combinations  shown  above  are  available  and  any  entry  other  than  an 
integer  between  one  and  eight  will  result  in  an  error  message  with  a 
request  for  a  new  entry- 

RUN  DURATION .  The  run  duration  is  the  time  for  which  data  will  be 
acquired  and  processed  once  execution  begins  It  can  be  set  by 
entering  an  "R"  followed  by  a  carriage  return,  and  then  the  desired 
duration  should  be  entered.  It  should  be  entered  in  decimal  hours,  and 
followed  by  a  carriage  return. 

Any  of  the  four  items  can  be  changed,  but  all  four  do  not  have  to 
be  changed.  When  all  of  the  desired  changes  have  been  made,  entering  a 
MQ"  followed  by  a  carriage  return  will  terminate  the  set-up  routine, 
and  the  MFDM  command  line  will  be  printed  as  follows: 
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ENTER  S(ETUP,  E(XECUTE,  0(UTPUT,  Q(UIT  -  ?<CR> 

This  completes  the  SETUP  procedure.  The  next  procedure  Is  EXECUTE 
EXECUTE .  The  execute  procedure  is  initialized  from  the  MFDM 
command  line  by  entering  an  "E"  followed  by  a  carriage  return.  The 
following  message  will  then  be  printed: 

THE  RUN  DURATION  WILL  BE  XX  HOURS 
DO  YOU  WANT  TO  CHANGE  IT?  Y/N<CR> 

where  the  XX  represents  the  current  run  time,  which  must  be  greater 
than  zero.  If  the  time  needs  to  be  changed,  this  is  accomplished  by 
entering  a  "Y"  followed  by  a  carriage  return.  The  new  duration  should 
then  be  entered.  If  the  time  does  not  need  changing,  entering  an  "N" 
followed  by  a  carriage  return  will  result  in  the  initiation  of  data 
processing.  The  following  message  will  be  printed  at  the  terminal: 

DATA  BEING  PROCESSED 

Upon  expiration  of  the  run  tine,  data  processing  will  terminate,  and 
the  following  lines  will  be  printed: 

DATA  PROCESSING  COMPLETE 

ENTER  S (ETUP,  E(XECUTE,  0(UTPUT,  Q(UIT  -  ?<CR> 

The  data  is  now  ready  to  be  printed,  and  this  is  accomplished  by  the 
output  routine. 

OUTPUT.  This  routine  prints  out  a  table  of  the  cumulative  damage 
factor  vs  temperature  for  the  frquency  range  selected,  an  example  of 
which  is  shown  in  figure  xx.  The  table  is  designed  to  be  printed  on  a 


96-column  printer*  In  order  to  be  able  to  print  all  of  the  data  in  96 
columns,  the  data  has  been  divided  into  two  categories  depending  on 
temperature*  The  first  set  of  tables  covers  a  temperature  range  of  -50 
degrees  Fahrenheit  to  +100  degrees  Fahrenheit,  in  25  degree  increments, 
printing  41  lines  per  page  until  all  the  frequency  intervals  have  been 
printed.  Then  the  printing  process  is  repeated  for  a  second  temper¬ 
ature  range  of  +75  degrees  Fahrenheit  to  +225  degrees  Fahrenheit. 

After  printing  the  CDF  table,  the  MFDM  command  line  is  displayed, 
and  the  system  is  ready  to  repeat  the  test  run  using  the  same  config¬ 
uration  (enter  an  "E"),  to  be  re-configured  for  a  new  test  (enter  an 
”S"),  or  to  print  another  copy  of  the  CDF  table  (enter  an  "0"). 

Stopping  Program  Execution 

When  all  of  the  testing  has  been  performed,  the  system  should  be 
shut  off.  This  is  accomplished  by  removing  the  disk  switching  the  DC 
power  switch  to  the  position,  and  then  turning  off  all  AC  power  off 

Summary 

This  completes  the  instructions  on  the  operation  of  the  MFDM  bench 
model  system.  The  procedures  to  configure  the  system  and  to  run  the 
programs  have  been  discussed*  If,  for  any  reason,  the  system  does  not 
respond  In  the  correct  way,  remove  the  disk,  shut  down  the  power,  and 
check  the  system  configuration.  Then  re-start  the  MFDM.  If  the  system 
still  fails  to  operate  check  Chapters  V-VII  for  the  operation  of  the 
programs,  in  order  to  determine  the  problems. 
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tion  of  a  microprocessor  based  system  to  acquire  vibration  and 
temperature  information,  and  computes  the  cumulative  damage  factor 
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verr.uc  frequency  ar.d  temperature  in  real-time.  The  cumulative 
damage  factor  is  computed  for  one  of  eight  user  selectable  fre¬ 
quency  ranges  and  for  a  temperature  range  of  -50  to  +225  degrees 
Fahrenheit,  and  is  stored  in  a  non-volatile  memory. 

The  designed  system  used  an  LSI-11  microcomputer  system  con¬ 
sisting  of  an  LSI-11  processor,  an  analog-to-digital  converter,  a 
real-time  clock,  -PROM,  a  non-volatile  memory,  and  a  serial  inter¬ 
face.  The  software  developed  for  this  system  controls  the  data 
acquisition,  the  computation  of  a  fast  Fourier  transform  (FFT), 
the  power  spectrum  of  vibration  data,  and  the  subsequent  computa¬ 
tion  of  the  cumulative  damage  factor.  The  data  acquisition  and 
FFT  modules  were  coded  in  assembly  language, and  the  other  modules 
were  coded  in  FORTRAN.  The  output  of  the  system  consists  of  a 
tabic  containing  the  cumulative  damage  as  a  function  of  frequency 
and  temperature.  A  bench  model  system  was  constructed  and  tested, 
and  a  proposal  for  a  flight-v/orthy  system  is  made. 
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